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SECTION  I 
INTRODUCTION 


Sonic  fatigue  has  become  a  recognizable  and  persistent  structural  design 
problem  over  the  past  twenty-five  years.  Although  not  usually  a 
catastrophic  problem  in  terms  of  human  lives,  it  has  resulted  in 
structural  failures  adversely  affecting  maintenance  costs,  mission 
effectivity  and  often  requiring  major  structural  redesigns.  Sonic  fatigue 
problems  have  been  characterized  by  a  significant  degree  of  inherent 
unpredictability  that  has  so  far  denied  the  structural  designer  the 
precise  analytical  tools  that  are  available  in  other  areas  of  structural 
analysis.  These  limitations  have  been  accompanied  by  the  need  for  minimum 
weight  designs  in  increasingly  severe  and  varied  acoustic  environments. 
This  situation  has  led  to  the  development  and  application  of  semi- 
empirical  design  techniques  based  on  Miles'^'  single  degree-of -freedom 
approach  in  combination  with  experimental  data  from  full-scale  airplane 
tests  and  laboratory  sonic  fatigue  tests.  References  2,  3  and  4  have  used 
such  techniques  to  develop  design  nomographs  for  various  types  of 
structures.  References  5  and  6  present  much  of  this  work  as  part  of 
overall  sonic  fatigue  design  guides. 

These  existing  design  methods  have  been  developed  for  metal  structures. 
However,  recent  advanced  composite  materials  development  has  led  to  a 
wide-spread  aerospace  application  of  nonmetal  lie  structures,  in  the 
interests  of  cost  and/or  weight  savings.  The  most  notable  of  these 
materials  to  date  is  graphite-epoxy.  Although  there  have  been 
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investigations  into  the  sonic  fatigue  resistance  of  graphite 
structures^),  there  are  no  sonic  fatigue  design  methods  available  for 
these  materials  that  are  comparable  to  those  currently  available  for  most 
metal  structures.  Consequently,  it  is  difficult  for  the  designer  to 
translate  the  potential  weight  savings  of  graphite  structures  into  a 
practical  reality  with  the  necessary  level  of  assurance  against  sonic 
fatigue  failures.  The  primary  purpose  of  the  program  described  in  this 
report  was  to  remedy  this  by  developing  a  semi -empirical  sonic  fatigue 
design  method  for  both  flat  and  curved  graphite-epoxy  stiffened-skin 
panels . 

The  program  comprised  three  phases:  analytical,  experimental  and  the 
development  of  a  design  method.  The  analytical  approach  consisted  of 
incorporating  composite  laminate  analyses  into  finite-element  computer 
methods  in  order  to  determine  the  static  and  dynamic  response 
characteristics  of  a  range  of  graphite-epoxy  stiffened-skin  panels.  These 
panels  were  3  x  3,  4  x  3  and  6x3  arrays,  with  various  laminate 
thicknesses,  stiffener  spacings,  radii  of  curvature  and  ply  orientations 
represented.  The  experimental  phase  consisted  of  fabricating  and  sonic 
fatigue  testing  a  range  of  test  panel  configurations  corresponding  to 
those  subjected  to  analysis.  Sonic  fatigue  testing  was  carried  out  in  a 
"progressive-wave  tube"  with  the  panels  being  subjected  to  random  acoustic 
loading  at  grazing  incidence.  The  panels  were  instrumented  with  strain 
gauges  and  flush-mounted  microphor.es,  and  data  taken  over  a  wide  range  of 
sound  pressure  levels.  The  sonic  fatigue  test  program  was  augmented  by 
performing  shaker  tests,  with  random  loading,  on  sections  of  skin- 
laminates  in  order  to  develop  random  fatigue  curves. 

The  design  method  phase  of  the  program  attempts  to  relate  the  analytical 
results  and  the  test  data  in  order  to  provide  a  semi-empirical  design 
method.  Measured  random  strains  are  compared  to  those  calculated  from 
Miles'  equation,  using  the  analytically  determined  static  strains  and 
frequencies  as  inputs.  The  test  results  were  also  compared  to  values 
obtained  from  the  AGARD  nomographs^),  with  density  and  elastic  modulus 
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values  modified  to  reflect  the  graphite-epoxy  laminates.  Finally, 
"multiple  stepwise-regression"  analysis  techniques  were  used  to  develop 
empirical  relationships  between  the  measured  strains  and  frequencies,  and 
various  combinations  of  panel  configuration  parameters  and  finite  element 
analysis  results.  From  these  regression  analyses,  a  set  of  design 
equations  was  developed  and  a  design  nomograph  constructed.  The  design 
method  is  presented  as  a  self-contained  unit  (Section  IV. 5),  allowing  it 
to  be  utilized  independently  of  the  remainder  of  this  report.  A  worked 
example  is  also  presented.  Figure  1  shows  the  program  phase/task  flow 
di aaram. 


Figure  1.  Program  Phase/Task  Flow  Diagram 


SECTION  II 
ANALYTICAL 


1.  INTRODUCTION 

This  section  describes  the  analytical  work  performed  in  support  of  the 
program.  A  description  of  general  sonic  fatigue  theory  is  given  in 
Section  II. 2.  The  analytical  approach  consisted  of  generating  a  complete 
set  of  elastic  properties  for  each  composite  laminate  used  in  the  program; 
these  properties  were  then  used  as  inputs  to  both  the  preliminary  analyses 
and  the  finite-element  solutions.  The  preliminary  analysis  consisted  of 
using  Miles'  equation  and  Reference  5  to  calculate  natural  frequencies  and 
dynamic  stresses  for  each  of  the  proposed  test  panel  configurations.  This 
was  done  in  order  to  ensure  that  their  expected  response  characteristics 
were  compatible  with  the  expected  sonic  fatigue  test  envelope.  The 
finite-element  analysis  consisted  of  constructing  a  series  of  coarse  and 
fine  grid  finite-element  models,  and  using  the  NASTRAN  computer  program  to 
generate  a  set  of  natural  frequencies  and  static  strains  to  be  used  as 
inputs  in  determining  acoustically  induced  dynamic  strains.  An  additional 
set  of  natural  frequencies  was  generated  using  equations  developed  by 
Lin(®) . 

Table  1  lists  the  panel  configurations  used  in  this  program.  An 
analytical  comparison  was  also  made  between  Z  and  J  type  stiffeners. 

Figure  2  presents  a  flow  chart  of  the  analytical  work. 
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Figure  2.  Flow  Chart  of  Analytical  Work 

2.  GENERAL  SONIC  FATIGUE  THEORY 

The  central  problem  in  sonic  fatigue  analysis  is  the  calculation  of 
tne  vibratory  stress  levels  in  structural  panels  subjected  to  the  random 
acoustic  excitation  associated  with  jet  engine  noise,  and  then  to  predict 
the  resulting  fatigue  life.  Since  the  structural  loading  is  random 
(Gaussian),  the  structural  response  is  also  random  and  multimodal  in 
nature.  It  also  follows  that  the  amplitude  distribution  of  the  random 
response  must  be  taken  into  account  in  order  to  determine  corresponding 
fatigue  lives. 

The  complete  response  of  a  complex  structure  to  a  random  noise  field  can 
be  fully  described  by  an  equation  developed  by  Powell However, 
Powell’s  theory  is  too  cumbersome  to  be  used  in  everyday  design  and 
requires  input  data  that  is  never  available  in  the  design  stage  of  a 
vehicle.  Irt  order  to  simplify  the  theory  to  the  level  of  practical  use, 
the  following  assumptions  are  made: 
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(1)  Only  one  mode  of  vibration  contributes  to  fatigue  failure, 
and  that  this  mode  is  the  fundamental  mode  of  the  individual  panel  bays. 
This  mode  is  usually  assumed  to  be  the  fundamental  fully-fixed  mode  or  the 
fundamental  in-phase  mode,  in  which  adjacent  bays  vibrate  in-phase  with 
each  other,  putting  the  panel  stiffeners  into  bending.  Full  scale  tests 
on  aircraft  have  shown  this  assumption  to  be  generally  true. 

(2)  The  vibratory  mode  shape  is  identical  with  the  static 
deflected  shape  of  the  panel  when  subjected  to  a  uniform  static  press:' -e. 

(3)  The  acoustic  pressure  is  exactly  in-phase  over  the  whc : 
panel.  This  assumption  is  reasonable  for  jet  noise  excitation  of 
typically  sized  aircraft  panels.  It  may  not  be  valid  for  boundary- layer 
excitation. 

(4)  The  power  spectral  density  of  the  acoustic  pressure  is 
constant  over  the  frequency  range  near  the  fundamental  natural  frequency 
of  the  panel.  It  is  also  assumed  that  the  whole  of  the  energy  represented 
by  the  acoustic  spectrum  level  at  the  frequency  of  the  assumed  inode  of 
vibration  is  used  to  excite  that  mode. 

These  assumptions  simplify  the  structural  response  equation  to  the  form 
developed  by  Miles' 

Mean  square  stress  o^t)  =  ^n'^^n^o"  (*) 

where  r,  is  the  damping  ratio  of  the  fundamental  mode, 

often  assumed  to  be  typically  0.017  (Reference  b). 

fn  is  the  natural  frequency  of  the  assumed 
fundamental  mode  in  Hz. 

G(fn)  is  the  spectral  density  of  the  acoustic  pressure 
at  the  frequency  fn. 


u0  is  the  static  pressure  at  the  point  of  interest  due 
to  a  unit  uniform  static  pressure  over  the  whole 
or  the  panel . 

This  equation  forms  the  Oasis  of  most  design  oriented  sonic  fatigue  work 
to  date,  including  trie  nomogr  iphs  presented  in  Reference  5.  Many  of  the 
simplified  sonic  fatigue  design  methods  assume  fully-fixed  panel  edges  in 
the  calculation  of  f'n  and  c<0>  In  this  program,  these  values  were  to  be 
determined  from  the  finite-element  solutions,  using  actual  boundary 
conditions. 

The  usual  estimating  procedure,  using  Miles*  equation  is  as  follows: 

a.  Estimate  the  fundamental  natural  frequency  of  the  panel, 
usually  assuming  fixed  edges.  Reference  5  provides  an  appropriate 
nomograph  for  this  purpose. 

b.  Obtain  the  acoustic  spectrum  level  at  the  estimated  frequency. 
NOTE:  The  spectrum  level,  L(fn),  is  the  square  root  of  the 

spectral  density  G(fn).  Since  the  acoustic  spectrum 
level  corresponds  to  the  acoustic  energy  in  a  1-Hz 
bandwidth,  acoustic  data  expressed  in  other  bandwidth  form 
must  be  converted  to  the  spectrum  level  usirg  the  following 
relationship 

L  =  Sound  Pressure  Level  -  10  Log^g  (fg-fi)  (2) 

where  f2  and  fj  are  the  upper  and  lower  frequency  limits,  respectively,  of 
the  given  bandwidth. 

c.  Calculate  cr0.  Reference  10  gives  a  simplified  equation  for 
the  maximum  static  stress  in  a  fully-fixed  panel. 
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d.  Calculate  a(t)  using  Equation  i,  assuming  c=  0 .017 . 

e.  Determine  sonic  fatigue  life  using  specially  generated  random 
fatigue  curves.  Reference  6  contains  examples  of  random  S-N  curves. 

NOTE:  Random  fatigue  curves  can  be  developed  from  conventional  cyclic 

ful ly-reversed  flexural  fatigue  curves.  This  is  accomplished  by 
applying  Miner's^)  cumulative  damage  law  to  the  kayleioh 
distribution  function  for  peak  amplitudes  in  a  Gaussian  process. 

3.  COMPOSITE  LAMINATE  ANALYSIS 

Rohr  has  several  computer  programs  available  tor  analysis  of  composite 
laminates.  These  analytical  techniques  range  from  large  general  purpose 
programs  down  to  simplified  procedures  used  on  the  desk  computers. 

The  primary  general  purpose  program  developed  for  laminate  property 
analysis,  called  COMPOSITE,  calculates  the  laminate  elastic  and  strength 
properties  for  a  specified  laminate  layup.  The  program  may  also  analyze  a 
laminate  with  up  to  five  different  materials  in  the  layup,  hence,  is 
useful  in  determining  the  properties  of  hybrid  laminates.  Options  for 
1  iree  failure  criteria  are  also  included  within  the  program  and  can  be 
used  to  assist  in  determining  laminate  failure  modes.  The  laminate 
analysis  can  be  performed  for  combinations  of  in-plane  and  bending  loads. 

The  COMPOSITE  program  was  developed  with  several  additional  features  for 
the  analysis  of  laminate  properties.  If  the  laminate  fails  under  the 
specified  load,  one  program  feature  will  remove  the  failed  plies  from  the 
layup  and  recalculate  the  laminate  elastic  and  strength  properties.  This 
feature  is  useful  in  evaluating  nonlinearities  due  to  ply  failure  and 
determining  to  what  degree  the  laminate  with  the  failed  plies  removed  can 
sustain  the  load.  The  laminate  stiffness  matrix,  suitable  for  direct 
input  into  the  NASTRAN  finite-element  program,  is  also  computed  and  is 
part  of  the  output. 


The  program  can  also  calculate  the  buckling  coefficients  for  flat  laminate 
panels.  By  selecting  the  options  and  inputting  panel  size-  the  buckling 
coefficient  car.  be  determined  and  displayed  either  in  tabular  or  graphical 
printout.  A  material  data  bank  is  also  incorporated  into  the  COMPOSITE 
program.  Material  laminate  properties  are  stored  within  the  program  and 
may  be  called  by  identification  number  for  a  laminate  analysis.  This 
feature  saves  time  in  setting  up  the  computer  deck  and  provides  consistent 
properties  for  use  on  a  regular  basis.  The  laminate  properties  in  the 
data  bank  may  be  updated  as  necessary  to  reflect  current  data. 

Material  property  data  lor  all  the  skin  and  stiffener  laminates  were 
calculated  and  tabulated  based  on  carpet  plots  in  Reference  12.  Rohr  has 
used  its  "COMPOSITE"  computer  program  and  performed  laminate  property 
tests  to  confirm  selected  data  points  in  Reference  12.  Tables  2  through  6 
list  the  laminate  properties  generated  and  used  as  inputs  for  the  finite 
element  analyses.  Values  shown  in  parentheses  are  computer  generated 
values  used  to  check  those  obtained  from  Reference  12.  Elastic  properties 
were  computed  for  the  skin,  skin/stiffener  attach  flange,  stiffener  web, 
and  the  stiffener  free  flange  with  unidirectional  reinforcement.  Modulus 
values  are  times  106  (lb/in2). 

Effect  of  Stacking  Order 

One  of  the  advantages  of  composite  materials  is  the  capability  to  tailor 
structural  properties  by  dictating  the  number  and  orientation  of  plies. 

T he  in-plane  strength  and  elastic  properties  (Ex,  Ey,  GXy)  of  the  laminate 
can  be  readily  determined  for  specified  orientation  patterns  through  the 
use  of  computer-  programs  or  "carpet  plots."  These  procedures  are 

I  I  o  \ 

documented  in  the  Air  force  Composite  Design  Guide^'  and  other  sources. 

The  elastic  properties  are  customarily  used  in  the  structural  finite- 
element  programs,  such  as  NASTRAN  (see  Section  11.5). 
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Skin  - 
(0,  +  45)s 


Skin  +  attached  stiffener 
flange  - 


St’  ffener,  web  - 


Stiffener,  free-flange  -  10.2 


Modulus  values  are  in  units  of  10”  lb/in.'. 


TABLE  3 

LAMINATE  PROPERTIES  FOR  SHAKER  SPECIMENS  2,  5  AND  6 
AND  PANELS  b,  e,  f,  g,  h,  i,  k  AND  m 


TABLE  4 

LAMINATE  PROPERTIES  FOR  SHAKER  SPECIMEN  3  AND  PANEL  c 


TABLE  6 

LAMINATE  PROPERTIES  FOR  SHAKER  SPECIMEN  7  AND  PANELS  n,  p,  q  AND  ■. 


Laminate 

E,V 

E 

X 

Gyx 

V 

yx 

V 

GTT~ 

E 

y 

G22 

E.v 

G21  =  G12 

V  Ey 

G33 

G 

.yx 

1  -v  V 

XJL.J'X 

1  -  V  V 

xy  yx 

V 

Skin  - 
(0,  90) s 

9.4 

(9.4) 

9.4 

(9.4) 

.65 

(•6b) 

.05 

(.04) 

.05 

9.42 

9.42 

.47 

.65 

(.65) 

Skin  +  attached  stiffener 
flange  - 

4.8 

4.8 

3.25 

.4 

.4 

5,71 

5.71 

2.29 

3,25 

Stiffener,  web  - 

2.4 

(2.3) 

2.4 

(2.3) 

4.5  . 
(4.5) 

.76 

(.76' 

.76 

5.68 

5.68 

4.32 

4.5 

(4.5) 

Stiffener,  free-flange  - 

10.2 

3.0 

2.4 

.61 

.17 

11  .38 

3.35 

1.93 

2.4 

NOTE:  Modulus  values  arc  in  units  of  106  I  b  /  i  n .  2 . 


However,  for  this  investigation  where  composite  structures  are  exposed  to 
a  sonic  environment,  additional  composite  properties  are  desired.  A 
laminated  structure  subjected  to  a  bending  load  whether  applied  by  a  sonic 
or  structural  source  requires  the  use  of  the  inertia  or  bending  stiffness 
properties.  For  laminates,  the  bending  stiffness  is  defined  by  the 
matrix.  The  matrix  is  computed  from  the  individual  ply  properties 
transformed  from  the  specified  orientation  to  the  desired  stiffness 
direction.  The  ply  location  from  the  center  of  the  laminate  is  also  taken 
into  consideration.  The  matrix  is  therefore  written  in  short  notation 
as: 

D.  =  --  1  (0  )  (h3  -  h  3)  (3) 

ij  3  K=1  j '  K  k.  "K-l  1 


The  Qjj  matrix  is  the  in-plane  stiffness  of  each  ply  and  h^,  h^_^  provides 
the  geometric  location.  The  summation  provides  the  bending  stiffness  of 
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the  laminate.  The  position  of  the  ply  in  the  laminate  therefore  will 
dictate  the  stiffness. 

The  effect  of  ply  position  or  "stacking  order"  on  the  laminate  stiffness, 
(D^j),  can  be  determined  by  using  the  COMPOSITE  computer  program.  The 
matrix  results  for  different  ply  stacking  orders  are  shown  in  Table  7.  In 
the  table,  the  stiffness  is  in  the  laminate  0  degree  direction,  the 

d12  is  the  Poisson  effect,  D22  is  the  laminate  90*  stiffness  and  D6b  is 
the  in-plane  shear  stiffness.  Even  the  quasi-isotropic  layup 
(-  45°/90"/0*)2S  has  different  values  in  the  laminate  orthogonal 
directions. 

The  variability  of  the  factors  indicates  that  stacking  order  has  an 
effect  upon  the  performance  of  composite  panels  subjected  to  an  acoustic 
environment.  As  an  example  of  the  stacking  order  effect,  the  natural 
frequency  of  simply  supported  composite  plate  is  of  the  form 

W  »  -~2~  2  K  (4) 

P1 

where  K  =  D„  (f)  +  2  (D1?  +  2  D66)  +  D22  and 

p-j  =  mass  density. 


Since  the  stacking  order  affects  only  the  factor  "K",  its  value  was 
tabul  '.ted  for  various  laminate  layups  in  Table  8.  For  a  sixteen  ply 
lami  ate,  the  stiffness  factor  has  a  12  percent  variation  depending  upon 
the  stacking  order. 


Complications  arose  in  trying  to  quantify  the  effects  of  stacking  order  on 
panel  response.  This  is  discussed  in  Paragraph  II. 5. a. 


TABLE  7 


EFFECTS  OF  STACKING  ORDER  AND  PLY  ORIENTATION  ON 
BENDING  STIFFNESS  MATRIX 

LAMINATE  "D.."  MATRIX 

*  J 

GR/EP  3501/AS  TYPE 

E1  =  17  x  106  E2  =  1.7  x  106  G  =>  .65  x  10S 


Orientation 

No . 

PI  ies 

Dn 

°22 

D12 

. 

°66 

<*4S/902)2s 

16 

246.01 

518.74 

171.32 

187.87 

(M5/02)2s 

16 

518.74 

246.01 

171 .32 

187.87 

(i45/90/0)2s 

16 

351  .7 

413.06 

171.32 

187.87 

(i45/02)4 

16 

549.49 

235.85 

161 .03 

177.58 

(n/90/t45)s 

8 

66.407 

40.48 

5.051 

6.80 

(0/90/  ±45  > 

24 

1793.1 

1093.0 

136.4 

183.6 

(i45)8s 

16 

265.17 

265.17 

202.78 

216.76 

(02/i458)s 

20 

1046.3 

343.28 

219.18 

246.5 

D1 6  =  °26  =  0 


=  7  1  ^ij^K  ^hK  '  hK-T^ 


ij 


K=>1 


Where  (T. .  is  the  transformed  ply  property 

’  J 

and  hK>  is  the  distance  of  the  ply  surfaces 


from  the  reference. 


TABLE  8 

EFFECTS  OF  STACKING  ORDER  ON  NATURAL  FREQUENCY  FACTORS 


Number  of 


Layup  K  Plies 


.4514  16 

.5072  16 


(±45/90/0)-  .4738  16 


.5075  16 


(0/90/±45)s  .149  8 
(0/90/±45)3s  .7745  24 
(±45)8s  .4543  16 


(0-/±45„) 


.6531 
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4.  PRELIMINARY  SONIC  FATIGUE  ANALYSIS 


Preliminary  sonic  fatigue  analysts  were  performed  in  support  of  the 
design  of  the  sonic  fatigue  test  panels,  given  in  Table  1.  These  analyses 
were  made  to  ensure  that  the  application  of  these  panels  in  acoustic 
environments  appropriately  spanned  the  full  range  of  aircraft  application. 
It  was  also  necessary  to  ensure  that  the  sonic  fatigue  resistance  of  the 
test  panels  was  within  the  available  progressive-wave  tube  test  envelope. 
It  is  important  in  a  sonic  fatigue  test  program  to  obtain  a  good  spreaG  cf 
response  characteristics,  and  to  obtain  some  sonic  fatigue  failures  out  in 
the  10®  to  107  cycle  range,  without  having  too  many  panels  fail  either  toe 
quickly  or  not  at  all.  The  AGAKO^®^  sonic  fatigue  desiqn  nomographs  were 
used  in  this  analysis,  with  the  results  being  modified  to  take  account  or 
the  elastic  modulus  and  density  values  for  the  appropriate  skin  laminates. 
The  results  are  shown  in  Table  9.  A  pre-test  evaluation  of  the 
progressive-wave  tube  indicated  that  endurance  testing  would  oe  best 
carried  out  in  the  160  to  165  dB  overall  sound  pressure  level  range, 
corresponding  to  acoustic  spectrum  levels  in  the  130  to  150  dB/Hz  range. 
The  results  show  a  good  spread  in  both  predicted  frequencies  ana  rms 
stresses.  They  also  show  that  almost  all  of  the  panels  could  be  expected 
to  fail  at  an  acoustic  spectrum  level  of  150  dB/Hz. 

5.  FINITE  ELEMENT  SOLUTIONS 

The  general  sonic  fatigue  theory  described  in  Section  11.2  utilizes  as 
inputs  the  static  stresses  or  strains  due  to  a  uniform  unit  pressure  loan 
and  the  natural  frequency  of  the  fundamental  in-phase  stringer-bending 
mode.  These  stress  and  frequency  inputs  were  determined  for  each  of  the 
panel  configurations  given  in  Table  1,  using  a  variety  of  finite-element 
models  in  conjunction  with  the  NASTRAN  computer  program.  NASTkAN  is  a 
general  purpose  finite-element  digital  computer  program  especially  suited 
for  the  analysis  of  large  complex  structures.  Its  ability  to  handle  a 
large  range  of  problems  has  resulted  in  its  adoption  throughout  the 
aerospace  industry.  This  wide  acceptance  and  versatility  are  the  primary 
reasons  for  the  selection  of  NASTRAN  as  the  fundamental  analytical  tool  . n 
this  program.  The  uniform  pressure  load  condition  is  widely  used  in  sonic 


18 


TABLE  9 


PRELIMINARY  SONIC  FATIGUE  A 


Panel 

Config¬ 

uration 


b&e 


Acoustic 

Spectrum 

Level 

dil/Hz 


130 

140 

150 


130 

140 

150 


130 

140 

150 


130 

140 

150 


Stringer 

Spacing 

_  k_(An_L 


Skin 

Thickness 

. _t.  (in)  . 


Radius 

of 

Curvature 

L.liinl. 


.033 


.044 


.044 


.065 


Flat 


Flat 


Flat 


Flat 


130 

140 

150 


.044 


30 


130 

140 

150 


.044 


60 


130 

140 

150 


.044 


90 


130 

140 

150 


.044 


Flat 


130 

140 

150 


.044 


90 


130 

140 

160 


.044 


Hat 


RESULTS 


Fully-  ~~ 

Fixed 

Frequency 

160 

RMS 

Stress 

.  (J  wjjl! _ 

20,900 

66,200 

209,000 

Laminate 

Ftu 

jjb/jfdl 

70,000 

202 

12,800 

40.400 

128,000 

64,000 

247 

14,200 

44,800 

96,000 

142,000 

308 

MOD 

23,400 

74,000 

70,000 

fatigue  work  since  its  resulting  displacement  field  closely  resembles  the 
fundamental  in-phase  mode  shape. 


I 

L 


1 

r 

i 

\ 


»  ft 


The  accuracy  of  finite-element  solutions  is  highly  dependent  upon  the 
element  size  and  the  applied  boundary  conditions.  Consequently, 
considerable  effort  was  put  forth  in  the  determination  of  each.  This 
involved  many  iterations  before  arriving  at  optimum  model  configurations. 
Due  to  the  nature  of  the  test  panel  designs,  i.e.,  relatively  massive 
stiffeners  interfacing  with  thin  plates,  some  difficulty  was  experienced 
in  generating  the  in-phase  mode  from  the  dynamic  models.  Eventually,  well 
defined  stringer-bending  in  phase  modes  were  obtained  for  all  but  two  of 
the  panel  configurations  (the  two  30-inch  curved  panels,  f  and  s,  being 
the  exceptions).  However,  the  mode  shapes  for  the  stiffer  panels 
exhibited  excessive  substructure  deflections,  resulting  in  low  frequency 
estimates.  This  conditioning  problem  was  successfully  overcome  in  the 
static  analysis. 

Finite  element  models  were  also  constructed  for  the  shaker  specimens 
described  in  Section  III. 4.  Computed  natural  frequencies  gave  good 
agreement  with  the  shaker  test  results,  and  are  given  in  Paragraph  II. 5. d. 


a.  Analytical  Approach  -  In  order  to  provide  static  and  dynamic 
analyses  in  sufficient  detail  to  support  the  development  of  a  semi- 
empirical  sonic  fatigue  design  method,  finite-element  models  were 
constructed  to  represent  each  of  the  panel  configurations  shown  in 
Table  1.  Initially  it  was  believed  that  relatively  coarse  grid  models 
would  be  sufficient  for  the  dynamic  analysis.  Consequently,  models 
comprising  2-inch  plate  elements,  with  bar  elements  representing  the 
stringers  and  frames,  were  constructed.  These  models  were  used  to  generate 
a  set  of  static  and  dynamic  solutions.  Although  primarily  intended  as 
dynamic  models,  they  also  provided  a  good  starting  point  for  the  static 
analysis.  The  material  properties  used  in  these  and  subsequent  models 
were  determined  from  the  Kohr  composite  laminate  properties  program 
"COMPOSITE,"  described  in  Section  1 1. 3. 
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The  dynamic  results  from  the  2-inch  grid  models  appeared  to  be 
satisfactory  for  the  3x3  panel  arrays;  however,  it  was  decided  to  use  a 
finer  grid  for  the  panels  with  smaller  bays  (4x3  and  6x3  panel 
arrays).  Models  comprising  1-inch  plate  elements  were  therefore 
constructed  and  a  second  set  of  results  was  generated.  The  mode  shapes 
and  natural  frequencies  generated  by  the  1-inch  and  the  2-inch  models  were 
in  close  agreement.  However,  some  difficulties  were  encountered  with  both 
sets  of  models  in  identifying  the  desired  in-phase  stringer-bending  mode. 
Further  modeling  refinements  did  not  result  in  significant  improvements  in 
the  dynamic  solutions,  consequently  the  1-inch  coarse  grid  quarter  models 
were  used  for  the  dynamic  analyses. 

As  expected,  neither  the  2-inch  nor  the  1-inch  coarse  grid  models  provided 
the  necessary  detail  for  the  static  analysis,  particularly  in  high  stress 
gradient,  areas.  A  set  of  1/2-inch  models  was  constructed  and  another  set 
of  static  solutions  obtained.  Accuracy  at  the  skin-stiffener  interfaces 
was  still  considered  inadequate.  All  of  these  coarse  grid  models 
represented  a  quarter  of  each  panel  array,  as  shown  in  Figure  3. 

It  was  then  decided  to  represent  the  center  bay  portion  of  each  panel  with 
a  fine  grid  model,  also  shown  in  Figure  3.  Because  of  the  lack  of  symmetry 
of  zee  stiffeners,  these  fine-grid  models  included  the  stringers  on  both 
of  the  long  sides..  The  maximum  stress  in  stiffened  skin  panels  occurs  at 
the  center  of,  and  normal  to,  the  longer  edges.  Detailed  accuracy  is 
therefore  of  great  importance  in  these  areas.  This  was  accomplished  by 
representing  the  zee  stiffeners  as  a  series  of  plates  (thereby  creating 
three-dimensional  models)  rather  than  as  simple  bar  elements.  In 
addition,  the  grid  size  was  optimized  by  computing  the  bending  moment 
distributions  from  three  specially  constructed  small  plate  (6  in.  by 
4  in.)  models  employing,  in  turn,  1-inch,  1/2- i rich  and  1/4-inch  grid 
sizes.  The  results  were  then  compared  with  hand  calculations  using 
Timoshenko^).  Figure  4  shows  the  results  of  the  comparison.  It  can  be 
seen  that  1/2-inch  grid  size  provides  accurate  results  at  the  panel 
center.  Even  the  1-inch  grid  has  reasonable  accuracy  at  the  center  of  the 
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Figure  4. 


Comparison  of  Bending  Moment 
Different  Size  Elements  and 
from  Timoshenko  (13) 


Distributions  Using 
PI  ite  Solution  Taken 
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pane!.  However,  the  gradient  at  the  panel  edges  requires  a  very  fine  grid 
in  order  to  achieve  reasonable  accuracy.  A  compromise  between  structural 
accuracy  and  practical  constraints,  such  as  computer  size  limitations, 
resulted  in  1/4-inch  elements  being  selected  to  represent  the  skin- 
stiffener  interface  regions.  This  results  in  an  11-14  percent 
underestimate  in  computed  bending  moments  at  the  panel  edges,  compared  to 
the  values  calculated  using  Timoshenko.  Although  a  maximum  grid  size  of 
0.875-inch  was  used  in  some  non-critical  skin  areas,  all  strains  used  in 
the  development  of  the  design  method  (Section  IV)  were  taken  from  1/4-inch 
elements.  There  are  practical  limitations  in  combining  radically 
different  element  sizes  within  one  model.  In  order  to  limit  the  number  of 
grid  points  to  within  manageable  proportions,  the  smaller  elements 
(1/4-inch)  must  have  higher  aspect  ratios  than  the  larger  elements 
(1/2-inch).  Unfortunately,  analytical  accuracy  deteriorates  with 
increased  aspect  ratio  (above  unity).  An  element  aspect  ratio  of  3:1  is 
considered  the  maximum  for  reasonable  accuracy. 

In  order  to  obtain  boundary  conditions  for  the  fine  grid  center  bay 
models,  a  cubic-spline  computer  program  was  written  to  interpolate  the 
displacement  and  rotation  fields  along  the  interfaces  with  the  2-inch 
coarse  grid  quarter  models.  This  method  assumes  the  deflected  shape  which 
minimizes  potential  (strain)  energy.  Conventional  “beam  theory"  shows 
this  energy  to  be  proportional  to  the  integral,  with  respect  to  the  arc 
length,  of  the  square  of  the  curvature  of  the  spline.  The  accuracy  of 
this  approach  was  verified  using  the  previous  1/2-inch  grid  model. 
Displacement  data  at  2-inch  intervals  on  the  1/2-inch  model  were 
interpolated  to  obtain  intermediate  displacements  at  1/2-inch  intervals. 
These  interpolated  displacements  were  within  1  percent  of  the  actual 
results  from  the  1/2-inch  model.  Initially  it  was  thought  that  the  fine- 
grid  boundary  conditions  could  be  adequately  described  using  displacement 
data  only.  However,  when  this  was  attempted  with  the  1/2-inch  model,  it 
was  found  that  resultant  stress  and  displacement  fields  were  not 
sufficiently  accurate.  Consequently,  it  was  decided  to  also  include  the 
interpolated  values  of  the  two  components  of  rotation  along  the  bounaary. 
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This  improved  I.he  accuracy  of  the  interpolated  stresses  to  within 
2  percent  of  the  1/2-inch  model  results. 

Ouring  the  subsequent  static  analysis  of  the  curved  panels,  it  was 
determined  that  additional  in-plane  displacements  were  needed  in  order  to 
fully  define  the  boundary  conditions.  Flat  panels,  under  normal  loading, 
do  not  undergo  axial  displacements  and  all  the  load  is  taken  in  pure 
bending.  However  curved  panels,  under  normal  loading,  experience  both 
hoop  and  bending  stresses,  requiring  the  application  of  in-plane 
displacement  boundary  conditions.  A  further  refinement  was  evaluated, 
which  was  to  apply  rotation  and  displacement  boundary  conditions  to  the 
out-of-plane  zee  stiffener  elements,  in  addition  to  the  skin  elements. 

The  effects  of  this  refinement  on  one  flat  and  one  curved  panel  was  less 
than  10  percent  and  have  not  been  included  in  the  results  in 
Paragraph  II. 5. c. 

A  comparison  was  made  in  both  the  static  and  dynamic  analyses,  between  Z 
and  J  stiffener  designs.  A  sample  calculation  using  one  of  the  dynamic 
quarter  models  showed  no  significant  differences  in  natural  frequencies 
nor  mode  shapes  between  fhe  Z  and  the  J  stiffened  panels.  No  further 
dynamic  analysis  of  the  »  stiffeners  was  performed.  Significant 
differences  between  the  Z  and  the  J  stiffened  panels  did  occur  in  the 
static  analysis  results,  and  are  given  in  Paragraph  II. 5. c. 

b.  Effects  of  Skin  Ply  Stacking  Order  -  Attempts  were  made  to 
evaluate  the  significance  of  these  effects  on  the  computed  static 
stresses.  The  previous  analyses  utilized  Rohr's  "COMPOSITE"  computer 
program  to  generate  laminate  elastic  properties,  which  are  not  dependent 
upon  stacking  order,  leaving  the  bending  stiffness  (El)  to  be  computed  by 
NASTRAN  in  the  usual  manner. 
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Because  laminated  composite  materials  exhibit  orthotropic  properties,  it 
is  necessary  to  input  the  total  plate  constitutive  equation 


[Nl 

X 


(5) 


in  matrix  form  to  fully  describe  the  behavior  of  a  general  orthotropic 
plate.  Where  I'1  and  M  are  the  applied  loads  and  e  and  k  are  the  resultant 
strains.  The  A  and  0  matrices  define  the  extensional  and  bending 
stiffnesses  respectively.  The  nB"  matrix  defines  the  bending-extensional 
coupling  for  the  laminate.  From  a  practical  standpoint  this  term  is 
nearly  always  zero,  because  ply  orientation  and  stacking  orders  are 
selected  to  give  a  "balanced  symmetric"  layup  which  eliminates  bending- 
extensional  coupling.  The  bending  and  extensional  constitutive  equations 
can  therefore  be  shown  below: 

[N]  =  [A]  [e]  and  [M]  =  [D]  [k]  (6) 

For  the  isotropic  extensional  case,  NASTKAN  computes  the  constitutive 
equations 


A11 

=  a22  =  (— y)  t 

1  -V 

(7) 

A12 

B  vAn 

(8) 

A66 

=  (2"o'+xr)  t 

(9) 

from  the  data  supplied  on  the  PQUAD  quadrilateral  plate  element  card 
(thickness,  t)  and  on  the  MAT1  material  card  (E  and  v) .  The  remaining 
matrix  terms  are  zero. 


If  an  orthotropic  material  is  to  be  analyzed  for  axial  loading,  the  MAT2 
card  is  utilized  to  input  the  material  property  matrix  "G"  terms 


0 -v  V  p  1“  %  ’  etc-)- 

xy  yx;  xy  yx 


(10) 


The  complete  constitutive  equations  are  obtained  by  the  product  of  this 
matrix  and  the  material  thickness  which  is  again  input  on  the  appropriate 
PQUAD  card. 


For  the  isotropic  bending  case,  NASTRAN  computes  the  constitutive 
equations 


11  “22  '  12(l-v2 ) 

(ID 

12  =  vDll  =  U21 

(12) 

=  S  t3 

66  24  0  +V) 

(13) 

from  the  same  data  supplied  for  the  extensional  case.  If  no  additional 
information  is  supplied,  NASTRAN  will  also  compute  the  orthotropic  bending 
constitutive  equations  in  the  same  manner  (multiplying  the  "G"  matrix 
terms  by  the  appropriate  t^/12  term). 

However,  the  true  bending  stiffness  of  an  orthotropic  laminate  is  a 
function  of  the  laminate  stacking  order  in  addition  to  the  laminate 


The  constitutive  equation  for  bending  (the  "D"  matrix)  can  be  input  into 
NASTRAN  using  the  MAT2  card.  In  this  case  the  "D"  matrix  must  be  factored 
by  the  t^/12  term  because  NASTRAN  is  programmed  to  multiply  the  "G"  matrix 
by  the  t^/12  term. 
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The  Cosmic  version  of  NASTRAN  does  not  have  the  capability  to  accept  the 
complete  orthotropic  constitutive  equation  matrix  for  plate  elements. 

In  many  cases,  where  loading  is  primarily  axial  or  the  laminate  has  a 
large  number  of  plies,  the  inaccuracies  introduced  by  using  the 
extensional  "G"  matrix  to  compute  the  bending  constitutive  equations  is 
small  and  this  approach  has  been  used  with  reasonable  accuracy. 

Conversely,  if  the  panel  has  a  small  number  of  plies  and  the  loading  is 
primarily  in  bending,  then  the  bending  "G"  matrix  can  be  used. 

The  Rohr  laminate  analysis  program  COMPOSITE  outputs  the  extensional  "G" 
matrix  directly  in  addition  to  the  "A,"  "B"  and  "0"  matrices. 

In  the  case  of  the  structural  analysis  of  the  sonic  fatigue  panels,  an 
attempt  has  been  made  to  input  both  the  extensional  and  bending 
constitutive  equation  matrices.  This  was  done  using  the  "PQUAL)  1"  general 
quadrilateral  element  property  card  which  is  primarily  utilized  for  the 
analysis  of  sandwich  structures.  This  property  card  allows  for  separate 
input  of  membrane  (extensional),  bending  and  shear  properties. 

The  extensional  information  required  is  the  material  identification  anci 
plate  thickness.  The  extensional  constitutive  matrix  can  be  input  by 
identifying  a  "MAT2"  material  property  card  containing  the  appropriate  "G" 
matrix . 

The  input  data  defining  the  bending  properties  are  the  material 
identification  and  the  area  moment  of  inertia  per  unit  width  (I)  of  the 
quadrilateral  element.  NASTRAN  is  programmed  to  calculate  the  isotropic 
constitutive  bending  equations  using  the  input  values  of  I  (t^/12)„  the 
elastic  constants  and  the  appropriate  numerical  values.  NASTRAN  is  also 
programmed  to  utilize  the  "I"  value  in  the  computation  of  the 
bending  stresses. 


To  obtain  the  orchotropic  constitutive  bending  equations,  the  bending 
material  was  defined  on  a  "MAT2"  material  property  card  containing  the  "D" 
matrix  factored  by  1/1  (12/t^).  The  complete  "D"  matrix  was  then  obtained 
internally  in  NASTRAN  using  the  computational  procedure  defined  above. 

The  results  of  thi:  investigation  were  inconclusive,  requiring  further 
study  beyond  the  schedule  of  this  program.  No  correlation  was  established 
between  the  stacking  order  and  the  effect  on  computed  stresses. 

c.  Static  Analysis  -  The  analysis  of  the  sonic  fatigue  panel 
configurations  in  this  program  required  the  construction  of  detailed 
finite-element  model*'  to  accurately  predict  the  panel  responses  to  a 
uniform  1  lb/in2  applied  pressure  load.  The  iterations  involved  in 
arriving  at  optimum  model  designs  are  discussed  in  Paragraph  II. 5. a.  The 
panel  conf igurations  are  given  in  Table  1,  and  the  location  of  the  finite- 
element  models  relative  to  the  entire  panel  arrays  is  given  in  Figure  2. 

Although  there  are  2C  panel  configurations  in  Table  1,  the  final  analysis 
results  in  this  section  are  limited  to  18  configurations.  Panels  "e"  and 
"m"  were  eliminated  prior  to  the  final  computer  runs.  Because  of  the  many 
iterations  involved  in  obtaining  the  final  analytical  results  and  the 
consequent  effects  that  this  had  on  program  schedule,  it  became  necessary 
to  limit  the  final  analysis  to  those  panels  to  be  used  in  the  development 
of  the  design  method.  Panel  "e"  was  intended  to  evaluate  honeycomb  beam 
stiffeners  and  panel  "m"  was  primarily  intended  as  a  data  link  to  sonic 
fatigue  test  panels,  comprising  a  single  row  of  bays,  typically  used  in 
previous  sonic  fatigue  programs.  This  program  utilizes  panels  comprising 
three  rows  of  bays.  Geometric  similarities  allowed  the  remaining  16 
configurations  to  he  represented  by  10  fine  grid  center  bay  models. 

Figures  A-l  through  A-6  in  Appendix  A  show  six  of  these  models.  The 
three-dimensional  modeling  of  the  stiffeners  is  clearly  seen,  as  is  the 
smaller  qrid  spacing  at  the  panel  edges  and  along  the  center  line  of  the 
center  bay. 
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The  2-inch  coarse  grid  models  were  used  to  generate  boundary  conditions 
for  the  fine  grid  center  bay  model,  utilizing  the  cubic  spline  program 
described  in  Paragraph  II. 5. a.  The  unit  uniform  pressure  load  was  then 
applied  to  the  fine  grid  models,  generating  a  series  of  stress 
distributions.  Although  the  subsequent  design  method  would  be  based  on 
maximum  edge  stresses  and/or  center  bay  stresses,  it  is  desirable  to  know 
the  displacement  patterns  and  stress  distributions  over  the  entire 
surface.  This  was  accomplished  by  plotting  isopleths  of  the  desired 
quantities.  These  are  shown  in  Figures  A- 7  through  A-18  in  Appendix  A. 
Static  deformations,  out-of-plane  displacements  (Z  direction)  and  stresses 
in  the  "y"  and  "x"  directions  are  given  for  panels  “b,"  "d"  and  "f."  The 
plots  are  consistent  with  expected  structural  behavior  and  show  the 
stiffeners  to  provide  good  edge  restraint.  The  stress  contours  show  the 
high  gradients  that  exist  at  the  skin-stringer  interfaces,  demonstrating 
the  need  for  accurate  modeling  in  these  areas.  These  plots  were  generated 
prior  to  some  model  corrections  and  the  introduction  of  in-plane  boundary 
conditions  (see  Paragraph  II. 5. a);  consequently,  the  stress  magnitudes  on 
the  plots  do  not  all  correspond  to  the  tabulated  stresses  shown  later  in 
this  section. 

In  determining  stress  magnitudes  at  critical  locations,  it  was  noticed 
that  the  curved  panels  exhibited  large  stress  differences  on  opposite 
faces  of  the  skin  elements,  indicative  of  significant  axial  stresses. 

This  is  logical  following  the  application  of  the  in-plane  boundary 
conditions  to  the  skin  elements,  described  in  Paragraph  II. 5. a.  Figure  b 
shows  the  locations  of  the  stresses  given  in  Table  ID.  Stresses  on  both 
skin  faces  are  given  for  all  the  curved  panels.  They  are  also  given  for 
the  flat  panels  at  locc:  ions  5  and  6.  Location  A  is  at  the  center  of  the 
center  bay  and  location  6  is  the  maximum  edge  stress.  The  "y"  direction 
is  across  the  narrow  span  and  is  therefore  the  critical  direction.  The 
results  show  the  flat  panels  to  be  in  pure  bending  and  also  show  the 
extent  of  axial  stresses  occurring  on  the  curved  panels.  However,  during 
sonic  fatigue  testing,  back-to-back  strain  gauges  gave  equal  and  opposite 
readings,  indicative  of  pure  bending  on  both  curved  and  flat  panels. 
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Location  of  Points  Used  in  Table  10 


ABLE 


Consequently,  it  was  decided  to  separate  out  the  bending  and  axial  stress 
components  in  the  analytical  results,  and  to  use  only  the  bending  stress 
component  in  the  development  of  the  design  method.  Table  11  lists  the 
results.  Since  the  design  method  will  utilize  strain  values  rather  than 
stress  values,  the  corresponding  bending  strains  are  also  shown  in 
Table  11. 

In  order  to  provide  a  direct  experimental  comparison  to  the  analytical 
static  stresses,  a  static  test  was  performed  on  panel  "d."  The  panel  was 
mounted  in  the  same  fixture  that  was  used  during  sonic  fatigue  testing  and 
that  was  also  used  to  generate  boundary  conditions  for  the  coarse  grid 
models.  A  uniform  pressure  loading  was  incrementally  applied,  from  1  to 
7-lb/in^  using  an  air  bag.  Strains  were  measured  at  each  load  increment 
using  strain  gauges.  Back-to~back  gauges  were  used  at  the  panel  center  to 
measure  a;:ial  strains  in  addition  to  bending  strains.  The  strain  response 
appeared  to  be  nonlinear,  with  stresses  increasing  approximately 
50  percent  for  a  doubling  of  load.  However,  the  back-to-hack  strain 
gauges  gave  readings  within  2  percent  of  each. other,  indicating  pure 
bending.  It  had  been  anticipated  that  any  nonlinearities  in  structural 
response  would  show  up  as  membrane  (axial)  stresses.  No  explanation  is 
offered  for  this  occurrence,  and  no  evidence  of  nonlinear  response 
occurred  during  sonic  fatigue  testing.  Because  of  the  nonlinear  response, 
the  comparison  with  the  analytical  results  varied  depending  upon  which 
load  magnitude  was  used.  At  7  lb/in^  the  center  bay  stresses  were  within 
20  percent  of  the  analytical  value.  The  higher  the  load  the  better  the 
comparison.  Figure  6  compares  the  analytical  results  with  the  one 
1  lb/in^  and  7  lb/in^  test  values.  Although  the  results  from  the  7  lb/in^ 
load  correspond  more  closely  to  the  analytical  results  at  the  panel 
center,  than  do  the  1  lb/in^  results;  the  reverse  is  true  at  the  panel 
edges.  Another  puzzling  aspect  of  this  static  test  was  that  the  biaxial 
strain  relationship  at  the  panel  center  was  markedly  different  during  the 
static  test  from  both  the  analytical  results  and  from  the  sonic  fatigue 
test  results.  In  the  static  test,  the  strains  in  the  long  direction  were 
very  small  (10  percent)  compared  to  the  strains  in  the  short  direction. 
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The  analytical  results  gave  a  corresponding  ratio  of  2:1,  which  is  a  more 
reasonable  ratio.  The  sonic  fatigue  test  results  gave  ratios  a  little  less 
than  2:1.  Thus,  the  static  result  seems  inconsistent  with  both  the 

analytical  and  the  sonic  fat  tout  results,  in  addition  to  appearing  to 

be  less  logical.  Conversely,  the  static  test  results  showed  the  edge 
stresses  to  be  approximately  85  percent  to  90  percent  of  the  center  bay 
stresses,  whereas  the  analytical  results  showed  the  same  edge  stress  to  be 
approximately  three  times  the  corresponding  center  stress. 

Under  fully-fixed  edge  conditions,  the  edge  stress  should  be  twice  the 
center  stress,  thus  the  static  test  result  appears  more  logical.  The 
sonic  fatigue  test  results  showed  a  corresponding  ratio  very  close  to  the 
static  test  results.  It  is  surprising  that  the  analytical  results  would 
produce  a  higher  stress  ratio  between  the  edge  and  center  stresses  than 
one  would  obtain  under  fully  fixed  edge  conditions.  In  summary,  the 
static  test  gave  a  logical  relationship  between  center  and  edge  stresses, 
but  an  unexplained  relationship  between  biaxial  stresses  at  the  panel 
center;  whereas  the  finite-element  results  gave  a  logical  relationship 
between  the  center  panel  biaxial  stresses,  but  a  surprising  relationship 
between  center  and  edge  stresses.  The  sonic  fatigue  test  results,  which 
are  more  tyDically  plagued  with  inconsistencies,  gave  logical 
relationships  for  both  biaxial  and  center-to-edge  stress  ratios. 

A  set  of  analytical  results  was  generated  using  a  J  type  stiffener  in 
place  of  the  Z  stiffeners,  for  stiffener  design  comparison  purposes. 

Since  the  stiffeners  in  this  program  were  adhesively  bonded  to  the  skins, 
the  J  configuration  offers  twice  the  bonded  footprint  area  on  the  skin 
than  does  the  Z  configuration. 

The  static  analysis  utilizing  the  J  stiffeners  was  accomplished  in  much 
the  same  way  as  with  the  Z  stiffeners.  Portions  of  the  previous  finite- 
element  models  were  utilized,  except  for  areas  near  the  stiffeners  which 
were  modified  to  incorporate  the  additional  f'ange  of  the  J  design. 
Identical  boundary  conditions  from  .  .:e  2- inch  coarse  grid  model  were  used 


and  a  unit  pressure  load  was  again  applied.  The  results  for  panel  "b"  are 
shown  in  Figure  7.  As  expected,  the  stress  distributions  for  the  two 
stiffener  designs  are  quite  similar  across  the  majority  of  the  panel,  with 
the  J  stiffener  effecting  a  20  percent  stress  reduction  at  the  panel 
center.  The  major  difference  occurs  at  the  panel  edge,  where  the 
additional  attach  flange  of  the  J  stiffener  significantly  reduces  the  peak 
stress  by  avoiding  the  abrupt  stiffness  change  at  the  attach  radius  of  the 
Z  stiffener. 

d.  Dynamic  Analysis  -  The  dynamic  analysis  of  the  sonic  fatigue 
panel  configurations  in  this  program  required  the  construction  of  finite- 
element  models  to  represent  a  quarter  of  each  panel  array.  The  primary 
purpose  of  the  analysis  was  to  determine  the  natural  frequency  and 
corresponding  mode  shape  of  the  fundamental  in-phase  stringer -bending  mode 
for  each  panel  configuration.  The  quarter  model  (shown  in  Figure  3) 
limits  modal  solutions  to  those  that  are  symmetric  or  antisymmetric  about 
the  panel  array  center  lines,  thereby  excluding  certain  intermediate  inodes 
that  are  not  of  interest  to  this  program.  The  quarter  model  does, 
however,  cover  all  the  bays  in  one  quadrant,  thereby  facilitating  the 
identification  of  a  stringer-bending  mode  in  which  all  bays  vibrate 
in-phase.  Skin  members  were  represented  by  the  NASTRAN  plate  element 
"CQUAD2."  Stringer  and  frame  members  were  inti  ally  represented  by  bar 
elements. 

Problems  were  encountered  in  identifying  an  in-phase  stringer-bending  mode 
for  certain  panel  configurations,  particularly  those  panels  having  greater 
stiffness  due  to  curvature  and/or  close  stringer  spacing.  In  such  cases, 
the  response  was  dominated  by  deflections  of  the  substructure  to  the 
extent  that  computed  frequencies  were  not  responsive  to  changes  in  skin 
thickness.  In  addition,  it  was  not  possible  to  distinguish  between 
overall  panel-array  modes  and  coupled  "bay"  modes.  Such  structural 
behavior  would  be  typical  of  panels  having  inadequate  stiffening  members, 
incapable  of  properly  serving  as  panel  breakers.  However,  both  the  static 
model  results  and  measurements  made  during  sonic  fatigue  testing  showed 
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the  stiffening  members  to  be  adequate  and  to  properly  break  up  the  overall 
panel  arrays  into  their  individual  bays.  This  modal  identification 
problem  was  therefore  assumed  to  be  related  to  the  finite-element  models, 
the  NASTRAN  plate  elements  or  even  f undament.al  analytical  problems 
associated  with  finite-element  techniques.  This  type  of  problem  is  not 
confined  to  this  program.  Previous  sonic  fatigue  programs  have 
reported^2)  similar  difficulties  regarding  the  dominant  behavior  of 
substructure  in  the  dynamic  analysis  of  skin  stringer  structures  using 
finite-element  techniques.  The  problem  is  compounded  in  curved  panels  by 
the  inherent  limitations  of  flat  finite-elements  to  represent  highly 
curved  structures.  NASTRAN  does  have  available  curved  plate  elements. 
However,  advice  from  several  sources,  including  AFFDL,  cautioned  against 
using  them. 

Various  attempts  to  solve  the  problem  were  undertaken.  NASTRAN  has  three 
dynamic  solution  methods  available:  "Inverse-Power,"  "Givens"  and  the 
"Determinant  Method."  The  "Inverse-Power"  method  was  being  used  when  the 
modal  identification  problem  was  encountered.  The  othp^  two  methods  were 
consequently  tried;  however,  the  results  from  all  three  methods  were 
strikingly  similar.  Many  of  the  mode  shapes  obtained  during  this  exercise 
were  observed  to  be  similar  to  those  generated  in  Reference 

Finite-element  methods  and  computer  programs  such  as  NASTRAN  are  known  to 
experience  mass  conditioning  problems  when  analyzing  thin  sheets 
reinforced  with  relatively  massive  stiffeners.  This  aspect  of  the  problem 
led  to  us^'ng  dynamic  models  similar  to  the  coarse  grid  quarter  models  usea 
in  the  static  analysis,  combined  with  representing  the  stiffeners  as  plate 
elements  as  in  the  fine  grid  three-dimeu  ^onal  models  used  in  the  static 
analysis.  This  resulted  in  a  significant  improvement  in  the  generation  of 
the  in-phase  stringer-bending  mode  for  all  but  the  two  30-inch  curved 
panels  (f  and  s).  These  ..vo  panels  failed  to  g\  ,erate  recognizable 
in-phase  modes.  Other  panels  that  had  previously  failed  to  generate  this 
mode  (d,  i,  ,j,  1  and  r)  now  produced  an  in-phase  mode,  but  with  excessive 
stiffener  deflections  and  at  unreasonably  low  frequencies.  The 
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compatibility  of  test  frequencies  with  fully  fixed  frequencies,  calculated 
using  Reference  5,  confirmed  that  the  problem  was  with  the  finite-element 
analysis  results.  Figures  A-19  through  A-24  in  Appendix  A  show  six  of  the 
ten  models. 


A  further  refinement  of  the  model  was  then  made.  Previously,  the 
stiffener  skin  interface  had  been  modeled  witn  a  series  of  single  elements 
whose  properties  were  composed  of  an  homogeneous  superimposition  of  the 
individual  skin  and  stiffener  flange  properties.  It  was  thought  that  the 
interface  between  the  massive  stiffener  element  and  the  thin  skin  element 
could  be  the  source  of  mass  ill-conditioning.  These  areas  were  therefore 
remodeled  with  the  frame  and  skin  elements  individually  represented. 
Connection  between  the  two  was  provided  through  the  use  of  multipoint 
constraints  (MPC)  that  enforce  displacements  of  equal  magnitude,  normal  to 
the  panel,  for  pairs  of  adjacent  grid  points.  This  required  some 
resequencing  of  grid  points,  which  resulted  in  a  significant  increase  in 
the  stiffness  matrix  bandwidth,  The  problem  was  overcome  by  using  a 
preprocessor  program  that  resequenced  the  grid  numbering.  The  results 
from  this  effort  were  disappointing,  however,  with  no  improvement  in  the 
dynamic  response  of  the  problem  panels. 


Another  area  of  concern  in  the  dynamic  analysis  was  the  sensitivity  of  the 
results  to  the  boundary  conditions  applied  to  the  test  panel  fixture 
frame.  During  sonic  fatigue  testing,  both  steel  and  aluminum  frames  were 
used  on  selected  panels.  Also,  changes  were  made  in  the  elastic 
restraining  forces  acting  on  the  panel -fixture  assembly.  Neither  of  these 
variations  influenced  the  dynamic  response  of  the  test  panels.  However, 
the  analysis  results  were  found  to  be  highly  dependent  upon  such 
variations.  It  was  also  noted  that  the  fixture  frame  had  much  greater 
predicted  deflections  from  the  analysis  than  occurred  during  testing.  The 
reasons  for  this  inconsistency  are  not  known.  It  was  decided  to  reduce 
the  influence  of  these  boundary  conditions  and  the  fixture  displacements 
by  modifying  the  finite-element  models  to  eliminate  the  out-of-plane 
motion  of  the  fixture.  This  resulted  in  changes  in  response  frequencies, 
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but  did  not  clarify  the  modal  identification  problems.  The  final  dynamic 
results  were  generated  with  this  fixture  motion  eliminated. 

Although  some  analytical  difficulties  remained  unsolved,  the  majority  of 
the  panels  produced  well  defined  in-phase  stringer-bending  modes,  many  of 
them  occurring  in  the  expected  frequency  range.  Figures  A-25  through  A-28 
in  Appendix  A  show  the  first  four  mode  shapes  for  panel  "b."  The 
frequency  progression  of  these  four  modes  is  interestingly  consistent  with 
the  elastic  properties  of  the  panel  and  its  boundary  conditions.  The  bay 
having  the  maximum  response  is  seen  to  shift  in  turn  from  that  of  least 
fixity  (center  hay),  to  that  bay  with  one  short  side  restrained,  to  that 
with  one  long  side  restrained  and  finally  to  that  bay  with  two  sides 
restrained  (corner  bay),  with  increasing  frequency.  It  is  also  clear  that 
Figure  A-25  is  the  desired  in-phase  stringer-bending  mode,  occurring  at 
171  Hz.  A  list  of  the  complete  set  of  dynamic  solutions  is  given  in 
Tabic  12.  No  solutions  arc  offered  for  panels  "f"  ar,d  "s."  The  first 
four  modes  obtained  for  panel  "f"  are  shown  in  Figures  A-29  through  A-32 
in  Appendix  A.  The  first  modes  are  first  order  modes  within  each  bay,  but 
all  contain  a  combination  of  in-phase  and  out-of ••phase  components.  The 
fourth  mode  (Figure  A-.T2)  shows  the  first  of  the  second-oruer  modes.  One 
last  attempt  was  made  to  force  a  first -order  in-phase  mode  by  forcing  an 
in-phase  displacement,  at  the  center  of  each  bay.  Under  this  condition  the 
moacl  did  not  generate  a  solution  (no  roots  were  found).  This  confirmed 
that  tiie  desired  mode  was  not  simply  being  missed  in  the  modal  search 
procedure,  hut  was  actually  nonexistent  within  the  analytical  framework 
presented  here. 

Early  in  the  program,  finite-element  models  were  constructed  to  represent 
the  shaker  specimens.  The  purpose  of  these  models  was  to  ensure  that  the 
shaker  specimens  were  designed  to  fall  within  the  test  envelope  of  the 
shaker  to  be  used  and  to  avoid  having  shaker  test  specimens  with  torsional 
and  bending  modes  too  close  together.  It  was  important  in  the  shaker  test 
program  to  avoid  exciting  torsional  modes.  This  analysis  also  served  as 
an  indication  of  the  accuracy  with  which  NASTRAN,  combi  nee  with  the 


TABLE  12 

TABULATED  RESULTS  OF  NATURAL  FREQUENCY  SOLUTIONS 


Panel 

Panel 

1st  Bending 
Frame  Free 
(In-Phase) 

Panel 

1st  Bending 
Frame  Fixed 
(In-Phase) 

Panel 

1st  Torsion 
Frame  Free 

Panel 

1st  Torsion 
Frame  Fixed 

Panel  Description 

a 

139 

- 

• 

.. 

3  x  3  Flat 

b 

171 

- 

177 

- 

3x3  Flat 

c 

179 

- 

187 

- 

3x3  Flat 

d 

- 

246 

271 

275 

3x3  Flat 

f 

- 

- 

464 

463 

3  x  3  R  =  30 

g 

398 

- 

285 

- 

3  x  3  R  *  60 

h 

318 

- 

236 

- 

3  x  3  R  s  90 

i 

305*, 332* 

570 

570 

6x3  Flat 

j 

- 

341 

611 

612 

6  x  3  R  =  90 

k 

302 

•  276 

312 

312 

4x3  Flat 

1 

347* 

317* 

520 

520 

4  x  3  R  -  90 

n 

94 

- 

- 

- 

3x3  Flat 

P 

219 

- 

160 

- 

3  x  3  R  -  90 

g 

- 

299 

330 

353 

6x3  Flat 

r 

348* 

- 

483 

6  x  3  R  =  60 

s 

. 

** 

372 

350 

4  x  3  R  =  30 

^Significant  Stringer  Movement 


laminate  properties  generated  by  the  COMPOSITE  program,  would  represent 
the  composite  laminates  used  on  this  program  for  the  sonic  fatigue  test 
panels.  Figure  8  shows  the  finite-element  model  used  to  represent  the 
shaker  specimens.  Table  13  gives  the  first  four  plate  bending  modes  for 
the  shaker  specimens  described  in  Section  III.  Hand  calculated  values 
using  Den  Hartog^4^  are  shown  for  comparison.  The  hand  calculated  values 
assume  the  zee  stiffener  to  represent  a  fully-fixed  support,  resulting  in 
slightly  higher  values  than  the  first  anti-phase  mode.  The  tact  that  the 
first  in-phase  mode  frequencies  are  higher  than  the  hana  calculated  values 
is  probably  due  to  stiffening  effects  of  the  zee  along  its  attach  flange, 
which  effectively  shortens  the  length  of  cantilevered  skin.  Figure  9 
shows  the  first  six  modes  for  shaker  specimen  type  2.  In  the  shaker 
specimen  analysis,  the  tern  "torsion"  refers  to  the  skins  twisting  out-of¬ 
plane  and  not  to  stringer  torsion.  These  results  show  the  torsion  and 
bending  modes  well  separated.  Table  14  shows  a  comparison  of  the  first 
four  skin  binding  mode  frequencies  with  measured  values  on  the  shaker 
table.  The  relatively  close  agreement  between  measured  and  calculated 
values  even  for  the  higher  order  inodes  is  indicative  of  a  sound  analytical 
approach. 

This  early  optimism  turned  out  not  to  be  fully  justified  when  analyzing 
the  more  complex  multi-bay  panels,  as  discussed  earlier  in  this  section. 

This  concluded  the  dynamic  analysis  using  the  f inite-elernerit  models, 
because  of  the  progressive  underestimation  of  computed  frequencies  with 
increasing  panel  stiffness,  these  computed  frequencies  are  thought  to  be 
unsuitable  for  use  in  developing  a  sonic  fatigue  design  method.  In  order 
to  present  alternative  frequency  prediction  techniques,  some  additional 
dynamic  analysis  was  performed  and  the  results  are  presented  in 
Section  II. 6. 
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Figure  8.  Finite  Element  Model  for  Shaker  Specimens 


TABLE  13 

NATURAL  FREQUENCIES  OF  SHAKER  TEST  SPECIMENS 


Panel 

Number 

NASTRAN  Generated 

Frequencies  (Hz) 

Hand  Calculated 
Frequencies 
(Den  Hartog) 

Anti  Phase  Modes 

In  Phase  Modes 

First 

Second 

First 

Second 

1 

40 

254 

61 

346 

49 

2  &  6 

58 

352 

77 

438 

62 

3 

56 

366 

89 

505 

74 

4 

81 

507 

122 

692 

97 

5 

73 

399 

74 

455 

77 

7 

33 

203 

45 

253 

36 
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Panel  Configuration  (Not  to  Scale) 


Mode  Shapes  (Not  to  Scale) 

No.  1  Plate  Bending  (f  =  58  Hz)  -  First  Anti-Symmetric  Mode 


No.  2  Plate  Bending  (f  =  77  Hz)  -  First  Symmetric  Mode 


No.  3  Plate  Torsion  (f  =  220.5  Hz) 


No.  4  Plate  Torsion  (f  =  352  Hz) 


No.  6  Plate  Bending  (  f  =  423.4  Hz)  -  Second  Symmetric  Mode 


Figure  9.  Mode  Shapes  of  Shaker  Specimen  No.  2 


TABLE  14 


COMPARISON  OF  MEASURED  AND  CALCULATED  FREQUENCIES 
FOR  SHAKER  SPECIMEN  TYPE  2 


6.  ADDITIONAL  FREQUENCY  ANALYS2S 

Because  of  the  unresolved  difficulties  with  the  dynamic  analysis  using 
finite-element  techniques,  it  was  decided  to  generate  a  set  of  solutions 
for  the  in-phase  stringer-bending  mode,  using  a  set  of  equations  developed 
by  Lin^).  Lin's  approach  utilizes  differential  equations  applied  to  a 
row  of  continuous  panels.  By  treating  the  skin  and  stringers  as  integral 
parts  of  the  structure,  the  method  utilizes  stringer  properties  in 
addition  to  skin  properties,  thereby  facilitating  accurate  comparisons 
between  different  stringer  properties  and  designs.  Differential  equations 
are  used  to  represent  flat  and  curved  panels.  The  flat  panel  equation  is 
derived  from  a  well  known  fourth  order  equation  of  motion,  applying 
Levy's^*5)  solution  and  appropriate  boundary  conditions  to  develop  the 
following  equation  for  the  frequency  of  the  in-phase  stringer-bending 
mode: 
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where 


and 

A  =  cross-sectional  area  of  stringer 

D  =  bending  stiffness  for  skin,  Lh^/12 ( 1  -v^) 

E  =  modulus  of  elasticity  of  skin  material 

Eb  =  modulus  of  elasticity  of  stringer  material 

In  =  moment  of  inertia  of  stringer  cross  section 

m,  n  =  positive  integers 

w>m  ■  natural  frequency  for  flat  continuous  panels,  radians/sec 

t  1  time,  sec 

b  =  width  of  individual  panel 

h  =  thickness  of  skin 

a  =  length  of  individual  panel 

p  =  mass  density  of  skin  material 

=  mass  density  of  stringer  material 

The  curved  panel  solution  is  obtained  by  expressing  the  strain  and  kinetic 
energies  in  terms  of  generalized  coordinates.  The  equations  of  motion  are 
then  derived  by  a  Lagrangian  formulation. 

The  results  of  this  analysis  are  shown  in  Table  15.  For  convenience, 

i  r-  V 

frequencies  calculated  from  AGARD  design  nomographs 3 1 ,  the  finite-element 
models,  Lin's  equations^8) ,  and  the  sonic  fatigue  tests  are  all  presented 
for  comparison  purposes.  It  is  not  known  why  Lin's  equations  failed  to 
generate  solutions  for  panels  "j"  and  "1."  It  is  Interesting  to  note  that 
Lin's  equations  gave  unexpectedly  high  frequency  values  for  the  30-inch 
curved  panels  (f  and  s);  the  same  two  panels  for  which  the  finite-element 
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TABLE  15 

COMPUTED  AND  MEASURED  NATURAL  FREQUENCIES 


Frequencies 

(Hz) 

Configuration 

Skin 

Thickness 
"t"  (in.) 

Strinq  iii- 
Spacing 
"b"  (in.) 

Radius  of 
Curva  t.ure 
“R"  .(In.)  j 

Agard^ 

Ful iy  Fixed 

Nastran 

,.in<8> 

Test 

d 

.033 

8 

1  '..it 

160 

139 

1  bl 

143 

b 

.044 

a 

flat 

202 

171 

177 

170 

d 

.066 

a 

Hal 

308 

246 

281 

340 

f 

.044 

0 

30 

82b 

- 

1 ,208 

505 

9 

.044 

,1 

60 

bl  0 

398 

583 

350 

h 

.044 

a 

90 

37b 

318 

363 

290 

i 

.044 

4 

i  iat 

6/4 

332 

558 

800 

j 

.044 

4 

90 

1124 

340 

- 

950 

k 

.044 

0 

1  ldl 

312 

30? 

292 

380 

1 

.  Obfc 

6 

90 

b31 

347 

- 

680 

n 

.02? 

1  iat 

120 

94 

101 

140 

p 

.022 

It 

90 

37/ 

219 

447 

180 

q 

.022 

4 

1  l,il 

409 

299 

363 

370 

r 

.033 

4 

60 

782 

348 

442 

780 

s 

.022 

. 

f> 

30 

888 

- 

1,366 

380 

analysis  failed  to  generate  solutions.  It  is  also  interesting  to  note 
that  of  the  three  sets  of  calculated  frequencies,  the  AGARD  nomograph 
results  showed  the  best  correspondence  to  the  test  results. 


SECTION  III 
EXPERIMENTAL 


1.  INTRODUCTION 

The  purpose  of  the  experimental  program  was  to  provide  the  empirical 
data  base  for  the  design  method  described  in  Section  IV.  This  phase  of 
the  program  consisted  of  designing  and  fabricating  a  range  of  shaker  test 
beam  specimens  and  "progressive-wave  tube"  ( F JT )  multi-bay  test  panels. 

The  multi-bay  panels  covered  a  range  of  stringer  spacinqs,  skin  laminate 
thicknesses  and  radii  of  curvature  typical  for  aircraft  appli cation.  The 
configurations  are  shown  in  Table  1.  They  were  instrumented  with  strain 
gauges  and  microphones  and  their  response  characteristics  measured  over  a 
wide  range  of  sound  pressure  levels,  before  being  tested  to  failure.  The 
shaker  tests  augmented  the  PWT  tests  by  providing  additional  random 
fatigue  data  for  the  composite  skin  laminates  used  in  the  multi -bay  panel 
designs. 

2.  TEST  SPECIMEN  AND  FIXTURE  DESIGN 

a.  Progressive-Wave  Tune  Test  Panels  —  Twenty-seven  multi-bay 
test  panels,  comprising  eighteen  configurations,  were  designed  ana 
fabricated  for  subsequent  sonic  fatigue  testing  in  a  progressive-wave 
tube.  Seven  of  the  configurations  had  duplicate  panels  and  one 
configuration  (m)  was  a  reference  panel.  The  reference  panel  provided  a 
data  link  to  a  set  of  existing  Rohr  test  panels,  which  were  also  tested  in 
this  program. 


M 


■  \K, 
1 

V; 


1  . 
r- 


\  > 

'  t, 


4 


These  existing  panels  included  an  aluminum  skin-stringer  panel  to  provide 
a  data  link  between  graphite  and  aluminum  panels.  The  configurations  of 
the  new  panels  are  listed  in  Table  1.  They  comprise  flat,  and  curved 
graphite-epoxy  skins,  ranging  from  4  ply  (.022")  to  12  ply  (.066"), 
stiffened  with  adhesively  bonded  graphite-epoxy  Z  stringers  and  longerons 
in  3  x  3,  4  x  3  and  6x3  panel  arrays.  One  configuration  (e)  had 
honeycomb  beam  stiffeners  and  two  configurations  (b  and  c)  had  additional 
panels  fabricated  with  J  stiffeners,  thus  facilitating  a  comparison 
between  different,  stiffener  designs.  Figure  B-l  in  Appendix  B  shows 
engineering  drawings  of  the  test  panels. 

The  overall  test  panel  size  was  kept  constant  (24-inch  by  36-inch)  in 
order  to  minimize  tooling  and  test;  fixture  costs.  Ten  of  the 
configurations  consisted  of  nine  8-inch  by  12-inch  equal  size  bays  in 
3x3  arrays.  Ten  configurations  were  flat  and  eight  were  curved,  with 
radii  of  curvature  of  30-inch,  60-inch  and  90-inch.  These  curvatures 
encompass  radii  ranging  from  small  aircraft  nacelles  through  to  wide  body 
fuselages.  Five  different  skin  laminates  were  used,  two  of  which 
(b  and  c)  had  the  same  numbe1"  of  plies  but  with  different  orientations. 
This  served  to  isolate  the  effects  of  ply  orientation  for  a  given  number 
of  plies.  With  the  exception  of  the  reference  panel  (m);  4- i rich ,  6-inch 
and  8-inch  stringer'  spacings  were  used  in  6  x  3,  4  x  3  and  3x3  arrays 
respectively.  Panel  e  had  honeycomb  beam  stiffeners,  utilizing  a 
nonmetallic  core  material.  This  is  a  lightweight,  low  cost  stiffening 
concept  whose  sonic  fatigue  resistance  relative  to  the  more  conventional 
Z  stiffeners  is  of  considerable  interest.  Thus,  with  27  panels,  a 
comprehensive  range  of  design  parameters  were  covered,  with  duplicate 
panels  of  some  conf igurations  provided  to  check  test  repeatability  and  to 
provide  more  reliable  fatigue  data  points.  The  panel  parameters  (stringer 
spacing  "b,"  laminate  thickness  "t"  and  radius  of  curvature  "R")  were 
varied  such  that  any  two  panel  responses  can  be  related  hy  varying  one 
parameter  at  a  time,  thereby  facilitating  a  quantitive  identification  of 
the  parametric  cause  of  the  difference  in  response. 
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The  reference  panel  (m)  comprised  a  single  row  of  bays,  duplicating  the 
panel  geometry  of  the  five  existing  Rohr  panels,  shown  in  Figure  8-2  in 
Appendix  B.  Panels  1  and  5  provide  a  direct  comparison  between  aluminum 
and  graphite  multi-bay  structures.  Panels  2  and  3  are  identical 
unstiffeneu  graphite  panels,  representative  of  a  current  nacelle  structure 
in  commercial  service  on  an  experimental  basis.  Panel  4  had  a  single 
honeycomb  stiffener. 

b.  Stiffener  Design  --  Three  types  of  stiffeners  were  evaluated 
in  this  program;  graphite  epoxy  Z  and  J  section  stiffeners  and  honeycomo 
beam  stiffeners  with  graphite  reinforced  caps.  The  program  concentrated 
on  the  Z  stiffeners  which  are  widely  used  on  aircraft  structures.  The 
J  stiffeners  were  included  because  of  their  ability  to  reduce  edge 
stresses  (compared  to  a  Z)  for  minimal  cost  and  weight  increase.  In 
aluminum  structures,  Z  stiffeners  are  inexpensively  formed,  whereas 
0  stiffeners  have  to  be  more  expensively  extruded  or  machined.  In 
graphite  structures,  however,  both  stiffener  types  are  similarly 
fabricated.  Consequently,  the  J  section  is  a  more  cost  effective  design 
in  graphite  than  in  aluminum.  The  honeycomb  stiffeners,  as  mentioned 
earlier,  were  included  for  their  low  cost  and  low  weight  advantages. 

Stiffener  details  are  shown  in  Figure  B-l.  The  Z  and  J  stiffeners  were 
constructed  from  -  45  deg.  graphite  epoxy  laminates  with  unidirectional 
fibers  buried  in  the  free  flanges.  The  number  of  -  45  deg.  laminates  were 
varied  from  configuration  to  configuration  in  order  to  provide  the 
appropriate  stiffnesses  for  their  respective  skins.  The  Z  stringers  were 
1-inch  deep  and  the  longerons  were  1 - 1/2- i rich  deep.  Their  stiffnesses 
were  designed  t.o  ensure  that  they  effectively  served  as  panel  breakers. 
This  was  accomplished  by  using  Reference  8  lu  calculate  their  fundamental 
'n-phase  stringer-bending  mode  frequencies,  and  comparing  the  results  with 
corresponding  fully-fixed  frequencies  of  individual  bays  calculated  using 
Reference  5.  If  the  stiffeners  have  adequate  bending  stiffness,  the 
frequency  of  the  stringer-bending  mode  will  approach  the  fully-fixed 
value.  The  results  of  this  comparison  can  be  seen  in  Table  15,  and  show 
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the  Z  stiffeners  to  be  effective.  At  stringer-longeron  intersections, 
continuity  of  the  attached  and  free  flanges  was  maintained  for  both 
stringers  and  longerons.  The  webs  of  the  stringers  were  also  continuous 
and  were  clipped  to  the  webs  of  the  longerons,  which  were  partially  cut 
away.  Titanium  clips  are  shown  in  Figure  B-l.  These  were  later  changed 
to  steel  following  some  clip  failures  during  the  first  sonic  fatigue 
tests.  For  the  honeycomb  stiffeners,  the  intersections  consisted  of 
honeycomb  core  splices  with  a  foaming  adhesive  locally  applied.  The 
graphite  caps  were  continuous.  At  the  edges  of  the  panels,  the  attach 
flanges  of  the  stiffeners  extended  under  the  test  fixture  frame  and  the 
upstanding  stiffener  webs  were  clipped  to  the  fixture  web.  Figure  10 
shows  a  photograph  of  a  honeycomb  beam  stiffened  pane!.  Examples  of  zee 
stiffened  panels  are  shown  in  Figures  29,  30  and  31.  For  cross  reference 
purposes  with  Section  II  of  this  report,  the  stringers  run  in  the 
X-direction  and  the  longerons  in  the  Y-direction. 


c.  Progressive-Wave  Tube  Fixture  Design  --  The  test  panels  were 
terminated  by  relatively  stiff  channel  sections,  shown  in  Figure  B-l. 

There  were  two  fixturing  approaches  considered  in  this  program.  One  was 
to  bolt  a  picture  frame/panel  assembly  rigidly  to  the  progressive-wave 
tube  (PWT) .  The  other  was  t.o  suspend  a  stiff  picture  frame/panel  assembly 
on  captive  wires.  The  former  approach  more  closely  approximates  fully- 
fixed  edge  conditions,  which  is  more  convenient  when  the  test  results  are 
to  be  compared  to  simplified  analyses  where  fixed  edges  are  usually 
assumed.  This  approach  has  the  disadvantage  of  the  panel  response  being 
affected  by  vibrations  in  the  PWT  itself.  The  latter  approach  eliminates 
any  response  interference  from  the  PWT  and  allows  the  fixture  to  be 
accurately  represented  in  the  finite-element  models.  Since  the  fixture  is 
constrained  only  spatially  in  the  PWT,  its  boundary  conditions  are  able  to 
be  accurately  represented  in  the  analysis.  The  latter  approach  was  chosen 
for  this  program.  An  additional  advantage  of  this  approach  was  that  it 
was  a  relatively  simple  matter  to  remove  a  panel  from  its  frame  and 
thereby  make  response  measurements  under  two  different  boundary 
conditions,  in  order  to  evaluate  their  effects.  A  set  of  steel  frames  was 


Figure  10.  Honeycomb  Beam  Stiffened  Panel  (e) 


added  to  the  program  as  back-up  f ixtur ing  should  the  aluminum  frames 
experience  any  sonic  fatigue  damage.  This  did,  in  fact,  occur  on  one  of 
the  curved  panels,  and  the  steel  frames  were  used  for  the  remaining  curved 
panels.  Some  panel  response  data  were  taken  in  both  tne  steel  and 
aluminum  frames  for  comparison  purposes.  The  test  fixturing  arrangement 
is  discussed  further  in  Section  III. 5,  and  can  be  seen  in  Figures  22  and 
24. 


d.  Shaker-Test  Specimens  --  Eighty-one  shaker  test  specimens  were 
designed  and  fabricated.  Figure  B-3  in  Appendix  B  shows  the  specimen 
details.  The  specimens  consisted  of  3-inch  by  10-inch  sections  of  skin 
laminates,  with  stringer  sections  attached  across  the  short  dimension. 

The  specimens  represented  each  skin-stiffener  combination  used  on  the 
multi-bay  test  panels.  A  set  of  specimens  having  the  stringer  riveted  to 
the  skin  was  included  to  provide  for  a  fatigue  life  comparison  with  bonded 
joints.  The  shaker  specimens  were  intended  for  ratigue  testing  the  skin- 
stringer  joints  in  order  to  develop  fatigue  curves  to  augment  the 
progressive-wave  tube  test  results.  This  objective  was  not  fully  realized 
due  to  some  adhesive  bonding  quality  problems  encountered  early  in  the 
program.  These  problems  and  their  effects  are  discussed  in  Sections  III. 3 
and  III. 4. 

The  fixturing  for  the  shaker  tests  was  originally  a  simple  tee  section, 
15-in.  long,  accomodating  five  specimens  at  a  time.  The  upstanding  webs 
of  the  Z  stringers  were  mechanically  fastened  to  the  upstanding  leg  to  the 
tee.  The  assembly  was  then  simply  bolted  to  the  shaker  table  for  testing. 
This  fixturing  was  later  changed  when  the  shaker  test  program  was  modified 
as  a  result  of  the  adhesive  bonding  problems  mentioned  above.  The  changes 
are  discussed  in  Section  III. 4. 

3.  TEST  SPECIMEN  FABRICATION 

Fabrication  of  the  sonic  fatigue  test  panels  involved  the  manufacture 
and  assembly  of  graphite  epoxy  skins  ana  stiffeners  fabricated  from 
Hercules  AS-3501  Pre-Preg.  The  skins  were  laid  up  and  cured,  on  a  flat  or 
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curved  tool  as  appropriate,  from  12-inch  wide  graphite  tape.  Cach  ply  was 
oriented  with  respect  to  a  reference  direction  in  order  to  build  up  the 
desired  panel  stiffness  properties.  After  cure  of  the  skin,  the 
stiffeners  were  attached  in  a  secondary  bond  cycle  using  3M's  AF147 
adhesive . 

Tii  Z  and  J  stiffeners  were  laid  up  and  cured  on  a  separate  tool.  Layup 
of  these  stiffeners  included  -  45  deg.  plies,  from  the  flanges  through  the 
web,  for  shear  stiffness  and  strength.  Additional  unidirectional  fibers 
were  added  to  the  fret  flanges  of  the  stiffeners  for  bending  stiffness  and 
strength.  Following  adhesive  bonding  assembly  with  the  skin,  the 
stringer-longeron  intersections  were  stabilized  with  angle  clips,  which 
were  mechanically  fastened  in  place. 

For  the  sandwich  stiffened  panels,  the  stiffeners  were  fabricated  in  place 
on  the  skins  by  cocuring  the  core-to-skin  ami  cap-to-core  bonds 
simultaneously  with  the  cure  of  the  cap.  Since  the  longerons  utilized 
deeper  core  than  did  the  stringers,  the  caps  for  both  were  continuous 
across  the  intersections.  The  only  tie  then  required  at  the  intersections 
was  a  foaming  core  splice  adhesive,  cocured  with  the  remainder  of  the 
stiffeners  as  described  ab^ve.  Fiber  orientation  in  the  caps  was 
primarily  unidirectional. 

All  panels  required  an  edge  bjildup  to  allow  for  mechanical  fastening  in 
the  test  fixture.  This  was  accomplished  as  part  of  the  layup  and  cure  of 
the  panel  skins.  Provisions  were  also  made  for  attachment  of  the  ends  of 
the  stiffener  webs  to  the  fixture. 

Fabrication  of  the  shaker  specimens  war.  basically  the  same  as  for  the 
sonic  fatigue  test  panels.  For  the  shaker  specimens,  however,  it  was  more 
efficient  to  fabricate  several  large  panels  and  subsequently  cut  them  into 
the  required  size  for  the  individual  specimens. 
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A  quantity  of  25  basic  tools  plus  2  rate  tools  were  designed  and 
fabricated.  Tooling  for  the  shaker  specimens  was  minimal,  with  only  two 
assembly  bond  jigs  and  one  2  section  layup  tool  required.  The  flat  sonic 
fatigue  panels  required  eight  tools.  The  remaining  16  tools  were  for  the 
curved  panels. 

Figures  11  through  17  show  photographs  of  the  tools  and  layup.  Figure  11 
shows  the  bond  jig  for  the  zee  layup.  The  zees  were  made  by  cutting  hat 
sections  in  half.  The  beam  in  the  center  of  the  tool  is  for  the  hat 
layup.  Around  the  edge  of  the  tool  is  a  rubber  tube  bonded  in  place  with 
RTV .  This  acted  as  a  vacuum  seal  for  the  rubber  bag  (Figure  13). 

Figure  12  shows  the  same  tool  with  the  graphite-epoxy  fabric  laid  down  f>r 
the  hat  section.  Figure  13  shows  the  tool  with  the  silicone  ruober  bag  in 
place.  This  was  a  reusable  bag,  which  effects  cost  savings  and  improves 
laminate  surface  condition,  as  compared  to  using  disposable  bags. 

Figure  14  shows  a  pair  of  cured  hat  sections.  Figure  15  is  an  assaubly 
tool,  which  was  used  to  locate  the  zees  on  a  skin  laminate.  The  two 
beams,  when  bolted  down,  locate  the  zee  section.  Figure  16  is  the 
-.  '.'responding  locating  tool  for  the  honeycomb  stiffeners.  Figure  17  shows 
a  layup  of  graphite/epoxy  prepreg  for  the  skin  elements. 

Coupon  tests  were  performed  on  each  layup  to  determine  resin  content, 
density,  fiber  content  and  void  content.  Two  1-inch  square  coupons  were 
used  for  each  layup.  Ultrasonic  C-scans  were  performed  on  each  specimen 
to  check  for  bondline  voids. 


During  specimen  fabrication,  problems  were  encountered  in  two  areas.  The 
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The  bond  tool  failed  to  generate  adequate  pressure  in  the  zee  radius  to  be 
adjacent  to  the  skins.  As  a  result,  this  radius  had  sporadic  areas  of 
surplus  resin  and  resin  starvation.  This  bond  pressure  deficiency 
resulted  in  inadequate  interlaminar  strength.  In  addition,  excess  resin 
areas  are  prone  to  surface  cracking,  which  in  turn  can  result  in  premature 
fatigue  crack  initiation  under  the  kind  of  severe  test  conditions  for 
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Figure  14.  Cured  Hat  Sections 


which  these  specimens  were  intended.  'lthough  there  were  no  reliable 
criteria  for  the  strength  requirements  in.  the  stiffener  radii,  the  fact 
that  zee  st.ffeners  are  not  symmetrical  sections  and  can  be  expected  to 
experience  some  rotation  under  random  acoustic  loading  conditions,  led  to 
them  being  rejected  as  unsuitable  for  sonic  fatigue  testing.  The  bond 
tool  was  then  modified  to  provide  more  effective  throw-in  blocks  in  the 
radii,  and  a  set  of  good  quality  stiffeners  was  then  fabricated.  Although 
this  problem  did  not  affect  the  sonic  fatigue  test  panels,  soma  of  the 
shaker  '■pecimens  had  already  been  completed  before  the  problem  was 
discovered.  The  effects  that  this  problen  had  on  some  of  the  shaker  test 
results  are  discussed  in  Section  III. 4. 

A  second,  and  more  serious,  problem  was  discovered  during  the  early  shaker 
tests.  The  first  specimens  tested  (see  Section  III. 5)  experienced 
premature  failure  in  the  adhesive  joint  between  the  stiffeners  and  the 
skins.  This  caused  considerable  concern,  since  the  fabrication  of  the 
shaker  specimens  had  been  completed,  and  fabrication  of  the  multi-bay 
panels  was  in  progress,  with  some  of  them  already  completed.  Visual 
examination  of  the  failed  adhesive  joints  revealed  excessive  porosity  In 
the  adhesive.  This  type  of  porosity,  consisting  of  a  large  number  of  very 
small  voids,  does  not  show  up  on  Otrasonic  C-scans.  In  addition,  there 
was  no  grau,  ite  fiber  pull-out  around  the  failed  joints.  Fabrication  of 
the  nulti-bay  panels  was  then  suspended,  and  a  thorough  investigation  of 
the  adhesive  bonding  problem  as  initiated.  The  investigation  centered  on 
an  ‘••■'fimination  of  the  bonding  process,  but  also  included  a  revaluation  of 
Af-tT7,  the  adhesive  selected  for  this  program.  AF 147  is  a  tough, 
elastomeric  adhesive  with  high  peel  strength.  It  is  used  extensively  on 
he  F16  airplane  and  has  good  strength  properties  over  the  temperature 
range  for  which  graphite-epoxy  structures  are  considered  suitable.  It  is, 
howc/er,  an  adhesive  that  is  particularly  susceptible  to  moisture  during 
fabrication 
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The  initial  investigation  of  the  bonding  process  and  the  condition  of  the 
adhesive  revealed  a  higher  moisture  content  than  was  considered 
acceptable.  An  additional  batch  of  shaker  test  specimens  was  fabricated 
after  additional  storing  of  the  adhesive  in  a  dessicator  in  order  to 
eliminate  any  moisture.  Some  of  these  specimens  were  statically  tested  by 
simply  pulling  the  skin  and  stringer  sections  apart.  A  10  percent 
increase  in  static  strength  was  obtained,  compared  to  the  original  shaker 
specimens.  Comparative  shaker  tests  were  then  carried  out.  Some  riveted 
specimens  were  added  to  this  comparative  test,  in  order  to  provide  a 
reference  to  which  the  bonded  specimens  could  be  compared.  The  results  of 
this  comparative  testing  were  disappointing.  Although  the  "dried" 
adhesive  produced  significant  fiber  pull-out  upon  failure,  the  fatigue 
life  did  not  significantly  increase.  Table  16  shows  the  results. 


As  a  result  of  this  adhesive  problem,  and  also  because  of  the  extreme 
importance  of  stiffener  and  skin  laminate  quality,,  a  thorough 
investigation  was  carried  out  to  determine  the  history  and  quality  of  all 
the  sonic  fatigue  pane,  components  and  assemblies,  including  aoditional 
assessments  of  adhesive  bond  quality.  The  following  tasks  were  performed: 


(1)  The  weight  percent  of  resin  and  void  percent  by  volume  of  all 
skins  were  determined  and  tabulated.  The  results  are  given 
in  Table  17.  Acceptance  criteria  for  this  program  required  a 
resin  content  of  27  percent.  Consequently,  skins  bl  and  cl 
were  rejected,  q  was  considered  marginal  at  26  percent. 
However,  the  flatwise  tension  tests  performed  on  the  rejected 
cl  skin  {see  Item  3  below)  produced  good  results,  indicating 


that  q  was  an  acceptable  laminate. 


(2)  Sections  of  2  stiffeners  were  cut,  mounted  and  photographed. 
The  .tees  made  on  the  original  tool  were  found  to  be  resin 
rich.  However,  the  zees  used  for  the  sonic  fatigue  test 
panels  were  found  to  be  representative  of  production  quality. 
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TABLE  16 

SHAKER  TEST  RESULTS  FOR  ADHESIVE  EVALUATION 
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(All  ,ype  2) 

Overall  RMS 
Stress  Level 
(lb/in2) 

Cycles  to 
Failure 

Riveted 

7,370 

234,000 

Bonded  -  dried 

a. 040 

3.000 

Bonded  -  dried 

6.432 

b.000 

Bonded  -  original 

Gauge  Lost 

5.400 

Bonded  -  original 

9.179 

5.400 

TABLE  17 

VOID  AND  RESIN  CONTENTS  OF  SKIN  LAMINATES 


Skin 

No.  Plies 

Average 
Resin  (ij 

Averaqe 

Void  U) 

Average 
Density  (1) 

a  1 

6 

2b 

.64 

1.61 

a-2 

b 

29 

.4/ 

1 .60 

b-1 

8 

24 

1  .20 

1.63 

b-2 

8 

.11 

-0- 

1.61 

c-1 

8 

74 

1.29 

1 .63 

« -2 

8 

29 

.1/5 

1 .60 

d 

77 

.31 

1 .63 

e-1 

8 

13 

.26 

1.58 

e-2 

H 

10 

.21 

1.60 

f-l 

8 

31 

.57 

1.59 

t  •  2 

8 

11 

.46 

!  .59 

9-1 

8 

79 

.5/5 

1  .61 

0-2 

K 

27 

.79 

1  .6'' 

h 

8 

2y 

.11 

1 .6- 

) 

a 

33 

.03 

1  .58 

j 

6 

27 

.155 

1  .63 

k-1 

8 

31 

-0- 

1 .60 

K  -l 

8 

31 

.765 

1  .58 

1 

1 2 

29 

.224 

1.61 

m 

8 

2H 

.16 

1  .62 

ri 

4 

28 

1.37 

1.60 

P 

4 

27 

.63 

1.63 

U 

4 

2b 

.58 

1.62 

r 

6 

29 

.36 

1.62 

s 

4 

27 

1 .595 

1.59 
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(3)  Flatwise  t  ision  tests  were  performed  on  the  rejected  cl  skin 
laminate  an  1  on  a  specially  fabricated  d  (3.2  ply)  laminate. 

The  results  are  given  in  Table  18.  They  show  the  c  laminate 
to  have  comparable  strength  to  the  d  laminate,  indicating  the 
percent  resin  content  criterion  to  be  conservative. 

(4)  A  range  of  shaker  specimens  were  Fabricated  using  different 

adhesives  and  different  processes,  for  comparison  vith  the 
original  specimens  The  adhesives  used  were  AF147,  the 
current  selection,  and  FM1000.  FK10G0  is  an  older  adhesive 
that  has  excellent  strength  properties  and  is  easy  to  use. 
However,  it  is  environmentally  susceptible  and  is  not  widely 
used  in  production.  It  is,  however,  an  excellent  reference 
adhesive.  Using  two  plies  of  AF147  was  also  evaluated. 

Table  19  shows  the  results  of  the  static  tests  and  Table  20 
shows  the  shaker  test  results.  The  AF147  was  found  to  have 
superior  static  strength,  but  FM1G00  did  better  in  fatigue. 

It  was  also  clear  that  a  second  ply  of  AF147  resulted  in  a 
significant  improvement  in  fatigue  life. 

Following  the  above  tests,  the  failed  static  and  shaker  test  specimens 
bonded  with  AF147  were  found  to  have  porosity  uniformly  dispersed  in  the 
weave  pattern  of  the  knitted  fabric  in  the  bond  line.  Additional  testing 
was  then  performed  in  order  to  determine  the  cause  of  this  porosity. 

These  tests  included  the  comparative  evaluation  of  (1)  solvent  wiping 
subsequent  to  grit  blasting,  (2)  no  solvent.,  just  dusting  with  a  clean  dry 
cloth,  (3)  an  evaluation  of  the  amount  of  vacuum  used  during  bagging  and 
curing,  (4)  oven  drying  of  composite  details  and  glass  cloth  (used  as  air 
bleeder)  and  (5)  evaluation  of  weight  loss  during  oven  drying.  All  of  the 
above  were  evaluated  through  lap  shear  testing  and  visual  examination  of 
failure  mode.  Volatile  contents  determinations  were  also  made.  The 
results  are  given  in  Table  21.  In  the  lap  shear  tests.  Process  2  gave  the 
highest  failing  stress,  but  more  importantly,  the  bond  line  porosity  was 
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TABLE  18 

FLATWISE  TENSILE  TEST  RESULTS  ON  SKIN  LAMINATES 


Specimens  from  Panel  c-1  (Flatwise  Tensile) 


Specimen  No. 


Failing  Stress  (psi ) 


3480 


Resin  content  23.6%  by  weight,  voids  '..29%  by  volume. 

Specimens  c-1  -1  and  c-1 -2  failed  between  surface  plies. 
Specimen  c-1-3  Failed  approxi  .lately  in  the  center  of  the 
laminate . 


Specimens  from  Panel  d  Noted  in  Item  3 
(Flatwise  Tensile) 


Specimen  No. 


Failing  Stress  (psi) 


Resin  content  32.8%  by  weight,  voids  0.21%  by  volume. 

All  failures  occurred  at  the  approximate  center  of  the 
laminate . 


COMPARATIVE  ADHESIVE  EVALUATION  -  STATIC  RESULTS 


filing 


Load  (Lbs.) 


Adhesive 


220 

205 
168 

206 
242 
186 
200 
202 
220 


Original  Lot  AF-14/,  dried  (43  hrs.) 
Original  Lot  AF -147,  as  received 
Second  Lot  AF- 147,  Dried  (120  hrs.) 


Second  Lot  AF-147, 
Second  Lot  AF-147,  2 
FM-ICOO 

Original  Lot  AF-147, 
Original  Lot  AF-147, 
Original  Lot  AF-147, 


ply,  dried 
plies,  dried 

cut  from  Panel  c-1 
cut  from  Panel  c-1 
cut  from  Panel  c-1 


TABLE  20 

COMPARATIVE  ADHESIVE  EVALUATION  -  SHAKER  TEST  RESULTS 


Time  To  Failure 

^Minutes) 

Adhesive 

12b 

FM-1000 

37 

AF-147,  1  plv,  dried 

b3 

AF-147,  2  plies,  dried 

29 

AF-147  -  original  group 

C9 


TABLE  21 

AF-147  ADHESIVE  PROCESSING  EVALUATION 


1.  One-half  inch  Overlap  Lap  Shear  Specimens  (Adherends  cut  from  Panel  c-1). 


Process  Ava.  Failing  Stress  (psi) 

1  3240 

2  3430 

3  3065 

Process 


1  Adherends  grit  blasted,  wiped  with  MEK  and  air  dried 

30  minutes;  25”  Mg  vacuum  used  during  bag  check  and  cure. 

2  Adherends  grit  blasted  and  wiped  with  clean  dry  cloth. 
Adherends  and  glass  breather  cloth  oven  baked  at  150°F 
for  45  minutes.  10"  vacuum  used  during  bag  check  and 
panel  vented  to  atmosphere  during  cure. 

3  Same  as  Process  2,  except  25"  Hg  vacuum  was  applied  to 
assembly  throughout  cure. 

Volatile  Content  Determination 

Four  adhesive  specimens  were  cut  from  the  roll,  placed  in  a  200nF  oven., 
withdrawn  at  the  noted  intervals  and  weighed. 


Spec.  No. 


Time  at  200° F  %  Weight  Change  (Decrease) 


1  15  mins.  0.53 

2  45  mins.  0.55 

3  90  mins.  0.49 

4  240  mins.  0.65 
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reduced  oy  approximately  75  percent  from  the  original  specimens. 

Processes  1  and  3  exhibited  excessive  porosity,  similar  to  that  of  the 
original  specimens.  Volatile  content  was  determined  to  be  within  normal 
limits  for  adhesive  films.  The  major  factor  in  the  porosity  problem  was 
concluded  to  be  the  amount  of  vacuum  used  durinq  bag  check  and  the  lack  of 
subsequent  venting  to  atmosphere  durinq  cure.  The  pulling  of  vacuum 
during  curing  combined  with  the  presence  of  slight  moisture  is  what  caused 
the  poor  bond  quality.  All  subsequent  assemblies  were  then  tabricateo 
with  the  AF 147  adhesive  system  (single  ply)  using  the  optimized  bonding 
process.  Subsequent  sonic  fatigue  tests  on  joints  utilizing  Process  l 
showed  order-of-magnitude  improvements  in  sonic  fatiuue  life  over  joints 
utilizing  Process  1.  It  is  interesting  to  note  that  a  major  improvement 
in  joint  quality  relative  to  porosity  and  random  fatigue  life  corresponded 
to  a  very  modest  improvement  in  static  strength. 
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This  adhesive  evaluation  underscores  the  crucial  importance  of  bond 
quality  in  a  program  of  this  type,  and  the  need  to  rectify  any  proDlems 
prior  to  fabricating  a  large  number  of  expensive  test  structures. 
Fortunately,  in  this  program,  the  problem  was  discovered  and  rectified 
prior  to  the  fabrication  of  most  of  the  multi-bay  test  panels.  Those  that 
had  already  been  fabricated  were  eventually  subjected  to  sonic  fatigue 
testing  (see  Section  III. 5),  and  tailed  prematurely  in  the  bonded  joints. 
They  were  subsequently  replaced  with  new  panels  and  successfully  tested. 
Unfortunately,  the  fabrication  of  the  shaker  specimens  had  already  been 
completed  prior  to  the  resolution  of  the  bonding  problem,  and  in  addition, 
many  specimens  were  used  in  the  process  of  achieving  a  solution.  This 
resulted  in  a  major  change  in  the  objective  of  the  shaker  tests.  It  was 
not  now  possible  to  use  the  existing  specimens  to  evaluate  the  skin- 
stiffener  joint  fatigue  proper. ies.  Instead,  they  were  used  to  evaluate 
the  skin  laminate  fatigue  properties,  which  were  unaffected  by  the  bonding 
problem. 
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NOTE:  It  was  later  learned  that  the  manufacturer  of  the  AF147  adhesive 
had  been  having  problems  with  air  porosity  in  this  adhesive,  and 
that  this  had  been  a  contributing  factor  in  the  bonding  problem 
experienced  by  Rohr.  The  manufacturer,  like  Rohr,  has  now  overcome 
the  problem. 

4.  SHAKER  TESTS 

The  shaker  test  program  was  originally  intended  to  provide  additional 
fatigue  life  and  mode  of  failure  data  on  the  skin-stringer  adhesively 
bonded  joints.  The  data  was  to  augment  the  sonic  fatigue  test  data  from 
the  multi-bay  panels.  However,  as  a  result  of  the  adhesive  bonding 
problems,  discussed  in  Section  III. 3,  tests  were  redefined  in  objective 
and  srope. 

Early  shaker  tests  revealed  poor  bond  quality  between  the  skin  and 
stringer  elements.  The  resulting  investigation  indicated  that  the 
ranaining  specimens  would  be  similarly  deficient.  Part  of  this  resulting 
investigation  consisted  of  performing  shaker  tests  on  some  of  the  original 
specimens  and  comparing  the  results  with  those  from  a. variety  of  new 
specimens  utilizing  different  adhesives  and  process  parameters.  In  this 
endeavor,  the  shaker  tests  proved  to  be  a  valuable  aid  in  both  discovering 
the  bond  problem  and  in  evaluating  solutions. 

One  of  the  major  justifications  for  shaker  testing  the  skin-stringer 
joints  was  the  belief  that  the  modes  of  failure  and  cycles  to  failure 
would  correlate  with  the  sonic  fatigue  test  results.  This  belief  turned 
out  to  be  fully  justified.  The  shaker  tests  that  revealed  the  poor  bond 
quality  were  characterized  by  rapid  failures,  with  stiffeners  completely 
delaminating  from  the  skins,  with  virtually  no  graphite  fibers  being 
pulled  from  the  skin  laminates  and  occurring  at  relatively  low  strain 
levels.  Some  early  sonic  fatigue  tests  on  panels  made  prior  to  the 
resolution  of  the  bonding  problem  displayed  the  same  failure 
characteristics.  Subsequent  testing  of  specimens  having  good  bond  quality 
resulted  in  considerable  skin  laminate  damage  prior  to  and  during  skin- 
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stringer  joint  failure  in  both  the  shaker  and  the  sonic  fatigue  tests. 
Additionally,  shaker  tests  performed  on  the  riveted  specimens  resulted  in 
failures  in  the  zee  radius  adjacent  to  the  skin.  This  same  mode  of 
failure  occurred  during  sonic  fatigue  testing  of  riveted  multi-bay  panels. 

The  specimens  were  ganged  together  in  groups  of  five,  mechanically 
fastened  through  the  stringer  webs  to  the  upstanding  leg  of  a  horizontal 
tee  bar,  and  subjected  to  1/3  octave  random  loading  centered  around  the 
specimen  response  frequency  and  tested  on  a  Ling  B290  shaker  having  a 
capacity  of  1,500  force  pounds.  This  method  was  specifically  intended  to 
primarily  load  the  skin-stringer  joint.  Following  the  discovery  and 
resolution  of  the  bonding  problem,  it  was  decided  to  use  the  remaining 
shaker  specimens  to  develop  random  fatigue  data  for  the  skin  laminates. 

In  order  to  accomplish  this,  the  test  fixture  was  modified  to  support  the 
skin  elements  as  cantilevers,  making  sure  that  the  skin-stringer  joint  was 
well  away  from  the  point  of  maximum  stress  on  the  skin.  The  skins  were 
mounted  in  tapered  blocks  in  order  to  avoid  abrupt  changes  in  stiffness, 
and  strain-gauged  at  the  point  of  expected  maximum  strain.  Testing  was 
then  carried  out  as  before,  with  1/3  octave  random  loading.  Figure  18 
shows  the  shaker  test  setup. 


Complications  arose  during  the  early  tests  due  to  the  high  strains 
required  to  cause  fatigue  failures  of  the  laminates.  Conventional  strain 
gauges  do  not  have  significant  fatigue  life  at  the  strains  required  to 
fail  the  graphite  laminates.  To  overcome  this  problem,  the  specimen 
holding  fixture  was  strain-gauged  and  tests  carried  out  to  establish  a 
relationship  between  the  fixture  gauge  and  the  specimen  gauges.  This  was 
done  using  static  loading  of  up  to  4,000  microinches/inch  on  the  specimen, 
and  noting  the  corresponding  fixture  gauge  readings.  Correspondence 
between  the  specimen  and  fixture  strain  gauges  was  also  determined  by 
apply  a  sinusoidal  load  up  to  a  specimen  strain  of  3,000  microinches/ inch. 

Finally,  similar  correspondence  was  also  established  using  low  level 
random  excitation.  The  fixture  gauge  was  found  to  read  approximately  1/20 
of  the  specimen  gauges.  Although  there  was  little  variation  from  specimen 
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Figure  18.  Shaker  Test  Set-Up 


to  specimen,  nevertheless  strain  conversion  factors  were  measured  for  each 
test. 

The  remaining  shaker  specimens  were  then  tested  to  failure.  The  results 
are  shown  in  Figures  19  and  20.  Figure  19  shows  actual  rms  strains  vs. 
cycles  to  failure.  The  curve  drawn  represents  minimum  values.  The 
numbers  in  parentheses  refer  to  the  specimen  types  given  in  Figure  B-3  in 
Appendix  B.  Since  the  elastic  modulus  varies  between  different  specimen 
types,  there  are  advantages  in  presenting  data  in  strain  form,  allowing 
users  to  apply  their  own  modulus  values.  Figure  20  shows  the  same  fatigue 
data  plotted  as  rms  stress  vs.  cycles  to  failure.  These  curves  are  used 
in  Section  IV  in  conjunction  with  the  progressive-wave  tube  test  results. 

5.  PROGRESSIVE-WAVE  TUBE  TESTS 

Sonic  fatigue  tests  were  performed  on  the  twenty-seven  panels  shown  in 
Table  1  and  existing  panels  1,  2,  4  and  5  shown  on  Figure  B-2  in 
Appendix  B.  The  tests  were  carried  out  in  a  progressive-wave  tube  (PWT) 
at  the  Acoustic  Test  Facility,  Rockwell  International  (Los  Angeles 
Aircraft  Division),  Los  Angeles,  California.  The  facility  is  powered  by 
four  Ling  EPT  200  transducers,  each  capable  of  generating  10,000  acoustic 
watts.  Sine  and  random  inputs  are  available  with  frequency  spectrum 
control  from  50  Hz  to  1,200  Hz.  Indefinite  endurance  tests  can  be  carried 
out  at  overall  sound  pressure  levels  of  167  to  168  dB.  The  main  test 
section  is  in  a  6-foot  by  1-foot  duct  cross-section,  capable  of  taking  two 
panels  simultaneously,  one  above  the  other.  The  PWT  has  an  acoustic  wedge 
termination  into  a  reverberation  room.  Rockwell  personnel  operated  the 
PWT.  All  instrumentation,  data  acquisition,-  signal  conditioning  and  data 
reduction  were  performed  by  Rohr  personnel  using  Rohr's  mobile  Vibro- 
Acoustic  Laboratory.  Figure  21  shows  the  Rockwell  facility  and  the  Rohr 
mobile  laboratory.  Figure  22  shows  the  test  section  with  panels 
installed.  The  main  purpose  of  the  tests  was  to  obtain  strain  and 
frequency  response  data  for  the  test  panels  under  random  acoustic 
excitation  at  grazing  incidence,  and  to  test  the  panels  to  failure,  using 
the  data  generated  to  develop  a  sonic  fatigue  design  method. 
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a.  Evaluation  of  Progressive-Wave  Tube  (PWT)  -  Prior  to  the 
actual  test  program,  a  series  of  measurements  were  made  in  the  Rockwell 
PWT  in  order  to  determine  the  maximum  overall  and  spectrum  acoustic  levels 
available,  without  noticeable  "clipping".  Amplitude  distribution  plots 
were  also  made  in  order  to  determine  if  the  acoustic  field  was  reasonably 
Gaussian.  Three  microphones  were  used,  in  a  vertical  spread  along  the 
center  line  of  the  test  panel  openings.  The  results  showed  acoustic 
spectrum  levels  of  around  140  to  145  dB/Hz  to  be  attainable  with  broad¬ 
band  loading.  Reducing  the  acoustic  loading  spectrum  to  1/3-octave  showed 
an  increase  in  maximum  acoustic  spectrum  levels  of  approximately  8  to 

10  dB.  Figure  23  shows  the  amplitude  distribution  function  for  the  center 
microphone  with  broad-band  input.  These  data  showed  the  Rockwell  facility 
to  be  suitable  for  this  program. 

b.  Instrumentation  -  All  the  test  panels  were  instrumented  with 
sufficient  strain  gauges  and  microphones  to  accurately  identify  dynamic 
strains,  mode  shapes  and  acoustic  loading.  The  center  bay  of  each  panel 
was  the  most  heavily  strain-gauged.  All  panels  hen  biaxial  gauges  at  the 
center  of  the  center  bay  and  adjacent  to  both  zees  on  the  longer  sides. 
Strain  gauge  and  microphone  locations  are  shown  on  Figure  B-l  in 
Appendix  B.  Sheet  1  shows  locations  for  panel  "b."  Referring  to  the 
numbering  system  for  panel  b,  the  other  test  panels  were  instrumented  as 
fol lows: 

Panels  a,  c,  d,  f,  g,  h,  n.  p  had  strain  gauges  at  positions  3,  4,  7,  8, 
10,  11,  18,  26  and  31.  In  addition,  panels  a,  b,  d,  t.  and  p  had  back-to- 
back  gauges  for  positions  3,  4  and  10. 

Panel  e  was  gauged  as  per  Figure  B-i,  sheet  2. 

Panels  i,  j,  q  and  r  had  gauges  in  positions  3  and  4  on  each  of  the  four 
center  bays,  plus  positions  7,  8,  10  and  11  on  one  center  bay,  plus 
position  31. 
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Panels  k,  1  and  s  had  gauges  at  positions  3  and  4  on  each  of  the  two 
center  bays,  plus  positions  7,  8,  10  and  11  on  one  center  day,  plus 
position  31. 

Panel  m  had  gauges  at  positions  1  through  13  on  one  of  the  center  bays, 
plus  a  gauge  corresponding  to  "4"  on  the  remaining  seven  bays. 

The  J  stiffened  panels  were  instrumented  as  panel  b. 

Back-to-back  gauges  were  used  on  selected  panels  in  order  to  separate  out 
membrane  and  flexural  strains.  Strain  gauges  were  also  installed  on  ttie 
test  fixture  to  check  for  unwanted  resonances. 


Small  (1/8-inch)  strain  gauges  were  used  to  provide  for  good  resolution. 
Larger  strain  gauges  result  in  excessive  strain  averaging,  particularly 
near  stiffeners  and  fixtures  where  there  are  high  strain  gradients.  Since 
these  locations  are  where  maximum  strains  and  fatigue  failures  occur,  good 
resolution  is  of  particular  importance.  Each  panel  had  two  flush-mounted 
microphones  installed.  "Kulite"  pressure  transducers  were  used.  Several 
panel/fixture  assemblies  had  extra  microphone  holes  provided  to  facilitate 
acoustic  measurements  on  all  four  panel  sides.  "Kul ite"  transducers  are  a 
strain  gauge  ype  microphone  and  therefore  used  compatible  signal 
conditioning  to  that  used  for  the  strain  gauges.  Their  high  natural 
frequency  (above  70  KHz)  and  low  mass  makes  them  especially  suitable  for 
mounting  on  vibrating  structures.  Three  B&K  condenser  microphones  were 
installed  and  monitored  inside  the  PWT  as  part  of  the  facility  operation. 


c.  Data  Acquisition  -  The  data  acquisition  consisted  of  twertty 
channels  of  strain  gauge  signal  conditioners,  coupled  through  a  patch 
panel  t.o  a  14-track  FM  tape  recorder.  Two  channels  were  set  up  for 
handling  microphone  (Kulite)  signals.  Prior  to  each  test  run,  insertion 
calibrations  for  all  data  channels  were  recorded  on  magnetic  tape.  These 
insertion  calibrations  consist  of  applying  a  calibration  resistor  in 
parallel  across  each  strain  gauge  to  simulate  a  known  compressive  strain. 
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Post  tost  calibrations  were  also  performed  as  a  check,  and  as  a  safeguard 
against  neglecting  gain  changes  made  during  test  runs. 

d.  Test  Procedure  -  Each  panel  was  installed  in  the  progressive- 
wave  tube  and  subjected  to  acoustic  loading  at  grazing  incidence.  The 
panel -fixture  assemblies  were  suspended  on  wires,  in  order  to  isolate  them 
from  PWT  vibrations  and  to  achieve  accurate  boundary  condition 
representation  for  comparison  with  the  analytical  results.  Figure  24 
shows  a  close-up  view  of  the  panel-fixture  installation  in  the  PWT.  Load 
cells  were  incorporated  into  the  wire  harness  supporting  the  panels.  This 
allowed  the  wire  tension  to  be  adjusted  identically  for  each  test  panel 
and  also  facilitated  dynamic  monitoring  to  ensure  that  there  were  no 
significant  resonances  in  the  panel  suspension  system.  The  turnbuck.le- 
pulley  arrangement,  seen  in  Fiqure  24,  automatically  centered  the  test 
panels  in  the  specimen  windows. 

The  test  procedure  for  each  panel  started  with  a  sine  sweep  from  bO  Hz  to 
1 ,200  Hz.  Thr  sine  sweep  was  used  to  identify  major  panel  resonances. 

This  was  followed  by  full  spectrum  (50  Hz  fo  1,200  Hz)  random  acoustic 
loading  from  140  dB  to  165  dB  iri  5  dB  steps  at  BO  second  intervals.  All 
strain  gauge  and  microphone  outputs  were  recorded  on  magnetic  tape 
throughout.  In  addition,  real-time  frequency  response  plots  were  rnaae  for 
one  key  strain  gauge  and  microphone.  Where  the  number  of  transducers  for 
a  given  panel  exceeded  the  14  channels  available  on  the  tape  recorder, 
these  runs  were  repeated  until  all  transducer  outputs,  including  the 
fixture  gauges  and  the  load  cell,  had  been  recorded  on  magnetic  tape.  The 
overall  strain  levels  from  all  the  strain  gauges  were  monitored 
throughout.  When  the.  strains  reached  levels  suitable  for  endurance 
testing,  the  random  response  check  would  not  proceed  to  the  next  acoustic 
load  level.  Careful  response  monitoring  is  of  particular  importance  in 
setting  the  test  levels  for  panel  endurance  runs.  The  strain  and  acoustic 
levels  measured  during  these  random  response  checks  form  the  data  base  for 
the  design  method  in  Section  IV.  Endurance  runs  were  then  made  at 
selected  sound  pressure  levels  until  panel  failure  occurred.  A  target  of 

82 


10  hours  exposure  time  was  set  as  a  maximum.  The  intention  in  the 
endurance  runs  was  to  avoid  rapid  or  protracted  failure  times  and  to 
obtain  a  good  spread  of  fatigue  life  data.  Panels  that  could  not  be 
failed  within  10  hours  at  165  dB  were  subsequently  subjected  to  reduced 
bandwidth  testing,  with  correspondingly  higher  acoustic  spectrum  levels. 
During  endurance  runs,  panel  gauges  were  continually  monitored  for  changes 
in  response  (frequency  or  strain  levels)  indicative  of  structural  failure. 
In  addition,  periodic  visual  inspections  of  the  panels  were  made.  First 
signs  of  visual  damage  were  noted. 

The  panels  often  had  slow  progressive  fiber  failures,  where  the  time  from 
first  visual  damage  to  major  damage  affecting  panel  response  was  several 
hours.  In  such  cases,  both  times  were  noted.  Major  damage  was  defined  as 
any  skin  damage  extending  through  the  laminate  thickness  or  fracture  or 
separation  of  stiffeners  from  the  skin  in  one  of  the  center  bays. 

In  addition  to  the  basic  tests  described  above,  the  following  additional 
tests  were  carried  out: 

(1)  Testing  identical  panels  in  steel  and  aluminum  fixture 
frames. 

(Z)  Testing  a  panel  without  a  fixture  fraim  and  also  rigidly 
bolted  to  the  PWT. 

(3)  Full  depth  vs.  panned  down  closures  on  the  honeycomb  beam 
stiffeners. 

(4)  Testing  with  and  without  the  stiffeners  clipped  to  the 
fixture  frame. 

(5)  Comparison  between  bonded  and  riveted  skin  to  stiffener 
joints. 
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(6)  Switching  panel  positions  in  the  two  test  windows. 

(7)  Measuring  response  on  one  panel,  with  other  test  window  open. 

(8)  Same  as  above  with  hard  wall  installed  in  other  test  window. 

e.  Data  Reduction  -  A  Spectral  Dynamics  Digital  Signal  Processor, 
Model  SD360  was  used  to  perform  all  spectral  analyses.  The  SD360  is  a 
self  contained  fast  fourier  transform  analyzer,  capable  of  displaying  and 
plotting,  in  real-time,  the  complex  relationship  of  two  signals,  both  in 
the  time  domain  and  the  frequency  domain.  An  analysis  range  of  1.2  KHz 
was  sel  cted,  corresponding  to  a  filter  bandwidth  of  2.16  Hz.  The  actual 
aliasing  filter  cutoff  was  960  Hz.  Overall  sound  pressure  levels  and  rins 
strain  levels  were  determined  by  converting  the  signal  to  a  d.c.  value 
proportional  to  its  instantaneous  rms  value,  integrating  over  a  20  second 
period,  and  reading  the  value  on  a  digital  voltmeter. 

Frequency  spectra  were  generated  for  all  microphones  and  strain  gauges  at 
the  endurance  test  sound  pressure  levels.  Spectra  were  also  generated  at 
each  sound  pressure  level  (140  dB  to  16b  dB)  for  selected  gauges.  The 
strain  gauges  selected  for  spectral  analysis  over  the  full  response  range 
were  the  center  biaxial  pair  (numbers  3  and  4)  and  the  gauge  near  the  zee 
radius  (number  10).  Overall  rms  levels  were  measured  for  all  transducers 
at  all  sound  pressure  levels.  Cross-spectral  density  measurenents  were 
made  between  corresponding  strain  gauges  on  adjacent  bays.  Integrated 
power  spectral  density  plots  were  made  for  some  panel  gauges  in  order  to 
determine  the  relative  contributions  of  individual  modes  to  the  overall 
rms  strain  value. 

f.  Progressive-Wave  Tube  Test  Results  -  The  first  panels  to  be 
tested  were  the  existing  panels,  shown  in  Figure  B-2.  Tables  22  and  23 
summarize  the  overall  rms  stress  levels.  Panels  1  and  5  have  the  same 
geometry,  and  offer  a  comparison  between  a  graphite  and  an  aluminum  panel. 
The  stresses  on  the  graphite  panel  (5)  ranged  from  50-75  percent  of  the 
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Stringer  Panel  5  -  Graphite  Skir.-Stringer  Panel  4  -  Single  Stiffener 


corresponding  aluminum  values  (1)  in  the  maximum  stress  direction  (short 
direction).  In  the  long  direction  (strain  gauge  15)  the  stresses  were 
comparable.  Panel  4,  which  has  four  times  the  bay  span  as  panel  5,  shows 
a  corresponding  stress  increase  of  from  2  to  4  times.  The  unstiffened 
panel  (3)  was  tested  both  as  a  flat  panel  and  with  a  curvature  of 
30-inches.  The  165  dB  data  points  may  not  be  valid  comparisons,  since  the 
panel  was  undergoing  extremely  large  deflections  at  this  high  load.  The 
165  dB  poi'its  may  also  represent  the  onset  of  failure.  The  center 
stresses  were  reduced  to  1/3  to  1/2  of  their  original  values  due  to 
cuivature,  whereas  the  edge  stresses  were  reduced  to  1/5  to  1/7  of  their 
original  values.  The  response  of  panels  1  and  5  were  plotted  against 
sound  pressure  levels  in  order  to  compare  the  degree  of  linear  response 
between  the  aluminum  and  graphite  panels.  Figures  25  and  26  show  the 
results.  The  dotted  lines  represent  linear  response.  These  graphs  show 
both  the  aluminum  and  the  graphite  panels  to  be  responding  in  a  linear 
fashion.  Figure  27  shows  the  graphite  skin-stringer  panel  following  sonic 
fatigue  failure.  Figure  28  shows  sections  of  honeycomb  stiffeners  with 
skin  laminate  fibers  still  attached.  The  time  to  failure  for  panel  5  was 
15  minutes  at  165  dB.  Panels  3  and  4  lasted  for  5  minutes. 

These  photographs  show  the  mode  of  failure  to  be  in  the  skin  laminate  at 
the  stiffener  locations.  This  shows  that  the  secondary  bond  between  the 
skin  and  stringers  is  superior  to  the  interlamina  bond  strength,  as  it 
should  be.  This  is  because  the  adhesive  strength  in  the  laminate  comes 
from  the  epoxy  matrix  material,  which  is  selected  for  criteria  other  than 
just  pure  strength.  Flow  characteristics,  for  example,  are  very  important 
when  laying  up  a  large  surface  area.  The  adhesives  used  to  bond  the  skin 
and  stringers  together  are  chosen  primarily  for  strength.  Consequently, 
extensive  fiber  pull-out  on  failure  is  indicative  of  good  bond  quality. 
This  mode  of  failure  also  indicates  that  the  flatwise  tension  strength  of 
the  laminate  may  be  a  critical  parameter  in  sonic  fatigue  resistance. 

This  is  a  property  that  is  not  commonly  measured  or  quoted  in  structural 
property  specifications  of  composite  laminates.  This  conclusion  has 
considerable  logical  appeal,  since  stress  concentrations  and  extra  inertia 


Figure  25.  Panel  1  -  Linear  Response  of  Aluminum  Skin-Stringer  Panel 


forces  are  experienced  by  the  surface  laminate  underneath  the  stiffeners 
during  dynamic  behavior.  Since  the  interlarnina  strength  will  be  the  same 
between  each  lamina,  failure  will  o^cur  at  the  first  bonded  interface. 

This  mode  of  failure  also  indicates  that  mechanically  fastening  the 
stiffeners  to  the  skins  may  result  in  longer  sonic  fatigue  life  than  using 
bonded  stiffeners.  This  is  because  fasteners  will  distribute  stresses 
across  the  whole  skin  laminate,  rather  than  just  into  the  surface 
laminate.  As  a  result  of  these  findings,  it  is  recommended  that  flatwise 
tension  tests  be  performed  on  composite  skin  laminates  in  the  future. 

Such  tests  may  provide  valuable  information  in  selecting  the  best  resin 
systems  and  adhesives  for  sonic  fatigue  critical  applications  of  advanced 
composites  structures. 

The  next,  and  most  important  phase  in  the  sonic  fatigue  test  program,  was 
to  test  the  multi-bay  panels  shown  in  Figure  B-l.  These  are  the  panel 
tests  upon  which  the  design  method  in  Section  IV  was  based.  The  first  of 
these  panels  were  fabricated  prior  to  the  resolution  of  the  adhesive 
bonding  problems,  discussed  in  Section  III. 3.  Some  of  these  panels  tailed 
prematurely,  with  the  stiffeners  delaminating  from  the  skins,  with  no 
fiber  pull-out  occurring.  The  response  data  for  these  panels  is 
unaffected  by  the  weak  bond,  but  the  times  to  failure  are  not 
representative  of  the  panel's  fatigue  lives.  The  panels  that  had  been 
fabricated  with  suspect  bonds  were:  al,  cl,  fl,  gl,  i,  kl,  n  and  q. 

Panel  al  gave  good  response  data  up  to  160  dB,  but  failed  prematurely  at 
166  dB,  with  very  slight  skin  damage.  Since  this  was  a  conf iguration  for 
which  there  was  a  planned  duplicate  panel  (a2)  yet  to  be  made,  there  was 
no  need  to  refabricate  al.  Panels  cl  and  gl  also  had  duplicate  panels 
scheduled,  and  it  was  decided  to  rivet  the  suspect  skin-stringer  joints  on 
panel  g  in  order  to  provide  a  comparison  between  bonded  and  riveted 
joints.  Panels  i  and  kl  failed  prematurely  with  no  skin  damage  occurring. 
These  panels  were  subsequently  riveted  back  together  and  retested  in  order 
to  provide  additional  response  comparisons  between  bonded  and  riveted 
specimens.  Panels  n  and  q  gave  good  response  data  and  failed  with 
significant  fibers  being  yulled  from  the  skin  laminate.  Panel  fl  gave 
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good  response  data  and  did  not  fail  after  9  hours  at  165  dB. 
Figure  29  shows  panel  al  following  sonic  fatigue  failure.  The 
skin-stringer  joint  areas  show  a  mixture  of  weak  bonding  (white 
areas)  with  no  attendant  fiber  pull-out  and  satisfactorily  bond¬ 
ed  areas  with  fiber  pull-out  occurring.  Figure  30  shows  panel  n 
following  sonic  fatigue  failure.  Here  the  suspect  bonding  proc¬ 
ess  does  not  seem  to  have  resulted  in  a  weak  joint,  and  the 
failure  shows  extensive  skin  laminat:  damage.  It  should  be 
pointed  out  that  skin  laminate  damage  is  a  desired  mode  of 
failure  and  represents  a  successful  test.  Figures  31  and  32 
provide  a  good  example  of  this  desired  mode  of  failure.  They 
show  the  front  and  back  faces  respectively  of  panel  p  follow¬ 
ing  sonic  fatigue  failure.  In  this  case  the  skin-stiffener 
bond  strength  and  the  skin  laminate  quality  are  well  demon¬ 
strated  by  the  even  distribution  of  the  failure  through  the 
entire  thickness  of  the  skin  laminate. 

Table  24  gives  the  overall  rms  strain  levels  for  those  panels 
whose  response  data  was  subsequently  used  in  the  development 
of  the  design  method.  Corresponding  response  spectra  are 
given  in  Appendix  C.  Omitted  from  this  table  are  those  panels 
whose  purpose  was  to  investigate  specific  effects,  outside  the 
main  design  method;  such  as  the  J  stiffened  panels,  the  honey¬ 
comb  stiffened  panels  (e),  panel  in  -  the  8x1  array  and  panel 
c  -  which  was  designed  to  investigate  the  effects  of  ply  ori¬ 
entation.  Strains  c re  given  in  microinches/inch  for  the 
following  strain  gauges:  3-  center  of  bay,  long  direction; 

4  -  center  of  bay,  short  direction;  and  10-edge,  short  direc¬ 
tion,  normal  to  longest  side.  Strain  gauge  10  gave  the  high¬ 
est  measured  strains  for  the  majority  of  the  panels,  and 
represents  the  location  of  maximum  interest  in  this  program. 

At  the  panel  centers,  where  the  response  strain  magnitudes  are 
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Figure  30.  Panel  n  after  Sonic  Fatigue  Testing 
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Figure  32.  Panel  p  after  Sonic  Fatigue  Testing  -  Front  Face 


comparable  in  both  directions,  it.  is  necessary  to  combine  their 
effects  using  the  relationship 
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and  v xy  and  VyX  are  from  Tables  2  through  6.  The  resulting  biaxial 
strains  are  also  given  in  Table  24. 


A  detailed  discussion  of  the  results  in  Table  24,  as  they  relate  to  the 

development  of  a  design  method,  is  given  in  Section  IV.  In  general,  the 

results  show  basic  logical  trends,  such  as  decreasing  strains  with 

increasing  skin  laminate  thickness;  and  increasing  strains  with  both 

increasing  stringer  spacing  and  increasing  radii  of  curvature.  There  are, 

however,  several  inconsistencies  in  the  data;  panel  a  is  identical  to 

panel  b  except  for  having  fewer  skin  plies,  yet  it  has  lower  edge  strains 

than  does  panel  b  at  the  higher  sound  pressure  levels.  Panel  n,  which  has 

even  fewer  skin  plies,  also  has  lower  strains  at  the  higher  sound  pressure 

levels.  However,  it  should  be  remembered  that  the  different  skin 

laminates  have  different  ply  orientations,  and  consequently,  different 

elastic  modulii.  Panels  a  and  d  have  an  elastic  modulus  of  7.5  x  10^ 

2  6  6 
lb/in  .  Panel  b  has  a  value  of  6.7  x  10  and  panel  n  has  9.4  x  10  . 

Although  this  program  compared  two  different  ply  orientations  for  the  same 

laminate  thickness  (panels  b  and  c),  this  variable  was  not  represented 

over  a  number  of  panels  sufficient  to  permit-  its  inclusion  as  a 

quantitative  variable  in  the  design  method.  Instead  typical  symmetric 

ply  orientations  were  chosen,  with  the  expectation  that  the  resultant 

design  methid  would  be  applicable  to  other  similar  laminates.  This 

limitation  should  be  remembered  if  radically  different  ply  orientations 

are  used  in  conjunction  with  the  results  of  this  program. 


Although  there  were  a  few  exceptions,  the  maximum  strain  response  on  the 
panels  occurred  at  t lie  center  of,  and  normal  to,  the  longest  bay  side, 
adjacent  to  the  radius  of  the  zee  (strain  gauge  10).  This  is  as  it  should 
be,  and  also  corresponds  to  the  finite-element;  static  analysis  results. 

The  measured  strains  did  not  linearly  increase  with  overall  sound  pressure 
level,  but  increased  at  a  lower  rate,  which  varied  from  panel  to  panel. 

It  is  not  entirely  clear  whether  or  not  this  is  indicative  of  nonlinear 
structural  response.  The  response  frequencies  of  the  in-ohase  stringer- 
bending  mode  for  the  panels  listed  in  Table  2d  were  given  in  Table  15, 
where  they  can  be  seen  to  correspond  quite  well  to  the  fully-fixed 
frequencies  calculated  from  Reference  5. 

fatigue  lives,  as  expected,  showed  consideraole  scatter.  Figure  33  shows 
the  fatigue  life  data  points  for  the  multi-bay  panels.  The  data  points 
are  shown  superimposed  on  the  shaker  test,  fatigue  data.  The  curve  shows 
the  strain  endurance  level  to  be  approximately  400  microiuches/ inch. 

Taking  a  conservative  line  through  the  data,  the  curve  for  the  ski ri - 
stiffener  joint  appears  to  be  approximately  42  percent  of  the  shaker  test 
curve  for  the  skin  laminate.  This  ratio  is  similar  to  that  for  riveted 
aluminum  skin-stiffener  panels.  The  lowest  strain  at  which  a  panel 
failure  occured  (excluding  the  defective  panels)  was  411  microinches/incn, 
and  that  appears  to  be  an  outlier  compared  to  the  other  data  points.  1  he 
next  lowest  failure  strains  were  444/446  microinches/inch,  occurring  at 
approximately  lQ^cycles.  Virtually  all  of  the  panels  (again  excluding  the 
defective  panels)  displayed  the  same  failure  mechanism.  The  first  signs 
of  failure  were  isolated  failed  skin  fibers  at  the  skin-stiffener  joints. 
The  number  of  failed  fibers  would  gradually  increase,  often  over  a  period 
of  several  hours,  without  having  any  effect  on  the  panel  response.  Only 
when  the  damaged  skin  fibers  had  extended  across  nearly  all  of  the  skin- 
stifferier  joints  was  a  change  is  response  detected.  This  would  usually  be 
closely  followed  by  a  major  failure  of  the  skin  laminate.  This  failure 
mechanism  presents  a  problem  in  defining  the  effective  fatigue  lives  of 
the  panels.  It  the  first  visible  sign  of  skin  damage  is  the  criterion 


101 


CYCLES  TO  FAILURE 


for  failure,  then  the  cluster  of  fatigue  data  points  shown  on  Fioure  33  in 
6  7  5  6 

the  10  -10  range  would  occur  in  the  5  x  10  to  10  range.  Also,  some 

panel  fatigue  data  points  shown  as  run-outs  had  very  slight  fiber  damage. 

These  effects  would  cause  the  assumed  endurance  level  to  drop  from  400  to 

300  microinches/inch.  From  a  structural  point  of  view,  failure  should  be 

defined  in  terms  of  significant  damage  or  a  reduction  in  load  carrying 

capability.  On  an  actual  aircraft,  however,  it  seems  likely  that  any 

structural  component  showing  visible  signs  of  damage  would  be  removed, 

even  if  the  damage  were  unlikely  to  propagate.  Several  sonic  fatigue  test 

panels  exhibited  a  small  number  of  fiber  failures  early  during  testing, 

but  did  not  experience  any  damage  propagation,  even  after  several  million 

more  cycles. 

Phase  and  cross-spectral  density  functions  were  generated  between 
corresponding  strain  gauges  in  adjacent  bays  in  order  to  identify  the 
stringer-bending,  in-phase  mode.  Figures  34  through  37  are  for  panel  r 
(6  x  3,  6  ply,  60-inch  radius).  Figure  34  shows  the  sine  sweep  at  the 
center  of  the  center  bay.  Figures  35  and  36  show  the  random  response 
spectra  for  the  centers  of  two  adjacent  center  bays.  From  these  spectra 
it  can  be  seen  that  the  major  response  modes  occur  at  350-430  Hz  and  at 
750-800  Hz.  Figure  37  shows  the  corresponding  phase  relationship  between 
the  adjacent  bays  (top  plot)  and  the  associated  cross-spectral  density 
function  (bottom  plot).  From  this  figure,  it  can  be  seen  that  the 
response  in  these  two  bays  is  coupled  at  the  major  response  peaks  (shown 
by  peaks  in  the  cross-spectral  density  function)  and  that  the  coupled 
response  peaks  at  360  Hz  and  400  Hz  are  18U  deg.  out  of  phase,  whereas  the 
response  at  760-780  Hz  is  the  in-phase  mode.  The  cross-spectral  density 
plots  also  assist  in  more  precisely  defining  the  coupled  mode  frequencies. 

Another  example  of  the  value  of  phase  and  cross-spectral  density  functions 
in  identifying  response  modes  is  shown  in  Figures  38,  39  and  40. 

Figures  38  and  39  show  response  spectra  for  adjacent  bays.  Figure  38 
shows  two  distinct  response  peaks,  at  240  Hz  and  400  Hz,  but  Figure  39 
only  shows  one  peak,  at  ?15  Hz.  From  these  two  plots  alone,  modal 
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Figure  39.  Random  Response  -  :Janel  s.  Gauge  4-2 
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identification  is  not  possible.  Figure  40  shows  the  corresponding  phase 
and  cross-spectral  density  plots.  The  phase  plot  shows  the  200  Hz  region 
to  be  out-of-phase.  But  since  the  response  spectra  (Figures  38  and  39)  do 
not  show  identical  peak  response  frequencies,  the  cross-spectral  density 
plot  is  required  in  ordi r  to  identify  the  frequencies  at  which  the 
response  is  coupled  between  the  adjacent  bays.  Then  it  can  be  clearly 
seen  that  the  out-of-phase  mode  occurs  at  200  Hz  and  the  in-phase  mode 
occurs  at  380  Hz. 

Table  25  compares  overall  rms  strain  levels  for  Z  and  J  stiffeners,  and 
also  between  a  quasi-isotropic  laminate  (b)  and  a  more  highly  oriented 
laminate  (c).  Strains  are  given  at  the  panel  centers  and  at  the  edges, 
over  a  range  of  overall  sound  pressure  levels.  Figures  41  through  48  give 
corresponding  strain  spectra  at  160  dB.  The  comparison  between  the  Z  and 
J  stiffeners  does  not  present  a  clear  picture.  At  the  panel  centers,  the 
strains  are  comparable  for  both  stiffener  types,  although  it  can  be  seen 
that  for  the  (b)  panel  at  160  dB  and  165  dB,  the  J.  stiffener  resulted  in 
higher  strains  than  did  the  Z.  Looking  at  the  spectra  on  Figures  41  and 
42  (note  scale  difference)  it  can  be  seen  that  the  higher  J  stiffener 
strains  are  due  primarily  to  the  response  peak  at  275  Hz.  The  fundamental 
170  Hz  peak  was  reduced  by  the  introduction  of  the  J  stiffener.  The  same 
comparison  for  panel  c  (Figures  43  and  44)  also  shows  the  J  stiffener 
effecting  a  reduction  in  response  at  the  fundamental  mode  frequency.  In 
this  case,  without  a  significant  increase  in  the  amplitudes  of  the  higher 
frequency  modes.  Corresponding  comparisons  of  the  edge  strains  show 
significant  response  reductions  on  panel  b,  but  not  on  panel  c.  The 
corresponding  spectra  for  panel  b  (Figures  45  and  46)  show  a  2:1  reduction 
in  the  peak  response  level  due  to  the  0  stiffeners.  Panel  c  also  shows  ?■ 
reduction  in  the  response  of  the  first  mode,  but  the  increased  response  in 
the  other  modes  results  in  an  increase  in  the  overall  strain  level.  In 
general,  it  is  clear  that,  compared  to  the  Z  stiffeners,  the  J  stiffeners 
resulted  in  a  lower  response  level  for  the  fundamental  stringer-bending 
mode,  but  had  a  tendency  to  stimulate  the  stringer-torsion  mode  at  280  Hz. 
Figure  49  shows  the  160-170  Hz  peak  to  be  in-phase  and  the  270-280  Hz  peak 
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Figure  47.  Edge  Strain,  Panel  c,  Z  Stiffener 
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to  be  out-of-phase.  The  J-stiffened  panels  did,  as  expected,  exhibit 
longer  fatigue  lives  than  dici  the  Z-stiffened  panels.  For  panel  b,  the 
J  stiffener  increased  the  fatigue  life  from  50  minutes  (for  the  Z  stiffener) 
to  7  hours.  For  panel  c,  the  fatigue  life  increased  from  1  hour  to  17 
hours.  This  large  increase  in  fatigue  life,  without  a  major  reduction  in 
response  was  due  to  the  increased  bonded  footprint  area  of  the  J  compared 
to  the  Z.  This  results  in  more  extensive  fiber  pull-out  from  the  skin 
laminate  upon  failure,  hence  there  is  slower  damage  propagation  and  longer 
fatigue  lives.  Based  on  these  results,  the  J  configuration  appears  to  be 
an  attractive  stiffener  concept. 

Table  25  and  the  corresponding  spectra  also  provide  for  a  comparison 
between  the  laminates  used  for  panels  b  and  c.  The  more  highly  oriented 
"c"  laminate  shows  a  significant  reduction  in  response  of  the  edge  strains 
(approximately  25  percent).  The  center  strains  are  similar  for  both 
laminates.  However,  the  reduction  in  edge  strain  does  not  occur  when 
using  the  J  stiffener.  The  corresponding  frequency  spectra  do  not  provide 
any  additional  information  on  the  response  difference  between  the  two 
laminates.  The  difference  between  the  two  laminates  involves  taking  t.wo 
of  the  eight  plies  running  in  the  1c.  i  bay  direction  (90  deg.,  in  the 
X  direction)  and  running  them  in  the  short  direction  (0  deg.,  Y  direct ioi- 
between  stringers).  Looking  at  the  overall  results,  it  was  concluded  that 
these  two  laminates  had  comparable  sonic  fatigue  resistance. 


Three  panel  configurations  (g,  i  and  k)  were  used  to  compare  the  response 
of  bonded  and  riveted  joints.  Table  26  summarizes  the  overall  rms  strain 
levels  at  the  center  and  edge  of  the  center  bay  of  each  panel.  The 
differences  between  strains  at  the  bay  center  do  not  appear  to  be 
significant.  There  is  a  tendency  for  the  edge  strains  to  be  a  little 
lower  on  the  riveted  panels  than  on  the  bonded  panels,  but  the  differences 
are  neither  large  nor  consistent.  Figure  50  shows  a  comparison  of 
response  spectra  for  panel  g  at  160  dB.  Although  the  overall  levels 
differ  (205  to  162  microinches/inch),  the  spectra  are  remarkably  similar. 
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Overall,  it  appears  that  the  riveted  joints  do  not  significantly  affect 
panel  response.  They  do,  however,  affect  the  mode  of  failure.  The 
riveted  panels  experienced  partial  failure  in  the  Z  radius  adjacent  to  the 
skin,  in  addition  to  skin  laminate  failures.  The  fatigue  lives  of  the 
riveted  panels  were  not  significantly  longer  than  for  the  bonded  panels 
with  good  quality  joints.  However,  observations  of  damage  propagation 
during  testing  indicated  that,  heavier  stiffeners  and  more  rigid  clipping 
at  the  stiffener  intersections  would  result  in  longer  fatigue  life  for  the 
riveted  panels.  This  is  because  the  riveted  panels  experienced 
significant  stiffener  and  clip  damage  prior  to  skin  failure.  This  was  not 
true  of  the  bonded  panels,  where  skin  damage  was  the  primary  mode  of 
failure.  Another  factor  to  bear  in  mind  is  that  a  slightly  substandard 
bond  is  difficult  to  detect  and  may  result  in  a  highly  premature  f'tigue 
failure.  A  slightly  substandard  rivet  joint  is  detectable  and  will  have  a 
less  severe  affect  on  fatigue  life. 

Figure  51  shows  the  strain  spectrum  corresponding  to  the  top  spectrum  on 
Figure  50,  without  the  stiffener  webs  clipped  to  the  fixture  frame. 
Although  the  overall  strain  level  increased  from  205  to  240 
microinches/inch,  due  to  removing  the  clips,  the  spectra  show  this 
increased  strain  to  be  predominantly  below  200  Hz.  This  indicates  that 
providing  proper  attachments  at  the  panel  boundaries  is  desirable  and 
reduces  the  low  frequency  overall  panel  motion. 

Early  testing  of  a  curved  panel  in  an  aluminum  fixture  resulted  in  a 
fixture  failure.  Steel  backup  fixtures  had  been  fabricated  for  such  an 
eventuality,  and  consequently,  it  became  necessary  to  determine  whether  or 
not  the  steel  versus  the  aluminum  fixtures  affected  panel  response. 

Figure  52  shows  the  results  of  this  comparison.  As  expected,  no 
significant  panel  response  effects  were  observed. 

Early  in  the  test  program,  several  side  experiments  were  performed  in 
order  to  ensure  that  testing  two  panels  simultaneously  would  not  produce 
any  unwanted  response  interrelationships.  One  of  the  panels  was 
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concurrently  tested  in  the  top  and  bottom  test  windows,  paired  with 
different  panels  and  also  paired  with  a  steel  plate  in  the  second  test 
window.  No  significant  effects  on  the  panel  responses  were  noted.  The 
second  test  window  was  left  open  at  one  point  and  the  only  noticeable 
effect  was  a  2-3  dB  drop  in  acoustic  levels  in  the  PWT.  When  the  driver 
outputs  were  increased  tc  bring  up  the  acoustic  level  to  that  previously 
used,  there  was  not  noticeable  change  in  the  response  characteristics  of 
the  test  panels  due  to  the  open  window. 

In  order  to  evaluate  some  boundary  condition  effects,  panel  d  was  tested 
for  response  under  three  different  edge  conditions.  First,  it  was 
suspended  in  the  test  window  without  a  fixture  frame,  i.e.„  with  free 
edges.  The  fixture  support  wires  supported  the  panel  at  the  four  corners 
only.  Then  the  panel  was  bolted  into  the  wall  of  the  PWT,  simulating 
fixed  edge  conditions.  Finally,  the  panel  was  supported  in  the  regular 
fixture  frame  and  suspended  on  the  fixture  support  wire,  as  used  cri  the 
remaining  panels.  The  effect  of  these  different  boundary  conditions  can 
be  seen  in  Table  27.  The  overall  rms  strains  show  that  the  free  edges 
result  in  lower  response  levels  than  when  using  the  test  fixture  or  fixed 
edges.  When  it  was  decided  to  perform  the  sonic  fatigue  tests  using  a 
picture-frame  fixture,  supported  on  wires,  it  was  hoped  that  the  response 
in  the  center  bay  would  be  the  same  as  for  fixed  edges.  The  overall 
strains  appear  to  confirm  this.  Figures  53,  54  and  55  show  corresponding 
edge  strain  spectra  at  165  dB  for  the  three  edge  conditions.  The 
stringer  bending  mode  occurs  at  340  Hz.  This  basic  mode  occurs  with  all 
three  boundary  conditions.  The  major  response  effects  of  the  boundary 
conditions  occur  at  160-180  Hz  and  at  250  Hz.  The  modes  at  these 
frequencies  do  not  occur  when  the  panel  is  freely  supported.  If  these 
extra  modes  were  due  to  vibrations  of  the  PWT  wall,  they  would  not  occur 
when  the  panel  is  in  the  fixture  frame,  supported  on  wires  unless  the  PWT 
v  ^rations  were  transmitted  through  the  wire  harness.  If  this  were  the 
case,  the  load-cell  in  the  wire  harness  would  detect  than.  Figure  56 
shows  the  corresponding  frequency  spectrum  for  the  load-cell.  As  can  be 
clearly  seen,  there  was  no  significant  dynamic  response  occurring  in  the 
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TABLE  77 

OVERALL  RMS  STRAIN  LEVELS  (MICROINCHES/INCH)  FOR 
PANEL  "d"  WITH  DIFFERENT  EDGE  CONDITIONS 


Sound 

Pressure 

Level 

(dB) 

Center  Strain  (Gauge  4) 

Edge  Strain  (Gauge  10)  j 

1 

Fixed 

Edges 

Edges 
Supported 
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Free 

Edges 

Fixed 
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23 

34 
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34 

145 

31 
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261 
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190 

296 

274 

165 

226 

372 

416 

281 

422 

425 

Figure  56.  Load  Cell  Spectrum 
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wire  harness  supporting  the  panel  fixture  assembly.  If  the  modes  at  the 
160-1SJ  Hz  and  250  Hz  were  associated  with  the  regular  fixture  frame,  they 
would  not  have  occurred  when  the  panel  was  bolted  to  the  PWT  wall  (for 
this  test  the  panel  was  first  removed  f  s  fixture  frame  and  then 
bolted  directly  to  the  PWT  wall).  The  rr<>.  -esponse  difference  between 
the  fixed  edge  conditions  and  the  regular  test  fixture  occurs  at  100  HZ 
and  below.  The  response  peaks  in  this  frequency  region  appear  to  be  due 
to  the  effects  of  the  PWT  wall.  Thus,  the  panel  support  system  chosen  for 
this  test  program  appears  to  be  structurally  similar  to  the  more  typical 
method  of  bolting  the  panel  to  the  PWT  wall,  without  picking  up  PWT 
dynamic  effects. 

Two  honeycomb  beam  stiffened  panels,  shown  in  Figure  B-l,  were  tested  and 
the  results  compared  to  an  equivalent  panel  with  Z  stiffeners  -  panel  b. 
The  results  are  shown  in  Table  28.  The  first  panel  (el)  had  panned  down 
closures  for  the  honeycomb  stiffeners,  as  shown  in  Figure  B-l.  This  panel 
failed  in  these  closures  after  1-1/2  minutes  at  165  dB.  The  second  panel 
(e2)  was  subsequently  modified,  replacing  the  existing  panned  down 
closures  with  full  depth  closures,  clipped  to  the  fixture  in  a  similar 
fashion  as  the  Z  stiffener  attachments.  This  panel  failed  after  5-10 
minutes  at  165  dB. 

Although  these  times  to  failure  seem  short,  t.he  maximum  strains  on  both 
panels  were  high,  such  that  if  superimposed  on  the  fatigue  curve  on 
Figure  33,  both  panels  appear  to  be  on,  or  slightly  above,  the  curve  drawn 
for  the  Z-stiffened  panels.  The  maximum  strains  on  both  el  and  e2  were 
higher  than  those  shown  in  Table  28  (which  are  presented  for  comparison 
purposes),  and  occurred  at  the  center  of  the  shorter  sides.  It  is  not 
known  why  the  honeycomb  stiffened  panels  had  different  maximum  reponse 
locations  compared  to  the  Z-stifferied  panels.  In  any  event,  the 
comparison  between  the  honeycomb  and  zee  stiffeners  was  inconclusive.  It 
does  seem  probable  that  the  honeycomb  stiffeners  would  have  been  more 
effective  if  the  core  material  had  had  greater  shear  strength. 
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Ideally,  a  sonic  fatigue  design  method  should  be  based  on  using  acoustic 
spectrum  levels  to  predict  strain  spectrum  levels  at  corresponding 
frequencies.  However,  Individual  response  spectrum  levels  usually  vary 
inconsistently,  compared  to  overall  response  levels,  making  them 
unsuitable  for  analysis  purposes.  An  example  of  this  is  shown  in 
Figure  57.  The  spectra  shown  are  for  the  same  strain  gauge  at  160  dB 
(top)  and  165  dB  (bottom).  The  overall  rms  strain  level  increased  from 
310  to  435  microinches/inch  due  to  the  5  dB  increase  in  acoustic  load. 
However,  the  strain  spectrum  level  at  the  major  response  mode  (170  Hz) 
actually  decreased  slightly  (from  82  to  68  microinches/inch) .  The 
increased  overall  strain  level  was  due  to  increases  in  the  strain  response 
at  other  frequencies.  Inconsistencies  of  this  type  make  it  impractical  to 
use  the  strain  spectrum  levels  in  the  development  of  the  design  method. 
This  leads  to  a  dilemma  in  the  treatment  of  random  strain  data.  Spectrum 
levels,  with  their  narrow  bandwidth,  often  vary  unpredictably;  whereas 
overall  levels,  with  their  wide  bandwidth,  include  response  that  does  not 
contribute  to  fatigue.  It  has  been  suggested  that  1/3-oct.ave  or  1-octave 
bandwidth  measurements  may  provide  the  necessary  stability  without  being 
influenced  by  superfluous  data.  Again  referring  to  Figure  57,  it  can  be 
seen  that  some  of  the  increase  in  overall  rms  strain  was  due  to  the  low 
amplitude  strains  in  the  250-1,000  Hz.  region.  This  is  quantitatively 
demonstrated  in  Figure  58.  Superimposed  on  the  power  spectral  density 
function,  is  the  integrated  power  spectral  density.  From  this  plot  it  can 
be  seen  that  the  major  response  peak  at  170  Hz  contributes  only  about 
20  percent  of  total  power  spectral  density  (mean  square).  Whereas  the  low 
amplitude  response  above  250  Hz  accounts  for  40  percent  of  the  total 
spectral  density.  Based  on  this  data,  it  might  be  thought  that  increasing 
the  sound  pressure  level  from  160  to  165  dB  would  not  necessarily  bring 
about  a  more  rapid  sonic  fatigue  failure.  However,  data  taken  during  the 
test  program  clearly  pointed  to  a  definite  relationship  between  the 
overall  rms  strains  and  fatigue  failures,  leading  to  the  conclusion  that 
some  low  amplitude  strains  that  might  be  thought  of  as  not  contributing  to 
fatigue  failure  do,  in  fact,  make  a  significant  contribution.  As  a 
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Figure  58.  Integrated  Power  Spectral  Density  for  Strain  Response 
for  Panel  b,  Gauge  4 


result,  overall  rms  strains  generated  in  the  PWT  test  program  were  used  in 
the  development  of  the  design  method,  described  in  Section  IV. 

Integrated  power  spectral  density  functions  were  generated  for  another 
purpose.  The  PWT  had  a  tendency  to  generate  an  acoustic  peak  in  the 
150-180  Hz  region  (see  Figures  C-l  through  C-6).  It  was  found  that  this 
peak  could  be  satisfactorily  controlled  by  adjusting  the  bias  voltage  in 
the  EPT  200  acoustic  drivers.  In  most  cases  the  peak  did  not 
significantly  affect  panel  response;  where  it  did,  integrated  PSD 
functions  were  used  to  quantify  the  effects.  Figure  59  shows  a  situation 
where  the  response  spectrum  makes  it  appear  as  if  the  acoustic  peak  has 
produced  a  major  response  peak  at  180  Hz.  However,  the  integrated  PSD 
curve  shows  that  only  14  percent  cf  the  energy  is  contained  within  that 
peak.  When  converted  to  rms  levels,  this  unwanted  peak  accounted  for  less 
than  10  percent  of  the  overall  rms  strain  level.  In  a  situation  where 
spurious  peaks  have  a  major  effect  on  response  levels,  the  integrated  PSD 
function  provides  a  quantitative  tool  for  subtracting  the  effects  out. 

The  strain  spectra  generated  during  the  progressive-wave  tube  tests 
(contained  in  Appendix  C)  were  used  to  determine  damping  ratios  for 
various  panels,  using  the  "half  power-point  bandwidth"  method.  The 
damping  ratio  is  found  from  the  following  relationship: 

c  -  & 

2f 

where  r.  -  Damping  ratio 

Af  =  Frequency  bandwidth  (Hz)  at  the  half-ampl itude 
point  of  major  strain  reponse  peak 
f  =  Center  frequency  of  major  strain  response  peak 

The  values  obtained  showed  considerable  scatter  and  did  not  show  any 
significant  correlation  with  strain  response  levels.  Consequently, 
damping  was  not  included  as  a  design  parameter  in  Section  IV.  Instead, 
the  use  of  a  typical  value  was  recommended. 
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Figure  59.  Integrated  Power  Spectral  Density  for  Panel  k.  Gauge  4-1 


Damping  ratio  values  obtained  ranged  from  0.02  to  0.08.  However,  the 
higher  values  occurred  when  more  than  one  response  mode  appeared  to  be 
contained  within  the  major  response  peak.  In  those  cases  where  the 
response  peak  was  clearly  a  single  mode  only,  damping  ratio  values  were 
typically  0.02  to  0.03.  These  damping  values  were  compared  to  values 
obtained  in  Reference  2  for  aluminur  skin-stringer  structures. 

Reference  2  quotes  values  of  0.010  to  0.018  using  the  "logarithmic 
decrement"  method.  However,  when  the  "half  power-point"  method  was  used 
on  response  spectra  in  Reference  2,  values  in  the  region  of  0.05  were 
obtained.  Based  on  these  observations,  it  was  concluded  that  the  damping 
characteristics  of  graphite  and  aluminum  panels  do  not  significantly 
d i f f er . 

Some  further  discussion  of  the  progressive-wave  tube  test  results,  as  they 
relate  to  the  development  of  the  design  method,  is  contained  in 
Section  IV. 
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SECTION  IV 

DEVELOPMENT  OF  DESIGN  METHOD 
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J. 
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1.  INTRODUCTION 

The  primary  objective  of  this  phase  of  the  program  was  to  utilize  the 

analytical  and  experimental  results  to  develop  a  practical  semi -empirical 

sonic  fatigue  design  method  for  graphite-epoxy  skin-stringer  panels. 

Measured  random  strains  from  the  sonic  fatigue  tests  were  compared  to 

those  calculated  from  Miles'  equation^,  using  as  inputs  the  static 

strains  calculated  from  the  finite-element  analyses.  The  test  results 

(5) 

were  also  compared  to  values  determined  from  the  AGARD  nomographs  for 
fully-fixed  edge  conditions.  Finally,  multiple  stepwise  regression 
analyses  were  performed  to  develop  empirical  relationships  between  the 
measured  strains  and  frequencies  and  various  combinations  of  panel 
configuration  parameter  and  finite-element  analysis  results.  From  these 
regression  analyses,  design  equations  were  developed  and  a  design 
nomograph  constructed.  A  worked  example  is  also  presented.  Section  IV. 5 
presents  the  design  method  and  nomographs  as  a  self-contained  unit, 
capable  of  being  used  independently  of  the  remainder  of  this  report. 
Appendix  C  contains  the  test  data  (overall  acoustic  and  strain  levels  and 
spectra)  used  in  the  development  of  the  design  method. 

An  early  problem  encountered  in  the  development  ot'  the  design  method 
involved  the  use  of  acoustic  spectrum  levels  as  loads.  In 
Paragraph  III .5 .f ,  reasons  for  not  using  strain  spectrum  levels  were 
given.  The  main  reason  is  that  although  the  response  spectrum  levels  show 
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logical  overall  trends,  the  individual  variations  from  one  data  point  to 
another  are  too  large  and  unpredictable  for  use  in  developing  a  design 
method.  Acoustic  spectrum  levels,  however,  while  exhibiting  a  certain 
degree  of  unpredictable  variation,  are  sometimes  consistent  enough  to 
facilitate  their  use  as  the  load  function.  However,  in  this  program  the 
acoustic  spectrum  levels  varied  in  such  a  way  as  to  invalidate  their  use 
as  a  regression  variable.  Figures  C-2  through  C-7  show  the  test  acoustic 
spectra.  These  spectra  correspond  to  flat,  1/3-octave  spectra. 
Consequently,  the  spectrum  levels  (1  Hz  bandwidths)  decrease  with 
increasing  frequency.  Since  panel  frequencies  increase  with  panel 
stiffness,  the  stiffer  panels  were  effectively  tested  at  lower  acoust  c 
spectrum  levels.  This  resulted  in  a  high  degree  of  interdependence 
between  the  acoustic  spectrum  levels  and  the  panel  configuration 
parameters,  thereby  violating  the  necessary  assumption  of  independent 
variables.  A  problem  similar  to  this  was  encountered  in  Reference  3.  In 
that  case  the  problem  was  overcome  by  dividing  the  load  into  the  dependent, 
variable  (measured  rms  strain).  When  that  was  attempted  in  this  program, 
the  resulting  regression  equations  showed  good  accuracy  and  satisfied  all 
the  usual  statistical  requirements  (F-values,  t-values,  Uurbin-Watson 
statistics,  etc.).  However,  when  these  equations  were  used  on 
combinations  of  panel  configurations  other  than  those  used  in  the  test 
program,  it  was  found  that  the  equations  were  numerically  dominated  by 
changes  in  the  acoustic  spectrum  levels,  and  were  not  sufficiently 
responsive  to  changes  in  panel  dimensions.  When  comparing  responses 
between  two  very  different  panels,  the  stiffer  panel  had  a  much  lower 
response  and  a  lower  acoustic  spectrum  'level  at  the  major  response 
frequency  than  did  the  less  stiff  panel.  The  regression  analysis  largely 
attributed  the  lower  response  to  the  reduced  acoustic  spectrum  level, 
rather  than  to  the  increased  panel  stiffness.  However,  a  review  of  the 
integrated  power  spectral  density  plots  showed  that  the  major  response 
peaks,  associated  with  the  pertinent  acoustic  spectrum  levels,  usually 
accounted  for  less  than  25  percent  of  the  overall  strain  response.  Thus, 
since  the  overall  sound  pressure  levels  did  not  vary  from  panel  to  panel, 
as  did  the  acoustic  spectrum  levels,  the  reduced  response  of  the  stiffer 
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panels  was,  in  fact,  due  largely  to  the  changes  in  panel  configuration 
paramenters.  Consequently,  regression  analyses  using  the  overall  sound 
pressure  levels  as  the  load  function,  resulted  in  acceptable  design 
equations . 

2.  SUMMARY  OF  ANALYTICAL  AND  EXPERIMENTAL  RESULTS 

(5) 

RMS  stresses  and  frequencies  calculated  using  the  AGARD  nomographs 
were  given  in  Table  9.  The  static  strains  and  stringer-bending  mode 
frequencies,  analytically  determined  from  the  finite-element  models  and 
the  NASTRAN  computer  program,  were  given  in  Tables  11  and  15, 
respectively.  Table  15  al|^listed  frequencies  calculated  from  the  AGARD 
nomograph,  Lin's  equations  and  those  measured  during  the  progressive- 
wave  tube  tests.  The  mis  strains  measured  during  the  progressive-wave 
tube  tests  that  are  pertinent  to  the  design  method  were  given  in  Table  24. 
Table  29  contains  calculated  and  measured  frequencies,  static  strains,  rms 
strains  calculated  using  Miles'  equation  and  measured  rms  strains.  These 
are  the  data  subsequently  used  in  the  regression  analyses.  For  the 
reasons  given  in  Paragraph  II. 5. d,  the  frequencies  computed  from  the 
finite-element  models  were  not  used  for  the  design  method.  The  strains 
calculated  using  Miles'  equation  have  very  high  values,  compared  to  the 
test  strains.  This  was  expected  and  is  due  to  using  the  overall  sound 
pressure  level  rather  than  the  spectrum  level  as  the  loal,  for  the  reasons 
given  in  Section  IV. 1.  The  only  purpose  of  generating  these  calculated 
strains  was  to  determine  if  there  was  a  consistent  relationship  between 
them  and  the  measured  strains. 

3.  REGRESSION  ANALYSIS 

Regression  analysis  is  a  statistical  method  for  investigating 
functional  relationships  between  variables,  based  on  sample  data.  It  is 
particularly  suitable  when  the  data  are  imprecise  and  there  is  a  need  to 
determine  optimum  relationships.  The  basic  approach  is  to  use  samples  of 
data  to  calculate  an  estimate  of  a  proposed  relationship  and  then  to 
evaluate  the  fit  using  statistics  such  as  "F"  and  "T." 


TABLE  29 

SUMMARY  OF  ANALYTICAL  AND  EXPERIMENTAL  RESULTS 
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5,195  56  47 

9,192  107  84 

16,287  172  133 

28,976  201  168 

51,559  293  258 


1,494  107  149 
2,647  184  262 
4,709  302  444 


The  regression  analyses  performed  in  this  section  utilized  a  modified 
multiple-stepwise  regression  computer  program.  Stepwise  regression 
involves  a  forward  selection  procedure  for  the  independent  variables,  with 
the  provision  for  eliminating  variables,  as  in  backward  elimination 
procedures.  The  program  analyzes  the  relationship  between  a  dependent 
variable  (measured  overall  rrns  strains)  and  a  set  of  independent,  variables 
(panel  configuration  parameters).  The  independent  variables  are  selected 
in  order  of  importance  for  entering  into  the  regression,  based  on  the 
reduction  of  sums  of  squares.  The  user  can  override  this  feature  and 
enter  the  independent,  variables  in  any  chosen  sequence.  The  program  has 
six  algebraic  transformations  available,  as  follows: 


Linear  y  =  a  +  b  *  +bx  +  .  bx 

112  2  n  n 


Log  y  =  a  b  x 
n  n 

Log  y  =  a  +  b  Log  x 
n  n 

2 

y  =  a  +  b  x  +  c  x 
n  n  n  n 

1/y  =  a  +  b  x 
n  n 

y  =  a  +  b  Log  x 
n  n 

y  =  a  +  b  /x 
n  n 

where  y  is  the  dependent  variable  and  x^,  and  ...  x  are  the 
independent  variables. 

For  each  variable  entered  the  program  computes  the  mean  values,  standard 
deviations  and  cross-correlation  coefficients.  For  each  step  in  the 
regression  analysis,  the  program  computes  and  lists  the  following: 


Sum  of  Squares  Reduced:  This  is  an  indication  of  the  amount 

of  points  summed  around  a  mean  value 
line  for  a  certain  step. 

Proportion  Reduced:  Indication  of  how  well  a  variable 

explains  the  regression  at  a  certain 
step. 

Cumulative  Sum  of  Indication  of  how  much  the  dependent 

Squares  Reduced:  variable  correlates  with  the  indepen¬ 

dent  variables  entered  at  that  point.. 

Cumulative  Proportion  Indication  of  how  well  the  indepen- 

Reduced:  dent  variables  explain  the  regression 

at  that  step. 

Multiple  Correlation  Indication  of  how  much  the  dependent 

Coefficient:  variable  correlates  with  the  indepen¬ 

dent  variables  entered  at  that  point. 

F-Value:  A  measure  of  the  scattering  of  values 

about  the  mean  accounted  for  by  the 
regression.  It  is  used  in  conjunc¬ 
tion  with  F-tables  to  determine  the 
degree  of  fit  of  the  regression 
equation . 

A  measure  of  the  dispersion  of  the 
observed  points  about  the  regression 
equation.  It  is  in  fact  the 
"standard  deviation  of  the 


Standard  Error  of 
Estimate: 


Regression 

Coefficients: 

Standard  Error  of 

Regression 

Coefficients: 

T-Values: 


Coefficients  of  the  regression 
equation. 

Indication  of  the  confidence  level 
for  the  regression  coefficients. 


Ratio  of  intercept  to  standard  error 
of  regression  coefficient.  It  is 
used  in  conjunction  with  T-tables  to 
determine  the  accuracy  of  the  corre¬ 
sponding  regression  coefficient. 


Table  of  Residuals: 


Difference  between  actual  and 
estimated  values  for  the  dependent 
variable. 


Durbin-Watson 

Statistic: 


Von-Neumann's  Ratio: 


Test  for  lack  of  autocorrelation 
between  error  terms.  A  bad  statistic 
is  indicative  of  an  independent 
variable  being  omitted. 

The  ratio  of  the  mean-square 
successive  difference  to  the  variance. 


The  following  is  a  description  of  the  sequence  of  regression  operations 
with  examples  of  compute,  program  outputs  used  to  develop  the  rms  strain 
nomograph  in  Section  IV. 5: 


The  input  data  were  of  the  form 


"rms  \ 
~SPL  ) 


=  F  (b,  t,  R) 


where 


Erms  is  the  overall  rms  strain  from  strain  gauge  10,  i.e., 
the  maximum  edge  strain. 

SPL  is  the  test  overall  sound  pressure  level  in  Ib/in 
corresponding  to  the  strain  value, 
b  is  the  stringer  spacing  in  inches, 
t  is  the  skin  laminate  thickness  in  thousandths  of  an 
inch,  and 

R  is  the  radius  of  curvature  in  inches. 

Table  30  lists  the  input  data  in  the  order  (left  to  right)  b,  t,  R  and 

(^SPr)  '  ^  bbe  b°ttom  of  the  table  are  the  means  and  standard 

deviations  of  the  input  data.  An  explanation  of  the  R  =  150  value  is 
given  later  in  this  section. 

The  program  then  computes  the  cross-correlation  coefficients  between  the 
variables,  giving  the  output  in  the  form  of  a  correlation  matrix.  This  is 
shown  at  the  top  of  Table  31.  Regression  then  proceeds  with  the  linear 
form  of  the  regression  equation:  y  =  a  +  t>nxn.  Each  variable  is  entered 
in  turn  and  a  corresponding  regression  coefficient  determined.  The 
program  also  computes  at  each  stage  the  parameters  shown  on  the  remainder 
of  Table  31.  Three  sets  of  statistics  are  shown,  one  for  each  of  the 
independent  variables.  The  bottom  set  represents  the  final  linear 
equation,  which  may  be  written: 

(~§?r)=  224,0  +  358,2  b  -51*65  t  +  11- 22 R  (19) 

The  program  then  used  this  equation  to  calculate  a  set  of  estimated  values 
for  the  dependent  variable  (G  rms  )  .  The  results  are  shown  in  Table  32. 

'“SPL  1 

As  can  be  clearly  seen  from  the  percent  deviation  column,  the  linear 
equation  does  not  have  acceptable  accuracy.  The  average  deviation  is 
given  as  40  percent. 


TABLE  30 


REGRESSION  INPUT  DATA 


TRANSFORMATION  CODE  =  L 

LINEAR  TYPE  CURVE.  OlNERAL  tOOAliCN  FURH  iS  Y  =  A  ■* 


INPUT 

CARL) 

ft 

I 

8  .  COO 

33 . 000 

130.000 

3506.800 

INPUT 

Card 

4 

2 

8. LOO 

44.0UL) 

4  30  »  OOG 

2665.500 

INPUT 

CARD 

ft 

3 

o  .  ODD 

tt. 030 

a  5  D  .  C  0  l 

2o28.800 

INPUT 

CARD 

n 

t 

o.OGG 

tt .  OOu 

1  u  Q  .  0  u  L‘ 

a 900. OOO 

INPUT 

LARD 

A 

3 

&.uOG 

tt. oco 

L3U.GQO 

2  19t.50Q 

INPUT 

CARD 

4 

t 

6. COO 

bb. 300 

150. DuO 

1155.200 

INPUT 

CARD 

t 

7 

D.  GOO 

66.000 

130. OOG 

1094.200 

INI  UT 

CARD 

V 

a 

U.uuu 

bb.OCO 

150. DOC 

1097.800 

INPUT 

CARD 

» 

9 

8  •  000 

6b. 000 

136.000 

1068.100 

INPUT 

CARD 

# 

10 

a. OCO 

bo.OliO 

1 30. ODC 

945.500 

1NPU  1 

CARD 

4 

1 1 

b.OOO 
t .  uO  0 

66 .000 

13U.  OC  u 

823.800 

INPUT 

Card 

ft 

12 

tt .ODD 

13O.GOO 

o86. 200 

INPUT 

ARD 

4 

13 

4.000 

tt. ODD 

l  ju . uoe 

VuV.000 

INPUT 

l.ARD 

ft 

It 

t.OOU 

tt.OOu 

130. uCD 

916.300 

INPUT 

CARD 

* 

15 

t.OOu 

tt. 000 

1 3u. OOu 

bio. 500 

INPUT 

CARD 

» 

1C 

t.OOC 

tt.OOG 

laO.OuO 

560.000 

INPUT 

CARD 

# 

17 

b.  DO  D 

tt. 0  DC 

l30.Duu 

i731.000 

INPUT 

CARD 

ft 

ID 

b  .  u  0  u 

tt.ODO 

15u.DOC' 

1605.700 

INPUT 

CARD 

4 

It 

u.LOu 

•♦4  .  u  0  0 

1 30. OoO 

1547.800 

INPUT 

CmRO 

4 

20 

8.000 

tt  •  0  uU 

1>0.00  j 

1292.600 

INPUT 

CARD 

ft 

2  1 

u.  GOO 

*.t .  LwO 

150. Ouu 

lul 7.200 

INPUT 

CARD 

• 

22 

b.OOu 

22.000 

1 30. 000 

ottt.ttOu 

INPUT 

Card 

ft 

tl 

8. COO 

22.OD0 

130.000 

5111.500 

INPUT 

CAKo 

* 

2.4 

a.ooo 

22.000 

15U.0U0 

3336.900 

INPUT 

CARD 

ft 

i5 

t.uOC 

22.000 

150.000 

1982.700 

INPUT 

CARD 

ft 

2o 

t.GOO 

22 . OOu 

130. OOl 

1782.70 0 

INPUT 

Cart 

ft 

2  7 

t.000 

22. OOu 

1 3L . OC  D 

1534.800 

INPUT 

card 

ft 

23 

t.COO 

22. OCO 

130.000 

l23b.2uu 

INPUT 

CARU 

ft 

2  V 

6 .  ODD 

tt . OGo 

3D .OOl 

7oa.9G0 

INPUT 

Card 

ft 

jU 

o.DOD 

tt.OOD 

.-<0.000 

576. VOD 

INPUT 

CARD 

ft 

31 

b.OOO 

44.000 

30.000 

527. 2D0 

INPUT 

CARD 

ft 

32 

8.000 

44.000 

30.u00 

476.100 

INPUT 

CARD 

ft 

33 

8.000 

44. Odd 

3u.0DC 

529.600 

INPUT 

CARU 

A 

34 

8.000 

44.000 

30. Ol  0 

479.400 

INPUT 

CARD 

ft 

35 

8.  DOC 

44.  L'OO 

bG.LOu 

1738.600 

INPUT 

LARD 

ft 

3b 

o  .000 

tt. 000 

bu.DOO 

1761*. 500 

INPUT 

CARD 

ft 

->7 

b.OOO 

4t, OOl 

6D.0CC- 

Ictb.bOO 

INPUT 

CARD 

ft 

38 

8.00  c 

tt. OOu 

bO.OCo 

1533.700 

INPUT 

CARD 

ft 

3V 

b.OOO 

tt. 000 

60. 00 v 

1379.600 

INPUT 

CARD 

ft 

to 

8.000 

tt.OOO 

oO.CO. 

1133.900 

INPUT- 

card 

CARD 

ft 

tl 

8.000 

tt.ODO 

VC. OOu 

18b5.500 

INPUT 

ft 

t*. 

8.000 

tt.OCO 

90. OOl 

1752 .700 

INPUT 

CARD 

ft 

t.T 

b.OOD 

tt. ulL 

tO.OCC 

1605.200 

INPUT 

CARD 

ft 

tt 

8.  OuG 

tt .ODD 

VO. OOu 

1 >96.500 

INPUT 

CARD 

ft 

t5 

0.  000 

bb. OoD 

Vu.fiUL' 

963.500 

INPUT 

CARD 

ft 

tb 

b.OOO 

oo .000 

90.000 

628.800 

INPUT 

CARD 

ft 

t7 

b.  CCD 

bb.OCO 

VO. 000 

742. 4C0 

INPUT 

CARD 

ft 

ta 

b.OOL 

bb.OOQ 

VD.uOO 

667.500 

INPUT 

CAfcD 

ft 

t9 

0.  000 

66.000 

90.000 

618.300 

INPUT 

Caro 

ft 

50 

6.000 

66. 0  CO 

90. OOU 

349.000 

INPUT 

CARL. 

ft 

51 

8.000 

22.000 

90.000 

3317. 200 

INPUT 

CARD 

ft 

j2 

o.0')0 

22 .Old 

9D.CU0 

2226.900 

INPUT 

CARD 

ft 

33 

8  .  uC  u 

22. Odd 

VO. OOO 

2746. 700 

INPUT 

CARO 

ft 

54 

4.  UOL1 

33.000 

60.000 

1186.300 

INPUT 

CARD 

ft 

33 

t.ODO 

33.000 

60.000 

1172.800 

INPUT 

CARD 

ft 

5b 

4.000 

33.  OOO 

60.000 

9B0.300 

INPUT 

Carl1 

tt 

37 

4.000 

33.000 

&O.OQO 

915.200 

INPUT 

CARO 

ft 

58 

tloOO 

33. OCu 

foO.GuO 

BbSIiGO 

INPUT 

CARD 

ft 

59 

b.OOO 

2. . OCu 

3U.000 

1744.200 

INPUT 

CARD 

ft 

fcC 

6.  OOO 

2 2.000 

jU.OOu 

1619.500 

INPUT 

CARD 

ft 

ol 

6.  ODD 

22 .OUU 

30.000 

I0O6.IOO 

INPUT 

CARD 

ft 

62 

b.OOO 

22.000 

30. OOC 

1532.400 

VARIABLE  MEAN  STANDARD 
NO„  DEVIATION 

1  6.61290  1.64322 

Z  42.22581  it.Dt.27V 

3  102 .09077  t7 .  4993t 

A  155b.958Lb  lu79 .0981 1 


TABLE  31 


REGRESSION  RESULTS  FOR  LINEAR  EQUATION 


CORRELATION 

ROM  1  b 
l.COOOO 

MATRIX 

t 

C. 22549 

R  . 

-0.1 5435 

/€rms\ 

rSHV 

0.32089 

ROM 

0.2*549 

1  .ODCOO 

0.14783 

-0 .50589 

ROM 

3 

-0.13435 

0.14783 

l.COUCO 

0.3166b 

ROM 

4 

0.32089 

-0.5C3&9 

D. 3 io6d 

1.00000 

VAR  1  AbLL  ENTERED . 2 


SUM  up  SCUAKc. S  KfcuULfcU  IN  IhlS  SUP....  l«j.i?Bf3a.i4S 
PROPORTION  RtUUt.hU  IN  This  STfcP .  C.*56 


CUMULATIVE  SUM  bh  SUUAKE3  RcbUCtU .  IB  1787 j6. 145 

CUMULATIV t  PROPORTION  REUUCtO .  0.25b  Of-  710316*0.831 


FOR  1  VARIABLES  ENTERED 

MULTIPLE  LOKRElATION  COEFFICIENT...  0.306 
(ADJUSTED  PGR  U.F.i. ..........  U.506 

F- VALUE  FOR  ANALYSIS  UE  VaRIaNCl...  Eu  .037 

STANCARL'  ERROR  uF  cSIIHAIL.. . VSo.SSj. 

(ADJUST  Lb  EUR  O.P.I... .  938.55 3 


VAKlAbLt  REORlSSION 

NUMBER  tuEFf ICILNT 

*  -37 .*3063041 

INTERCEPT  31*9. 03145833 


STD.  ERROR  bp 
Rill.  LOEFF. 
b. 19553 


COMPUTED 

T-VALUt 

-4.543 


VARIABLE  ENTEREO .  1 

SUM  OF  SOUAHES  REDUCE  0  IN  IhIs  STEP....  141369l9.o0l 
PRO  PUR  T  1UN  Kf:  UCED  IN  THIS  STcP .  C.19V 

CUMULATIVE  iUM  UP  SCUAREi  KEDUCtD. . ... .  31337557.746 

CUMULATIVE  PROPORTION  REDUCED .  0.453  OF  710316*0.831 

FUR  2  VARIABLES  ENTLREO 

MULTIPLE  CORRELATION  COEFFICIENT...  0 .675 

(ADJUSTED  FOR  D.F.I .  0.66C 

F-VALCt  FUR  ANALYSIS  OF  VARIANCE...  *4.654 
STANDARD  ERROR  up  ESTIMATE .  605.033 

(adjusted  Fur  Ci.f.)....  .......  biu.ds*. 

VARIABLE  REGRESSION  STU.  ERROR  OF  COMPUTED 

NUMBER  COEFFICIENT  RED .  COEFF.  T-VALUE 

2  -44.B35&BBJC  7.*5u49  -6.177 

I  300.9453 1125  04.76514  4.646 

INTERCEPT  1 4  60 .0 Soul.' 6 57 


VARIABLE  ENTERED . 3 


SyM,u^§t>lJ*RES  KEDULEO  IN  THIS  SltP.. 
PR  JPuRTlQN  REDUCE D  IN  THIS  STEP....*; 

CUMULATIVE  SUM  CIF  SwUARES  REDUCED.... 
CUMULATIVE  PROPORTION  REDUCED. . . 

FOR  3  VARIABLES  ENTERED 

MULTIPLE  CORRELATION  COEFFICIENT.. . 

(ADJUSTED  FOR  U.F.)..„, . 

^yftkUlrF?R  ANALYSIS  UP  VARIANCE... 

STANDARD  ERROR  uF  EST1MAT t . . 

(ADJUSltD  FOR  D.F. >..... . 


11*41  .’<886.655 
0.231 

43  76l’»44.<»0i 

0.686 


0.829 
0  .622 
4c. 331 
619 .691 
630 .066 


OF  710316*0. Q3I 


VARIABLE 

NUMBER 

2 

l 

3 

INTERCEPT 


REGRESSION 
COEFFICIENT 
-  51 .6548564C 
35b, 2*011376 
11,2 16629b 6 
224.01654015 


S  ft) .  ERROR  OF 
REl».  CUfcFF. 
3.6:>CV* 
40. 32664 
1.71503 


COMPUTED 
T-VALUt 
-9.141 
7.11B 
6.  ->40 


TABLE  32 


TABLE  OF  RESIDUALS  FOR  LINEAR  REGRESSION  EQUATION 


TA&Lt  ut-  RESIDUALS 


L  NO. 

V  VALUE 

V  ESTlhATL 

KfcOIUUAL 

1  OtViAl 

1 

3306.80000 

30  67 . 7  6  >3  3 

*39. 07443 

12.52006 

i 

26o5 .50000 

2499. 5t2l 3 

163.97787 

6 . cc6b9 

3 

282b  .80000 

2499.5c2.i3 

329.277a7 

11.64020 

4 

2900.0J000 

2499. 322.13 

400.47787 

13.8093b 

5 

2194.50000 

2499.52213 

—30  5.02213 

—13.89939 

6 

1155.20000 

1365.11529 

-207.91329 

-17.99621 

7 

1094.20000 

1503.11529 

— 2o«/»91  529 

-24.57643 

8 

1097.00000 

1  363  .  i  1  3c* 

-cos. 31529 

—24 . Ib79l 

9 

1  Cm-  h  .  1COOO 

1  505.  1  1 529 

-c9 5.01 329 

-27.O/.03? 

10 

945.50C 

1363.11529 

-4l 7.61329 

-44.10872 

11 

023  .eOOOO 

1363.115^9 

-539.31329 

— 65.46677 

1^ 

386 .26000 

IQbo • b09b  7 

-180.40967 

— 20.35 7o7 

13 

909.6Gu0o 

1366.60967 

-157.00967 

-17.26140 

14 

916.30000 

1 066.6096/ 

-130.30967 

-16.40398 

15 

816.50000 

1066.60967 

-250.10967 

-30.03193 

16 

680 .00000 

1066 • 609of 

-4B0.6G9o7 

-83.89822 

17 

1731.00000 

1763.06590 

-32.06590 

-3.00765 

18 

1605 . 70006 

1  7  63. 0659k/ 

-177.30390 

-11.04302 

19 

15*»7 .30000 

1  ?83.0t*59u 

—238.26590 

-13.20002 

20 

1/192  .60000 

1763.06590 

-490.46390 

-37.94414 

21 

1017.20000 

1763.C6590 

-7o5.fao590 

-75.2*157 

22 

6444 .60000 

3635. 92b97 

2bC 3.87103 

43.58353 

23 

511 1 .30000 

3635.92697 

147  3.57103 

28.86767 

24 

3836  .90060 

3635.9^o9? 

— _99. Oc89  7 

-8.96128 

25 

1902. 70000 

2c33.0i.j51 

-ccO. 31651 

-11.x 2194 

c6 

17o2.7(jOOO 

CC  03.0165  1 

— 4cb.31b3i 

-j-3.5  77  32 

27 

1534.60000 

c.203.0  lu51 

-oob. 21651 

-43.53709 

28 

1236.20000 

2  203. 01651 

—966.01651 

-18.20875 

29 

76b .90000 

1 153. 52665 

-584. 62664 

-50.02297 

3C 

376.90000 

1 1 53.  5/_uo4 

-  57 6.62664 

-99.95262 

31 

527 .2000 j 

1 153.52664 

—oco. 32664 

-lib. 80247 

3*. 

476.10000 

1 1  ->3.  5ct/64 

—67  7.4  1 064 

-142.20663 

33 

529.oC06J 

1 133. 3  coo4 

-oc  2.92bo4 

-1*.  7 .bli.92 

34 

479.4C0C-O 

1 1 35. 5/664 

—6 / 4. 12664 

-140.61882 

35 

175b .60000 

l99C.0c.53X 

260.37449 

15.27206 

36 

1786.50000 

1h90. 02531 

c96. 47449 

16.39527 

37 

1645 .60000 

1490  .Oc.931 

155.57449 

9.45397 

j/U 

1533.70006 

1*» 90. 04.551 

43.67449 

2.94766 

39 

l379.60GuJ 

1440 . 0  c  331 

-110.42551 

-8.00417 

46 

H33.VC.GuG 

1490.02351 

-336.12551 

-29.12952 

41 

1665.50000 

l6cv>.  5c-i  JO 

36.97562 

2 .06329 

4c 

1752.7OCu0 

lb2o. 52433 

-7 3.82436 

-4.21204 

43 

1605.2000j 

Ibcu. 5c438 

-22 1.32436 

-13.78796 

44 

1396.30000 

1 bc6 . 5 c43ti 

-430.02436 

-30.79301 

45 

965.50000 

—20. 33 too 

991.6386b 

102. 7279b 

46 

H28.B00CC 

—CO. 3566b 

855„13ftb& 

103.17793 

47 

742 .46000 

-2c. .  33b6d 

7t  b. 7366b 

103.54776 

48 

6<>7 .30066 

-4.6.  3  Jhbb 

b9  3 . b3o6B 

IX  J.9*3»7 

49 

6 1  o  . 3 0 0  66 

-2b.33oob 

o*4»638o8 

104.25986 

50 

549.00000 

-4.6.  33666 

37  3.33bo6 

104.7*737 

51 

3317 .20600 

2  9czc.  9  31c. /. 

3>4.  C.tib76 

xu. o7975 

52 

2226.9C&6J 

2962. 93122 

-736.03122 

-33.05183 

53 

2746.70000 

c962. 93122 

-216. ?3122 

-7.67240 

54 

1186.50000 

625. 3164b 

3b 1 .10352 

47.29739 

55 

1172 .oOOOO 

625.31646 

547.46352 

40.68175 

56 

980.50000 

OC3. 3 1146 

334.96352 

36.21 X 1 2 

57 

915.60003 

6c5. 31u4o 

■Lb  9. 8b 3 32 

3 1 .67434 

58 

b63.500wj 

623.31048 

lib. 16352 

27.58350 

59 

1744.20000 

1  373. 4 7  i  j 

170.72275 

9.7o603 

60 

1619.50000 

1573.47723 

46.0227a 

2.64179 

bl 

1606.1C000 

1 373.47725 

3c.bc 273 

2.03118 

02 

1532.40000 

1  *.73.4772.3 

-41.07725 

-2 . 6  CiL  58 

AVERAOL  4 

DEVIATION 

IS 

40.2  loC2 

DUKB IN— M A  1 3 UN 

STATliTU. 

lj 

l.l 1VB3 

VON  NtUMANN  *5  RATIO 

IS 

1. 13<j19 

lr»3 


The  program  repeats  this  analysis  procedure  for  each  of  the  algebraic 
transformations  listed  above.  They  can  then  be  compared  with  each  other 
and  the  best  one  selected.  In  this  case,  the  best  transformation  was  the 
Log-Log  form: 

Log  y  =  a  +  b  Log  x 
n  n 

Table  33  gives  the  mean  and  standard  deviations  of  the  transformed 
variables  and  also  the  cross-correlation  matrix. 

Table  34  gives  the  statistical  parameters  calculated  as  each  variable  is 
entered.  At  the  bottom  of  the  table  are  the  coefficients  for  the  final 
regression  equation: 


Log  =  3.0806  +  1.1045  Log  (b)-1.2069  Log  (t)  +  0.5519  Log  (R)  (20) 

Where  is  in  microinches/inch ,  b  is  in  inches,  t  is  in  thousandths 
of  inches,  R  is  in  inches  and  SPL  is  in  lb/in^. 


This  is  the  equation  used  to  generate  the  rms  strain  nomograph  in 
Section  IV .5. 

Table  35  gives  the  corresponding  residuals,  showing  an  average  deviation 
of  23  percent. 

There  are  many  aspects  of  the  regression  analysis  that  can  be  observed 
from  the  data  contained  in  Tables  30  through  35.  The  emphasis  in  this 
program  is  not  on  the  formal  statistical  tests,  but  on  relating  the 
regression  results  to  what  is  known  about  the  data. 

Referring  to  Table  34,  it  is  seen  that  as  each  successive  variable  is 
entered,  the  multiple  correlation  coefficient  and  F-value  increased,  while 
the  standard  error  dec> eased.  The  T-values  also  increased.  If  by  adding 
oris  of  these  variables  (b,  t  or  R),  or  an  extra  variable,  the  statistics 
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TABLE  33 


CROSS-CORRELATION  MATRIX  FOR  VARIABLES  USED  IN  FINAL 
REGRESSION  EQUATION 


TRANSFORMATION  CODE  =  2 

LOGARITHMIC  TYPE  CURVE.  GENERAL  EUUAT jluN  FORM  IS 


VARIABLE 

MEAN 

STANDARD 

NO. 

DEVIATION 

1 

0.80409 

0.12160 

2 

1.5979! 

0.16212 

3 

1.94624 

0.265bl 

9 

3.U>W 

0.26083 

LOG  Y  ~  A  *■  ft  (COG  Xt. 


CORRELATION 

b 

ROW  1 

1.00000 

ROW  c 

0.23993 

ROW  3 

-0.19990 

ROW  9 

0.26480 


MATRIX 

-t 

C .23993 
l.OUOOO 
0.15b46 


R 

— 0. 19990 
0.13C9O 
1.00000 
0.3668b 


("Sr) 

0.26980 

— O . 36220 


0 . 36886 

1 .oooco 
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TABLE  34 


REGRESSION  ANALYSIS  RESULTS  FOR  FINAL  REGRESSION  EQUATION 


VARIABLE  ENTERED . 2. 


SUM  OF  SL DARES  ReuUlEU  IN  THIS  SUP....  la.213 

PROPORTION  RE  DOLE  D  IN  THIS  SUP .  C.Slo 

CUMULATIVE  SUM  OP  SQUARES  RELULEO .  1.213 

CUMULATIVE  PROPORTION  REDULEt .  0..31G  OP 

PGR  1  VARIABLES  ENTERED 

MULTIPLE  CURRELATIJN  COEFFICIENT...  t .562 

(ADJUST  tO  PUR  La  F* ) • «>• 0 . 5  6  c 
P-VALUE  FOR  ANALYSIS  ClF  VARIANCE... 

STANDARD  ERROR  uF  fc ST la'tAT L . -.  u.209 

(ADJUSTED  FDR  U.P.) .  u.«:0'* 


VARIABLE 

NUMBER 

2 

INTERCEPT 


RlGRE  SSlUN 
OUtFF ICIENT 
-O.BtVtlSiA 
<»»5G:>32B12 


STD.  ERROR  jF 
REG-  LLEPP. 
O.lBbAB 


COMPUTED 
T -VALUE 
*S»t66 


3.333 


VARIABLE  EN1ERED . 3 


SUM  OF  SQUARES  Kc.JULlO  IN  THIS  STEP....  D.790 

PROPORTION  RLDUCLO  IN  THIS  ST  t  P .  O.rOG 

CUMULATIVE  SUM  vjF  SQUARES  REUUCtL. .  . .  .  .  2.003 

CUMULATIVE  PRuPGRlION  REOuCED .  0.>22  UP 

FOR  i  VARIABLES  cNTtRLD 

MULllPLt  CORRELATION  OOF  F  F1C  I  ENT  .  .  .  0.722 

(ADJUSTED  POK  O.F.) _ .......  0.717 

F-VALUF.  FOR  ANALYSIS  OF  VARIANLE...  32.202 

STANDARD  ERRuK  uF  estimate. ........  U.lVo 

(  AO  JUST  EO  FOR  D.P.) .  0.170 


VARIABLE 

NUMBER 

2 

3 

INTtRLEPI 


REGRESSION 
COEFFICIENT 
-0  .9«2‘‘.t»!’'U 
D  -tuD  :>i7tT 
i  «UCosGV<  ti 


STD.  ERROR  OF 
KEG .  Cut  F  r „ 
0.  lnlOfc 
O.O&VhO 


computed 

T-VALUE 
-O. V6S 
•j.DAO 


3.B3R 


VARIABLE  ENTERED.....  1 


SOM  UF  SQUARES  REDUCED  IN  THU  ilEP....  3  .OCC 

PROPORTION  REDUCED  IN  THIS  ST  c  P .  0.261 


CUMULATIVE  SUM  OF  SQUARES  REUuCEO. . . . . .  3.003 

CuMUEATl V  t  PROPORTION  REDUCED.  . .  0.7B3  OP 


PUR  3  VARIABLES  ENTERED 

MULTIPLE  lUKRELATIUN  COEFFICIENT...  O.USb 

(  AD  JUS  1  ED  FOR  D.P.) . .  U.oBG 

F-VALUt  FOR  ANALYSIS  OF  VaKIANlc...  b9.3bO 
SIANUAKU  ERROR  ur  L S i THAT t ........ .  0.120 

(ADJUSTED  FUK  U.F.) . .  0.12*. 


3.B3B 


VARIABLE 

NUMBER 

2 

3 

I 

INtEKCEPT 


RLGRtSSION 

luEFFICIENT 

-I .20t.502Sd 
Q  .6  ?1  (*.333 
1  .  S33C  u 
j .ObGbl 3 v 1 


STD.  ERROR  OP 
REG.  LutEF. 

0.0W6U 
0  aUtiiy  I 

6  .  l3*.Ab 


COMPUTED 
T-VALUt 
-12. U2 
6.901 
E..339 
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TABLE  35 


TABLE  OF  RESIDUALS  FOR  FINAL  REGRESSION  EQUATION 


SELL-CTION. 

•  •  •«  3 

TAULt  OF 

IU  SI  DUALS 

CAS  t  NO. 

Y  VALUt 

Y  LSI  I  HAT  t_ 

KfcSlDUAL 

4  ut V 1  AT  108 

1 

3306  .dOUOO 

*  193.39504 

11 1.40496 

20.28644 

c 

2665.50000 

1 975.39640 

ti'0»  1036J 

*5.840*1 

2828.80000 

1475.39640 

833.40360 

30.16840 

4 

2900.00000 

1975.39640 

924. 60360 

31.88288 

5 

2194.50000 

1975.39640 

219.10360 

9.98422 

6 

1155.20000 

1210. 958.31 

-33.75831 

-4. 3*672 

7 

1094.20000 

1210.951*31 

-116.75831 

—  1 0.6 TO  06 

b 

1097.60000 

l*10.4vS31 

-113. 15831 

-10.30773 

9 

1068. 1 OCuO 

12 10 . 4583V 

-142.85631 

-13.57499 

10 

943.50000 

1210.93631 

-263.43851 

-23.07597 

11 

323.»00u0 

1210,9.3631 

-387.15631 

-46.99664 

12 

886.20000 

918.60396 

-32.46396 

-3.66323 

13 

909.60000 

918.66396 

-9.06396 

•  -0.99643 

14 

916.30000 

910.66396 

-2.36396 

-0.25749 

15 

616.500u0 

918.66396 

-102. 16398 

-12.51243 

16 

3H0. 00000 

91b. 66396 

-338.66396 

-53.39034 

17 

1731.00000 

1437.66703 

293.34297 

16.94645 

1C 

1605.70000 

1437.65703 

168.04297 

10.46340 

19 

1647.60000 

1437.65705 

110.14*97 

7.11610 

20 

1*92.60000 

1437.65703 

-145.03703 

-11. *2211 

21 

10)7.20000 

1437.63703 

-420.43703 

—41.33475 

22 

6444.80000 

4560.03583 

ibtM).  7641^ 

29.24473 

23 

5111.50000 

4560.03383 

35 1.46417 

10.76370 

24 

33^6.90000 

4J60.03383 

-1**3.13583 

-36.65485 

2  3 

1982.7C000 

2  12v*.63ol9 

- 137.95U29 

-6.95610 

26 

17b2. 70000 

2120.63619 

-337.93614 

-Id. 93766 

27 

1534.80000 

2120.63819 

—50  5.838 19 

-38.17163 

28 

1236.20000 

2120.65819 

-684.45819 

-71.54653 

29 

768.90000 

812.59923 

-43.69923 

-5.68334 

30 

376.90000 

812.399.3 

-*55.69923 

-40.85617 

31 

527.20000 

612.599*3 

-?3 5.39923 

-54.13491 

3* 

476.10060 

b  1 4.  •  j  ^^^.3 

-33 6.549*3 

-7(i.  67  b*7 

33 

529.600C0 

812.394 o 

-20 2.999*3 

-53.43641 

34 

479. <*0000 

612.54923 

—3  3  3. 19923 

-09.50339 

35 

173b. 60000 

1191.30*04 

5o7. 29796 

32.25850 

36 

17B6. 30000 

1191.30.04 

395.14796 

33.31643 

37 

1643.60000 

1191.30204 

454.29796 

*7.60683 

38 

1533.70000 

1 191.30*04 

342.39746 

22.32496 

39 

1379.66060 

1291.3020s 

188.29796 

13.64874 

4u 

1153.90000 

1191.30.04 

-37.40204 

-3.24136 

41 

lbo  3.36000 

1  490.0  b-tli. 

375.41548 

20.12419 

4^ 

1732.70600 

1490.0640* 

262 .61398 

14.98351 

43 

1605.^0000 

1490.0840* 

115.11598 

7.17144 

44 

1396.50000 

1490.0840* 

-■*3.3840* 

-6.70133 

45 

963.50006 

604.79343 

300.70655 

31.14516 

46 

828.60000 

664.79345 

164.00655 

19.78844 

47 

742.40060 

664.79345 

77.60655 

10.45347 

48 

66  1.50000 

6o^ • 

*.70635 

O.sO JS  8 

49 

618.30000 

66*.  1  **!> 

-46.49345 

-7.51956 

50 

349.60000 

66*. 7  V366 

-1 1 3 .7  4343 

-*1.09170 

51 

3317.26060 

3439.7331b 

-1*2.53318 

-3-69331 

32 

2226.90000 

3439. 75318 

-1212.83318 

-54.46>85 

33 

2746.70000 

3439. 7331o 

-6V3. 03318 

-25.23148 

54 

1186.500U0 

78b. 99490 

40  ...50510 

33.92373 

55 

1172 .80000 

783.99490 

388.S0510 

33. 16\8t 

56 

980.30000 

783.99490 

196.30510 

*0.0*500 

57 

915.20000 

733.94490 

)3l.*0510 

14.33622 

38 

863.50000 

783. 99496 

79.50510 

9.20731 

59 

60 

1744.20000 

1619.50000 

1365.16479 
1365,.  16479 

379.01521 

*54.31321 

ikimi 

61 

1606 .100  uJ 

1365.16479 

*40.91521 

15.00001 

62 

1532.40000 

1365.16479 

167.21521 

10.91 198 
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IS 
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IS 
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IS 

1.26280 
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deteriorated,  it  would  indicate  that  that  variable  was  unwanted.  The 
regression  would  then  probably  be  repeated  with  that  variable  deleted. 

The  final  F-value  is  used  to  check  the  statistical  accuracy  with  which  the 
regression  equation  represents  the  data.  Referring  to  F-tables,  for  62 
observations  and  three  independent  variables,  an  F-value  of  greater  than 
five  corresponds  to  a  better  than  1  percent  level  of  significance. 

T-values  are  similarly  evaluated.  Here  a  T-value  of  2.4  corresponds  to  a 
1  percent  level  of  significance.  Although  F  and  T  values  are  important, 
it  was  found  that  in  this  program  many  completely  unsatisfactory 
regression  results  met  high  levels  of  statistical  significance. 
Consequently,  F  and  T  values  were  not  used  to  evaluate  the  effectiveness 
of  a  particular  regression  analysis,  but  only  to  ensure  that  high  levels 
of  significance  were  being  met.  What  was  found  to  be  desirable  was  that 
the  T  values  for  each  variable  be  comparable  in  magnitude.  It,  for 
example,  one  variable  had  ten  times  the  T-value  of  another  variable,  it 
would  follow  that  the  numerical  contribution  of  the  variable  having  the 
lower  T  value  tc  the  estimated  value  of  the  dependent  variable  would  be  so 
small  as  to  render  its  presence  useless. 

The  tables  of  residuals  provide  a  good  comparison  between  different 
regression  equations.  Table  3b  shows  the  Log:L'g  equation  to  be  a  much 
more  accurate  predictor  of  the  test  data  than  the  linear  equation,  whose 
residuals  are  shown  in  Table  32. 

£1 imin ation  of  Out! iers 

Tables  of  residuals  are  often  used  to  reject  data  points  as  outliers. 

Data  points  having  the  largest  percent  deviations  are  rejected  and  the 
regression  analysis  is  repeated.  This  invariably  results  in  a  more 
"accurate"  equation.  However,  this  is  a  more  effective  procedure  when 
dealing  with  data  about  which  little  or  nothing  is  known  quantitatively, 
such  as  public  opinion  type  surveys.  In  this  program,  however,  a  great 
deal  is  known  about  the  data.  In  such  cases,  apparent  regression  outliers 
must  be  evaluated  against  the  actual  test  data.  An  illustration  of  the 
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importance  of  this  is  given  in  Figure  60.  The  four  points  a,  b,  c  and  d 
represent  response  strain  values  for  four  different  skin  thicknesses.  The 

line  drawn  (T)  -  (T)  represents  a  computed  regression  relationship 

showing  strain  increasing  with  skin  thickness.  Based  on  a  table  of 
residuals,  data  point  "d"  appeared  to  be  an  outlier.  Elimination  of 

point  "d"  resulted  in  a  new  regression  line  (T)  -  (T)  ,  which  nad 

greater  statistical  accuracy  than  the  first  line.  Since  it  is  known  that 
strain  decreases  with  increasing  thickness,  it  can  be  seen  that  the 
regression  analysis  resulted  in  an  illogical  relationship.  In  addition, 
by  removing  an  outlier  on  the  basis  of  statistical  accuracy,  the  incorrect 
trend  of  the  regression  was  worsened.  By  plotting  the  data  prior  to 
regression,  and  knowing  that  strains  decrease  with  increasing  skin 
thickness,  it  is  obvious  that  data  point  “ a 1  is  the  main  outlier  and  not 
point  "d".  When  point  a  was  removed,  the  new  regression  line  was 
( 1 )  —  (T)  ,  which  shows  a  more  reasonable  relationship.  This  example 
illustrates  the  importance  of  checking  data  for  technical  inconsistencies 
prior  to  regression  analysis,  preferably  by  graphical  means;  and  also 
demonstrates  the  danger  in  allowing  statistical  decisions  t.o  replace 
technical  ones. 


Figures  61,  62  and  63  show  graphical  representations  of  measured  rms 
strains  versus  stringer  spacing,  skin  thickness  and  radius  of  curvature 
for  each  test  sound  pressure  level.  The  only  imposition  made  on  the  data 
prior  to  regression  was  that  response  strain  must  increase  with  increasing 
sound  pressure  level,  stringer  spacing  and  radius  of  curvature;  and  must 
decrease  with  increasing  skin  thickness.  No  prior  limitation  was  placed 
on  the  rate  of  change.  This  is  a  sensible  and  practical  approach  when 

drt  1  r,  « .  i  4-  U  ^  {mm  ^  ^  rl  ^  ».  iU  ™  J  ^  4 

ca  i  my  min  iumi  luuiyuc  tc.)t  uaua*  wiici  c  ouiatcu  iliuynai  uam  put  it  Li 

are  not  unusual.  Referring  to  Figure  61,  there  were  no  inconsistencies  in 
the  data  presented.  On  Figure  62,  however,  both  panels  a  (6  ply)  and  n 
(4  ply)  had  lower  strains  than  panel  b  (8  ply).  In  order  to  determine 
whether  panel  b  response  was  too  high  or  that  for  panels  a  and  n  were  too 
low.  Figures  61  and  63  were  referred  to.  From  these  graphs,  it  does  not 

appear  as  if  the  response  for  panel  b  was  too  high.  In  fact.  Figure  63 
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Fiqure  62.  Response  Comparisons  for  Different  Skin  Thicknesses 


h  ( R=*90") 


g  (R-80") 


f  (R-30") 


Figure  63.  Response  Comparisons  for  Different  Radii  of  Curvature 
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indicates  that  at  160  dB,  the  response  for  panel  b  was  too  low.  Based  on 
this  observation,  referring  back  to  Figure  62,  the  following  data  points 
were  rejected:  panel  a  at  160,  155  and  160  dB;  and  panel  n  at  155  and 
160  dB.  Referring  now  to  Figure  63;  panel  b  at  160  dB  was  rejected,  as 
was  panel  h  at  145  and  150.  S  ne  other  data  points  were  similarly 
rejected,  based  ‘Mi  comparisons  that  did  not  lend  themselves  to  graphical 
representation.  They  were: 


panel  g  @  160  dB 

panel  i  @  165  dB 

panel  j  @  140  dB  through  165  dB 

panel  p  0  155  dB 

panel  r  @  140  dB 

panel  s  0  140  dB 

There  was  a  tendency  tor  data  at  140  dB  to  be  generally  inconsistent  due 
to  the  low  response  levels.  With  the  possible  exception  of  panel  n,  the 
data  points  rejected  did  not  interfere  with  observations  regarding  the 
nonlinearity  of  the  panel  responses. 

In  developing  a  design  method  for  both  flat  and  curved  panels,  using 
regression  analysis,  it  was  necessary  to  determine  a  numerical  value  for 
the  radius  of  curvature  of  a  lat  panel.  Since  regression  analysis  deals 
entirely  with  numerical  values,  using  a  very  high  value  to  represent  an 
infinite  radius  of  curvature  must  be  avoided.  Orginally,  a  value  of 
R  =  10,000  inches  was  used,  resulting  in  a  very  small  regression 
coefficient  for  this  variable.  Since  the  response  differ  :nces  between  the 
flat  and  the  R  =  90-inch  curved  panels  was  much  less  than  between  the 
R  =  30-inch  and  the  R  *  60-inch  panels,  it  seemed  likely  that  a  much  lower 
number  than  10,000  inches  would  be  appropriate.  Various  numbers  were 
tried,  and  also  a  graphical  review  of  the  data  was  performeo.  It  was 
determined  that  R  =  150  inches  was  a  satisfactory  number  to  represent  the 
flat  panel  radius. 
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NOTE:  The  following  important  observation  was  made  during  these 
regression  analyses.  The  signs  (+  or  -)  of  the  regression 
coefficients  of  the  independent  variables  were  determined  within 
the  computer  program,  based  on  the  sign  of  the  cross-correlation 
coefficient  between  the  independent  variables  and  the  dependent 
variable.  When  using  a  Log  function  to  represent  an  independent 
variable  having  values  of  less  than  unity,  the  program  fails  to 
take  account  of  the  logarithm  of  the  value  changing  sign.  To 
overcome  this,  the  values  of  independent  variables  were  adjusted  to 
be  always  greater  than  unity.  In  this  program,  the  skin 
thicknesses  were  entered  x  10  . 

Durbin-Watson  Statistic 

It  will  be  noticed  that  the  tables  of  reciduals  quote  the  Durbin-Watson 
statistic.  In  regression  analysis  it  is  assumed  that  the  residual  errors 
are  an  independent  random  variable.  If  the  error  terms  are  not 
independent,  but  show  serial  correlation,  this  indicates  that  the  best 
possible  curve  fit  may  not  have  been  obtained.  It  also  means  that  the  F 
and  T  tests  may  not  accurately  reflect  the  true  confidence  levels  of  the 
regression  results.  The  Durbin-Watson  statistic  is  used  to  check  for 
serial  correlation  in  the  error  terms.  In  the  regression  analyses 
performed  here,  the  Durbin-Watson  statistic  indicated  that  there  was  some 
serial  correlation  in  the  error  terms.  However,  this  is  not  critical  for 
the  regression  analysis  performed  in  this  section.  As  mentioned  earlier, 
the  emphasis  here  was  more  on  uncovering  patterns  in  the  data  than  on 
formal  tests.  The  best  possible  curve  fit  was  sacrificed  partially  in 
order  to  ensure  a  technically  acceptable  result.  Also,  the  F  and  T  values 
obtained  here  were  considerably  in  excess  of  those  values  required  to 
assure  good  statistical  accuracy,  and  therefore,  the  validity  of  these 
regression  results  is  not  affected  by  small  changes  in  their  values. 

Many  regression  analyses  were  performed  prior  to  selecting  one  for  the 
final  design  method.  Some  resulted  in  unacceptable  equations,  others 
yielded  results  that  may  be  used  as  alternatives  to  the  one  chosen  here. 
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A  major  effort  was  made  to  use  nondimensional  independent  variables,  of 
the  form  used  in  AGARD^  nomographs.  Regression  equations  do  not  usually 
balance  dimensionally,  and  therefore,  care  must  be  taken  if  using  units 
other  loan  those  used  in  developing  the  equations.  Nondimensional 
parameters  overcome  this  objection.  However,  by  imposing  certain  fixed 
relationships  between  variables,  the  accuracy  of  the  equation  usually 
suffers. 

Other  regression  analyses  performed  included  using  as  an  independent 
variable  rms  strains  calculated  from  Miles'  equation,  with  the  frequency 
component  estimated  from  the  AGARD  nomograph  and  the  static  strain 
component  computed  from  the  finite-element  models.  These  estimated  strain 
values  are  given  in  Table  2S.  The  following  equation  was  obtained: 

Log  c  =  0.33  +  0.37  Log  c  calculated  (21) 

rms  rms 

The  calculated  and  measured  rms  strains  had  a  cross-correlation 
coefficient  of  0.73  and  satisfied  the  F  and  T  tests.  However,  the  average 
error  was  50  percent  with  maximum  errors  exceeding  2:1.  Table  3b  shows 
the  residuals. 

Since  the  above  regression  yielded  unacceptable  results,  the  next  step  was 
to  regress  directly  on  the  frequency  and  static  strain  components  used  in 
the  preceding  regression.  The  calculated  frequency  had  a  cross¬ 
correlation  coefficient  of  -0.64  with  the  measured  strains.  The  computed 
static  strains  had  a  corresponding  correlation  coefficient  of  0.83.  The 
average  error  was  32  percent.  Table  37  shows  the  residuals  for  the 
following  regression  equation: 


*  (1,825.3)  1.193  0 


rms 


1  ,.001 


(22) 
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TABLE  36 


TABLE  OF  RESIDUALS  FOR  REGRESSION  OF  MILES'  EQUATION 


TABLt  OB  Rfc  SIDilALS 
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30.5  1669  40.12  3  31 

6b. 13366  82.86636 

86.51272  175.28728 

1U -19993  331.00007 

44.23738 
55.6poI 
66.954  65 
74.46248 

AVfcRAOt  4  UtVIAIiUN 

IS 

49.53363 

UUKii  1N-VT  AT  SON  ST  A1 15 TIC 

IS 

0.89991 

VON  NtUKANN*  S  RATIO 

IS 

0.91«66 
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TABLE  OF  RESIDUALS  FOR  REGRESSION  USING  CALCULATED  FREQUENCIES 

AND  STATIC  STRAINS 


TABLE  OF  RESIDUALS 


CASt  NO • 
l 

Y  VALUE 

Y  tSTXMAlt. 

RLSaUUAL 

X  DEVIATION 

33O6.9OJU0 

i9  08. 7296  3 

598.17033 

17.05696 

■i 

t 

2665.30000 

2X63.16459 

300.33341 

18.77079 

1 

2626.80000 

41.63.164  39 

663.63341 

£:$»* 

4 

2900.00000 

2X65.16459 

734.83541 

5 

2194.30000 

2165.16459 

29.33341 

1.33677 

-37.36649 

6 

1135.20000 

1586.83772 

-431.65772 

7 

1094.20000 

1586.85772 

-492.65772 
-4  89.057  72 

-45.02441 

6 

1097.60000 

1506.83772 

-44.54069 

9 

1066.10000 

1366.837/4 
1386. B 3774 

-310.73  772 

-48. 368, '.7 

XU 

S’45.  >01)00 

-641.35772 

-67.83263 

XI 

823.80000 

1386.63/72 

-763,03772 

-32.92038 

12 

768.90000 

Sol .82038 

13 

576.90000 

601 .82038 

-2  24.9  2  038 

-38.98776 

1*» 

327.10000 

801.620  3b 

-274.72038 

-52.11922 

15 

476.10000 

601.82038 

-323.72038 

19 

329.70000 

U 01 .82038 

-2  72.12038 

17 

4  79  .300  00 

801 .82038 

-3  2a  ,320  38 

-67.22010 

14 

1/36 .60000 

1144.49731 

6 14.10269 

34.91998 

19 

1766.50000 

1149.94731 

642 .00269 

20 

1643.60000 

1144.49731 

301.10269 

21 

1 3  33 . / 0000 

1144.44731 

369.20269 

22 

1379.60000 

1144.49731 

233.10269 

23 

1133.90000 

1144.44731 

9,402  69 

A:\rn 

29 

1863 .30u00 

1432.89788 

4  i.*  .602 12 

25 

1732.70000 

1432.89788 

299.80212 

17.105  16 

26 

16U3. 20000 

1432.09786 

132.302  12 

9.46803 

27 

l.-*96. >1)000 

14  32  .6  97  88 

—  36  .39  7  Uo 

26 

886.200 00 

1025.81169 

-139.61169 

29 

9u9  .60000 

1023.61169 

-116. 21169 

-iil*?il! 

30 

416.30000 

1025.81169 

-109.51169 

31 

816.50000 

1023.61169 

-209 .31X69 

-25.63523 

-76.66408 

32 

360. 000 UO 

1023.01164 

—945 ,0 1169 

33 

1/31 .OOuQO 

1618.176  39 

32.82161 

-72.47839 

3.03151 

39 

1603. 700 UO 

16  78.17839 

-4. 31382 

34 

1341.8OO00 

16  78.17839 

-l  30.37  839 

-29.82968 

3o 

1292.60000 

16  78.17839 

-385.57  8  39 

37 

1017.20000 

1678.17839 

-660.978  39 

-64.98018 

36 

463.30000 

1010.29191 

-44.79191 

-4-63923 

3» 

828 .80000 

1010.29191 

-181.49191 

-21.89816 

-36.08458 

Hi) 

74 2  .HUUO'J 

1010.24191 

“267  .89.191 

41 

6o7.300U0 

1010.29191 

-342.79191 

-31.35439 

42 

618.30000 

1010.29191 

-391.991191 

-63.39834 

h3 

3h9  .1)0000 

1010.29191 

-461.291,91 

6.28539 

-43.55107 

ii:Mh 

44 

6944.8  0000 

4790.22223 

1634.3  77  77 

45 

3111.30000 

4790.22223 

321.277)77 

46 

3336.90000 

4790.22223 

-I453.32i.23 

47 

46 

>3X7 .20000 
2226.90000 

1373 .8 0590 

1  3  7  3  .  8  03  90 

1943.39410 
833.094 10 

**9 

27h6. 70000 

1373.60590 

1372.69410 

49.9c 3 40 
21.21739 

30 

1982 .7UOOO 

1362 .02263 

420.67713 

51 

X7U2 .70000 

1362.02285 

220.67713 

12.37062 

52 

1534.80000 

1562 .02283 

-27.222  85 
-325.82285 

■i.mn 

53 

1236.20000' 

2562.02285 

-26.35681 

54 

1X86.50000' 

864 .31602 

322.18198 

27.15398 

55 

56 

X 172 .80000 
9o0 .4  OOoO 

064.31802 
864  .318  02 

308.48198 

116.08198 

it:l88 

17 

his.^uoOO 

364.31802 

30.08 1 98 

5.55966 

38 

863 .30OU0 

o64.31oU2 

—0.8 18 02 

-0.09473 

59 

60 

1744.20000 

1619.50000 

936.99664 

612.296.64 

61 

1606.  lOO  OO 
1532.40000 

607.20336 

798.69664 

62 

807.20336 

725.19664 

AVER  AO t 

4  DkVIATIUN  IS 

UURE  1N-WAISUN  STATISTIC  IS 
VUN  NEUMANN'S  RATIO  IS 


uiimi 

tvMn 

31.31944 
X. 314.02 
X.  337  60 
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This  equation  also  results  in  major  differences  from  many  of  the  test  data 
points,  and  is  not  recommended  for  design  use.  It  is,  however, 
significantly  better  than  using  Miles'  equation.  It  may  be  used  as  a 
design  guide  for  irregularly  shaped  structures  for  which  there  are  no 
simple  dimensions  such  as  b,  t  and  R;  but  for  which  computed  static 
strains  and  frequencies  are  available. 

A  regression  analysis  was  performed,  combining  c  and  f  with  b,  t  and  R  as 

o 

independent  variables.  However,  the  results  were  no  better  than  using 
only  b,  t  and  R. 

Effects  of  Nonlinear  Structural  Response  /t  s 

The  regression  analyses  described  thus  far  have  all  used  as  the 

dependent  variable.  This  imposes  the  assumption  of  linear  response. 
However,  although  some  panels  show  linear  response,  there  was  a  general 
tendency  for  strains  to  not  quite  increase  linearly  with  sound  pressure 
level.  This  does  not  necessarily  mean  that  the  structural  reponse  is 
truly  nonlinear  .  Several  test  panels  had  back-to-back  strain  gauges 
installed,  and  none  detected  any  in-plant  (membrane)  strains.  There  are, 
however,  aspects  of  son;c  fatigue  testing  that  can  give  the  appearance  of 
nonlinearity,  such  as  "clipping"  at  high  sound  pressure  levels,  overall 
if:l  end  fixture  motion  which  may  not  be  fully  responsive  to  changing 
sound  pressure  levels,  for  these  reasons,  <>pd  because  the  degree  of 
ror ! ioearity  was  not  excessive,  the  design  method  proposed  in  Section  IV. 5 
is  based  on  linear  response.  However,  in  order  to  provide  quantitative 
information  on  the  degree  of  nonlinearity,  and  thereby  make  an  alternative 
design  equa1  ion  for  those  who  may  with  to  use  it,  a  regression  analysis 
was  performed  with  overall  sound  press. ire  level  us  an  independent, 
variable.  It  should  be  v-emembereu  that  some  of  the  less  stiff  panels  were 
not  tested  at  the  high  sound  procure  levels  that  the  stiffer  panels  were 
subjected  to,  and  as  a  result,  there  will  be  a  slight  bias  in  the 
regression  equation.  The  bias,  however,  will  be  small.  Tne  following 
equation  was  obtained: 


(23) 


Log  e  =  0.3528  +  1.0458  Log  b-1.1241  Log  t 
rms  3 

+  0.4994  Log  R  +  0.873  Log  (SPL  x  10  ) 

3 

(SPL  was  entered  times  10  in  order  to  keep  its  value 
above  unity  and  prevent  the  logarithm  froii  :hanging  sign) 


1 

\ 


This  equation  had  an  average  accuracy  of  22  percent,  slightly  better  than 
the  equation  selected  for  the  design  nomograph.  The  residuals  are  shown 
in  Table  38.  Table  39  shows  the  corresponding  correlation  matrix  and 
statistics  for  the  above  equation.  A  comparison  between  Tables  35  and  38 
show  the  nonlinear  regression  to  have  significantly  lower  residuals. 
Comparing  the  cross-correlation  matrices  for  the  linear  and  nonlinear 
regressions  (Tables  33  and  39)  shows  much  lower  correlation  between  tnj)S 
and  the  panel  dimensions  (b,  t  and  R)  than  between  and  the  same 

panel  dimensions.  This  would  indicate  that  there  was  significant 
correlation  between  SPL  and  the  panel  dimensions.  The  correlation  matrix 
in  Table  39  hows  significant  correlation  between  t,  R  and  SPL.  In  fact, 
t  and  R  show  higher  correlation  with  SPL  than  with  c|ins.  This  lack  of 
independence  between  the  "independent"  variables  is  called  multi- 
col  1  inearity,  and  is  a  condition  of  deficient  data.  The  presence  of 
nuilticollinearity  does  not  necessarily  invalidate  the  regression,  but  it 
is  a  signal  for  caution.  Performing  comparative  calculations  using  the 
linear  and  nonlinear  equations  indicated  that  the  degree  of  multi- 
collinearity  was  not  excessive.  In  order  to  illustrate  the  degree  of 
nonlinearity,  the  change  in  strain  due  to  6  dB  increase  in  sound  pressure 
level  was  found  to  be  times  1.83.  The  linear  equation  would  produce  a 
ratio  of  times  2.. 

It  should  be  pointed  out  that  the  nonlinear  equation  is  the  most  accurate 
representation  of  the  test  data  in  this  program,  and  may  be  quite 
acceptable  for  design  use.  However,  since  the  degree  of  nonlinearity  was 
not  high,  and  in  the  interests  of  some  conservatism,  the  linear  equation 
was  used  to  develop  the  design  nomograph  in  Section  IV. 5. 
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TABLE  38 


TABLE  OF  RESIDUALS  FOR  NONLINEAR  REGRESSION 


TarLF  Of  RESIftU«l-5 


CASE  NO. 

Y  VALUE 

Y  ES7l"ATF 

9FSTCMIAI 

1 

101. 

70000 

69.20*9? 

n.«i5oa 

? 

77. 

7(1000 

65.09009 

12.2*09! 

7 

1*7. 

IOOOO 

1 Uh . JU6H9 

78.7931* 

* 

?99  . 

«nooo 

1 78.93097 

86.899*7 

9 

79  7. 

7  0  0  0  0 

243.90069 

9*. 079*8 

9 

73. 

90000 

*1  .2392(1 

-7.73920 

7 

96  . 

90000 

68.66230 

-1  1.76230 

8 

1  0  1  • 

n  n  0  n  0 

1  12.791 19 

-11.091 19 

9 

17*« 

1  0000 

1«6. 16927 

-12. 06927 

1  0 

?7*  • 

POOOO 

307 .96920 

-33.959Pn 

1  1 

*?5. 

1  0  n  0  0 

60q.l3<»3n 

.8* .07*70 

1? 

22. 

70000 

*4. lc  1  3? 

-6.821 7? 

1  3 

7  V  • 

00000 

*8. *0631 

-1 6. *8871 

1* 

*6. 

90000 

74.  /U9*3 

-31. 280*7 

15 

7  7. 

90000 

1  31  .*6*98 

-53.86*88 

19 

153. 

90  00  0 

217.39669 

-83. 79869 

17 

?*  {  ■ 

*0000 

369.62835 

•  1  1  ?•  1  28  74 

18 

9  1  . 

oonoo 

*  1 . 169*9 

9.87*5* 

10 

82. 

90000 

6n.539‘-7 

?*. 390*7 

20 

1  5  1  . 

*0000 

1 1 • . (o919 

76.81 085 

?1 

pc  0  . 

0000(1 

189.038,39 

9* . 19381 

2? 

*0,0. 

10000 

307. 30873 

92.70127 

23 

595. 

*0000 

908. 22395 

87.1780* 

?* 

9*. 

1  0000 

9  0 .*0*09 

3.89501 

25 

?p6. 

70  00  0 

227.5*300 

5«.1570o 

?6 

*95. 

5  0  0  0  0 

379.27999 

P9.P270* 

27 

7?0  . 

9onoo 

f/V.  90287 

98,31717 

PH 

?9. 

70  0  0  0 

31,609** 

•9,809** 

?9 

*  f. 

.700*0 

5p.«0?*8 

-6,182*8 

70 

n*. 

70000 

86.33?6* 

—  2 • 0  326* 

71 

1  73, 

1  0000 

1*?. 2*5*9 

-9. 1*5*6 

3? 

1  9H. 

20000 

239.22*91 

-67.02*51 

77 

9U. 

20000 

*8. 1*9*2 

2.05058 

3* 

8  3. 

90  non 

80. 16768 

3.77275 

79 

1  *9. 

*0*00 

131 .4***9 

1  0  .  *  7  9  5 1 

79 

?1  0. 

70000 

217. 36*60 

-6.66*90 

77 

pq*'. 

o  0  0 11 0 

364. **s*fc 

-6* . A  *8*  * 

88 

p«. 

onnno 

23. 95p19 

*.3*781 

79 

*  3. 

loono 

3 9.3603* 

3.71098 

*0 

9**. 

70000 

6-., 60*97 

3, *95*5 

*1 

1  0«- 

80000 

109.7  7*67 

2.02513 

*2 

170. 

70000 

179.60850 

*.771^0 

*7 

PH  J« 

70000 

29?. 007  1  2 

-8.7071 3 

** 

1  «9. 

90000 

1*1.76310 

*6,11890 

*9 

P96 

80000 

238.0688? 

?9„73**P 

*9 

30?. 

oonoo 

38 7. *7117 

-Hi. *7117 

*7 

96. 

20000 

109.60064 

-1 3 ■ 66  080 

*8 

1 1  5. 

nnnoo 

182.415*7 

-6*.l 15*7 

*9 

P9  2 

70000 

301.00705 

-*  8.30708 

50 

57. 

sonno 

98.6779* 

-1 1.1779* 

91 

9  2 

7oono 

1  1  * . 3*  71 9 

_? 1  .8*71  * 

9? 

1*1. 

POOOO 

168.1 7093 

-*9,97057 

97 

?0  1 

90000 

31  (I  .  03608 

-1 08 .53008 

9* 

9  >• 

Too  00 

*9.07*66 

15.8257* 

99 

1  "  7 

90  0  00 

79. *4169 

7?. *0871 

59 

150, 

80000 

12*. 363*3 

35. *1657 

57 

26  5 

*0000 

206.66943 

.*S‘?129* 

“1  77 

**6 

n  0  ij  *i  D 

Da  •  ’  '  C  r  ” 

lnj««  1  M  > 

%  DFVIATION 
ll. 61758 
15.8*72* 

?9. 77195 
33.1064* 
l7.QO.-31? 

:$  o’.UUl 

-1  i  .872*6 
s6. 93238 
-12.275*2 
T1 9.T6613 
_3D.5«8H7 
-6l „9?103 
,6*.5\*29 

-A  1  ,<<3*30 

-*5.3?i*4«» 

1 9.383*0 
?6. 22225 
?5 .5026* 

?5 . A 96*5 
S>3.1620? 
1*.64}S9 
6.83163 
20.35547 
19.1971* 
>3.6*379 
-?2. 60*63 
-10.9143? 

.zhxm 

-39.84«10 
*. 08*81 

3.99083 

7.796*0 

-3.16J08 

l:t?m 

1.8*133 
|.9?3*7 
-3, 073*6 
?*. 1*999 
1 1 .1 0679 
-?6.83T«* 
-l*.?OOSl 
-87.9*1600 
-19.11WA 
Tl9,*7990 
-•>3. 78184 
133.26525 
-S3.«At>04 
2*. 6*8194 

IWIBI 

la  ^  <  .  lN4 

?  10-3 


99 

90 


91 


9? 


o'j./nooo 

1  *  ^  .  n  0  n  (((I 

r6l.n0000 

***.*0000 


{~.t£ 380 

126.  '25*7 

21 5 .04*09 

350.  '2465 


1 2.*76?o 
?0.?7*33 
*9. 70801 
93.67018 


8VEP*r.r  *  nFVJATYoil  T9  21.61903 


r>U>  *>  iN-WATSON  97 A T T 5Y r C  T9  0.86H77 


VO'.  MFIMANH*?  o  ATTf>  18  0.88301 
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TABLE  39 


REGRESSION  ANALYSIS  RESULTS  FOR  NONLINEAR  REGRESSION 


couffi  attpN  ma1*TX 


00* 

b 

t 

R 

SPL. 

1  rms 

1  , 0  0  n  0  0 

n.?3YS7 

-0 . 14990 

-0.03740 

0.14100 

RO* 

? 

0 , ?\Sc J 

l .noooo 

fl  ,  1SH4A 

0 . 1 A7?1 

-0 . 1 *0?A 

00* 

7 

0 . 1 

0 . \ SH4  A 

1 .uuflim 

-n , Hon  I  0 

n . n?S07 

DOW 

4 

0  .  (1774Y 

0. 1«7?I 

-D.FtiOlO 

1 .OOnop 

o.tO07*i 

BOW 

S 

n,n;ii  '> 

-n. 1 ho?a 

0 .  u<fHf)7 

I1.7007F 

1  ,  OnOOO 

V»WlAO|F  fNlrpH',,,,.  7 

SllH  Of-  sOiiAhf*;  PFOPCFO  TN  THIS  STFR....  0„«7n 

ppOhootton  npnm.rn  in  tm7S  5.UH. ,  o.m* 

riiHiii  srrvf  Sum  nr  souanrs  Hti-ut-c.fi .  7.477 

■"umih  attvf  Hpr.HuoT T tin  . .  0.014  of  a.ONS 

for  •  v4nl‘'f«  t.s  FMTf«rri 

MUITJPLF  CopHtLATlON  eOEF  F  IC  itNT  ,  .  .  il.OSf. 

(AoJ'isiH/  FOB  D.F . . .  .  O.RS* 

F-VAjUF  FUp  ANALYSIS  OF  VAWIANCt  .  .  .  IM.KB* 

STANnfiOn  LrpUp  or  FSttmATF..-* .  o.nn 

(AdJuctHi  fop  f'.t,) . .  n.117 

VAP5AOIF  PFOPKSSTON  b'D,  fHl-tUP  OF  rOMPIITFn 

NUMPFP  (.OFFFIl  TFNT  '<EG,  COEFf.  t-waiiif. 

4  “.A73041«H  0.0370?  .Z3.**? 

?  -]  .  1?40a?7?  0.0S4AS  -11. *7* 

1  l.n4S7*r,A3  0.1ZZ«S  »,SI  ■» 

0  n.AC'tfcnm  0.0SHO7  »,4t7 

INTFor.’-'P1  n  .  as?  7os  S  7 
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Regression  analyses  were  also  performed  on  the  test  frequencies  and  the 
static  strains  computed  from  the  finite-element  models;  against  the  panel 
dimensions  b,  t  and  R.  It  was  originally  hoped  that  simplified  regression 
equations  could  be  developed  for  frequencies  and  static  strains,  which  in 
turn  could  be  used  as  inputs  to  a  regression  equation  for  rms  strain. 
However,  the  regression  equation  for  rms  strain  which  utilized  frequencies 
and  static  strains  as  inputs  was  not  particularly  accurate.  In  addition, 
it  was  not  possible  to  obtain  a  satisfactory  regression  equation  for 
static  strains  on  the  curved  panels.  For  the  flat  panels  only,  however,  a 
regression  equation  between  the  computed  static  strains  and  the  panel 
dimensions  produced  an  equation  with  an  average  accuracy  of  5  percent: 

t  ,  C921)(L.4105b) 

0  (lO18)* 


Regression  analysis  on  the  test  frequencies  against  the  panel  dimensions 
gave  satisfactory  results,  but  since  the  existing  AGARD  nomograph  also 
gave  satisfactory  results,  there  was  no  point  in  using  a  regression 
equation  in  the  design  method. 

4.  DESIGN  EQUATIONS 

Before  proceeding  with  the  final  design  method,  it  seems  useful  to 
fully  describe  the  regression  equation  used,  and  to  present  those 
alternatives  that  are  /iable. 

The  basic  design  method  utilizes  overall  sound  pressure  levels  rather  than 
spectrum  levels  as  applied  loads.  Consequently,  it  is  not  necessary  to 
determine  the  natural  frequency  in  order  to  determine  rms  strain. 

Response  frequency  estimates,  however,  are  required  in  order  to  estimate 
fatigue  lives. 
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a.  Calculation  of  RMS  Strains  -  Acceptable  equations  developed  in 
Section  IV. 3  follow: 

(i)  Equation  for  rms  strain,  used  to  develop  the  design 
nomograph  in  Section  IV. 5  is: 

£ 

L°g  -J-J  =  3.080614  +  1.104533  Log  b  (20) 

-1.206903  Log  (103.t)  +  0.551923  Log  R 

where  erms  is  in  microinches/ inch 
SPL  is  in  lb/in^ 
b  is  in  inches 
t  is  in  inches 
R  is  in  inches 

and  where  R  =  150  inches  for  a  flat  panel,  i.o.,  for 
150  <  R  <  ,  use  150. 

NOTE:  (1)  Because  the  equation  is  in  Log:Log  form,  six 
significant  figures  are  required  in  order  to 
maintain  reasonable  accuracy. 

(2)  This  equation  assumes  linear  response,  i.e., 
6  dB  represents  a  doubling  of  strain. 

(3)  Strain  levels  were  quoted  rather  than  stress 
levels  in  order  to  allow  elastic  modulus 
values  for  alternative  skin  laminates  to  be 
used  in  determining  stress  levels. 

(4)  Ail  logarithms  are  to  the  base  10. 
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This  equation  may  be  written  in  the  form: 


(8 


.3234  x  ID'10)  10<SPL'20>  b1-1045  R0-5519 


'rms 


,1.2069 


(25) 


where  SPL  is  in  decibels  (dB). 

(ii)  Nonlinear  Response  Equation  (from  equation  23): 

-Q-  ,1.045b  .0.4994  in0. 04365. SPL 

(1.394  x  10  _J_b _ K _ 10 _ 

is  =  "  tl .  1241  ~ 

where  the  variables  are  defined  as  in  (i)  above. 

NOTE:  This  equation  is  slightly  more  accurate  than  (i), 
but  may  be  slightly  unconservative  at  very  high 
sound  pressure  levels 


Effects  of  Aspect  Ratio 

In  this  program,  the  long  side  of  the  individual  bays  was  kept  constant  at 
■  =  12  inches.  The  short  side  varied  from  4-inches  to  8-inches,  thus  the 
aspect  ratio  varied  from  1.5  to  3.  In  order  to  include  the  effects  of 
varying  aspect  ratio,  the  stress  nomograph  in  Reference  5  was  used.  RMS 
stresses  for  various  aspect  ratios  were  calculated,  and  the  ratio  effects 
applied  to  the  strains  calculated  in  this  program.  Within  the  levels  of 
accuracy  of  the  AGARD  stress  nomograph,  and  also  considering  the 
accuracies  of  the  regression  equations;  there  were  no  significant  changes 
in  strain  response  at  aspect  ratios  above  1.5.  Based  on  the  combined 
effects  of  aspect  ratio  on  static  strains  and  frequencies,  the  dynamic 
strains  at  a/b  =  1.5  and  a/b  =  •*>  are  within  approximately  5  percent  of 
each  other.  Consequently,  the  design  nomograph  was  constructed  assuming  a 
common  response  for  all  aspect  ratios  above  1.5.  Aspect  ratio  lines  for 
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a/b  -  1.2  arid  1.0  were  then  superimposed  based  on  stress  ratios  determined 
from  the  AGARD  rnis  stress  nomograph. 

Effects  of  J  Stiffeners 

The  static  analysis  showed  that  J  stiffeners  resulted  in  significantly 
lower  edge  strains,  compared  to  using  Z  stiffeners.  The  ratio  of  J  to  1 
being  0.71.  The  measured  data  gave  a  corresponding  ratio  of  0.7  to  0.8 
for  configuration  b.  Configuration  c  did  not  show  any  strain  reduction 
due  to  the  J  stiffeners,  but  it  did  show  a  major  increase  in  fatigue  life. 
Based  on  these  observations,  it  is  recommended  that  the  calculated  strains 
in  this  design  method  be  factored  by  0.8  when  using  J  stiffeners.  This  is 
believed  to  be  a  slightly  conservative  factor. 

b.  Calculation  of  Natural  Frequencies  -  Reference  5  was  used  to 
calculate  the  fundamental  fully-fixed  panel  frequencies.  The  equations 
and  nomograph  are  included  in  Section  IV. 5. 

This  method  was  derived  for  typical  metals,  and  assumes  typical  values  for 
elastic  modulus  (E)  and  material  density  (p).  It  is  necessary  to  modify 
these  values  when  using  composite  materials.  Graphite/ex.poxy  laminates 
have  a  density  of  0.055  lb/in^.  The  elastic  modulus  varies  with  ply 
orientation.  In  this  program,  elastic  modulus  values  are  given  in 
Tables  2  through  6. 

Reference  8  is  recommended  for  use  if  stiffener  properties  are  to  be  used, 
rather  than  assuming  fully-fixed  edges.  For  flat  panels  only,  Reference  5 
contains  a  simplified  '.net hod  based  on  Reference  8.  However,  both  the  flat 
and  curved  panel  equations  in  Reference  8  lend  themselves  to  programming 
on  modern  desk  top  computers. 


5.  DESIGN  METHOD 

This  section  contains  a  semi-empirical  method  for  estimating  rms 
strains  and  natural  frequencies  for  curved  and  flat  graphite-epoxy  skin- 
stringer  panels,  in  order  to  predict  their  sonic  fatigue  lives.  A  random 
rms  strain  versus  cycles  to  failure  curve  is  presented  for  bonded  skin- 
stiffener  joints.  A  worked  example  is  also  presented. 

The  design  equati  >n  and  corresponding  nomograph  for  rms  strain  was  based 
on  Z  stiffeners.  However,  they  can  be  readily  factored  to  allow  for  the 
use  of  J  stiffeners. 

a.  Estimation  of  RMS  Strain  -  The  RMS  strain  nomograph  is  shown 
in  Figure  64.  It  is  based  on  equation  25  from  Paragraph  IV. 4. a: 


_  (8.3234  x  10~1Q)  10^-SPL/2Q)  b1'1045  r0-5519 


,1.2069 


where  t  =  Maximum  rms  strain  at  panel  edge  due  to  random 
rms  ,  -6 

acoustic  loading  (10  in/in) 

SPL  --  Overall  sound  pressure  level  (dB) 

a  =  Panel  length,  between  longerons  (inches) 

b  =  Panel  width,  between  stringers  (inches) 

t  «  Skin  laminate  thickness  (inches)  (also  given  on  nomograph 

as  number  of  plies) 

R  =  Radius  of  curvature  in  !lb"  direction  (inches). 

R  =  150  in.  is  the  maximum  value  to  be  used  in  the 
response  equation,  and  is  valid  for  flat  panels  and 
all  R  >  150  in. 


Equation  is  valid  for  a/b  >1.5. 

For  a/b  =  1.2,  factor  equation  by  x  0.849. 
For  a/b  =  1.0,  factor  equation  by  x  0.744. 
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b.  Estimation  of  Fully-Fixed  Natural  Frequency 


From  Reference  5, 


f  =  V.K.—  (27) 

b2 

where  f  =  frequency  of  fundamental  fully-fixed  natural  frequency  (Hz). 


V  =  (Ey/p)1/2/200,000 


(28) 


where 


E  is  the  elastic  modulus  in  the  "y"  direction,  i.e., 

"la"  direction  (lb/in  ).  Obtained  from  Tables  2 
through  6  for  laminates  used  in  this  program.  For  other 
laminates.  Reference  12  may  be  used. 


o  =  density  of  skin  laminate.  If  expressed  in  units  of 
Ib/in  ,  it  must  be  divided  by  386.4.  For  graphite- 
epoxy  laminates, 


p  = 


K  is  obtained  from  Figure  65  for  given  b,  t  and  R. 


b,t  and  R  are  defined  as  in  Paragraph  IV. 5. a,  except 
that  true  values  are  lo  be  used  for  all  R,  including  °° 
for  flat  panels. 


c.  Estimation  of  Fatigue  Life  -  Estimated  sonic  fatigue  life  is 
obtained  by  reading  number  of  cycles  to  failure  (N)  from  Figure  66  for  rms 
strain  (‘"rms)  calculated  in  Paragraph  IV. 5. a.  The  number  of  cycles  to 
failure  is  converted  to  life  in  hours  by  the  relationship  . 
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Figure  66.  Random  Fatigue  Curve  for  Bonded  Skin-Stiffener  Joint 
RMS  Strain  vs.  Cycles  to  Failure 


d.  Calculation  Procedure 

(1)  Estimate  nns  strain  (e^  )  from  Figure  64  or  from 
equation  25,  for  given  SPL  (dB),  b,  t  and  R,  (For  flat 
panels  and  other  R  >  150,  put  R  =  150-in.). 

If  J  stiffeners  are  used,  multiply  £rms  by  0.8. 

(2)  Calculate  V  from  equation  28. 

(3)  Estimate  K  from  Figure  65  for  given  b,  t  and  R.  Use 
R  «  for  flat  panels. 

(4)  Calculate  frequency  (f)  from  equation  27. 

(5)  Estimate  number  of  cycles  to  failure  (N)  from  Figure  66 
for  estimated  erms  from  (1). 

(6)  Convert  cycles  to  failure  to  fatigue  life  in  hours,  using 
calculated  frequency  from  (4). 

e.  Worked  Example  -  A  curved  panel  is  made  up  from  an  8  ply  skin 
laminate,  having  a  ply  orientation  of  (0,-  45,  90)s,  which  corresponds  to 
t  =  .044  in.  and  Ey  =  6.7  x  10^  lb/in. The  panel  has  Z  section 
stiffeners,  with  bay  dimensions  of  a  =  12  in.  and  b  =  8  in.  The  radius  of 
curvature  is  60  in.  The  overall  sound  pressure  level  is  165  dB. 

(i)  From  Figure  64,  the  corresponding  rim.  strain  <;rms  -  615 
microinches/ inch. 

(ii)  V  =  (Ey/o)1/2/2Q0,00Q,  p  =  .000142 

1/2 

7200,000  =  1.0861 


*-'«»’** . 


(iii)  In  order  to  determine  K  from  Figure  65,  first  calculate: 
a/b  =  -  1.5 


bl 

Rt 


=  24 


From  Figure  65,  K  =  0.68  x  10  . 


( iv)  From  equation  27 : 


f  =  V.K. 


7 


f  =  (1 .0861 ) (0.68  x  10°) (.044) 

82 


=  508  Hz 

(v)  from  Figure  66, 


Cycles  to  failure  N  =  8  x  105 


for  eras  -  615. 


(vi)  Fatigue  life  =  jtztqTT 

=  0.44  hours 


NOTE:  The  estimated  rms  strain  of  615  p  in. /in.  corresponds  to 
a  measured  value  of  595  p  in-./in.  The  estimated 
frequency  of  508  Mz  corresponds  to  a  measured  value  of 
350  Hz. 


f.  Fan  Noise  as  Acoustic  Load  -  The  design  method  presented  here 
utilizes  the  overall  sound  pressure  level  as  the  design  load.  While  this 
is  adequate  for  the  broad-band  spectra  that  typify  jet  exhaust  noise,  it 
does  net  automatical ly  lend  itself  to  designing  for  the  inlet  and  fan  exit 
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acoustic  spectra  that  occur  due  to  fan  noise  on  high  bypass  ratio  engines. 
These  spectra  often  have  overall  sound  pressure  levels  that  are  dominated 
by  very  high  acoustic  spectrum  levels  occurring  at  the  blade  passage 
frequency  and  some  of  its  harmonics.  Since  the  blade  passage  frequency 
(typically  2-4  KHz)  is  usually  well  above  the  frequency  range  of  interest 
for  structural  reponse  (typically  50-1,000  Hz),  including  it  in  the 
overall  sound  pressure  level  to  be  used  in  this  design  method  may  result 
in  overly  conservative  designs.  In  order  to  deal  with  this  type  of 
acoustic  load  it  is  necessary  to  eliminate  these  high  frequency  peaks,  and 
develop  an  estimate  for  the  overall  sound  pressure  level  from 
50  Hz-1,000  Hz.  This  can  be  accomplished  in  one  of  the  following  ways. 

(1)  If  the  actual  measured  acoustic  data  are  available  (e.g.,  on  magnetic 
tape),  reanalyze  the  data  to  measure  the  overall  sound  pressure  level  with 
a  1,000  Hz  cutoff  filter  applied.  (2)  If  only  the  acoustic  spectrum  plot 
is  available,  then  sum  together  in  sequence  all  the  significant  peaks 
between  50  Hz  and  1,000  Hz  in  the  following  way: 


If  the  difference  (in  dB)  between  two  peaks  is  (any  bandwidth): 


0 

i  . 

2 

3 

B 

6 

7 

8 

9 

10 

- 1 

>10 

then  add  to 
the  larger 
peak  (dB) 

3 

_ 

1.8 

D 

1.0 

J 

0.8 

.. 

0.6 

0.5 

0.4 

0 

Example:  If  there  were  four  peaks  with  the  following  levels  (dB) : 


140,  134, 

\  / 

difference  =  6  dB 
add  140  +  1.0 
=  141dB 

\ 

difference 


139  and  139,  then 

\  / 

difference  ::  0  dB 
then  add  139  +  3 
=  142  dB 

/ 

1  dB 


add  142  +  2.5 
=  144.5  dB 


This  144.5  dB  is  the  overall  sound  pressure  level  to  be  used 
as  the  design  load. 
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6.  DESIGN  METHOD  COMPARISONS 

Graphite-epoxy  skin-stringer  structures,  of  the  type  evaluated  in  this 
program,  are  primarily  in  competition  with  similarly  configured  aluminum 
structures  for  application  on  both  military  and  civil  aircraft. 

Cost/weight  tradeoffs  between  graphite  and  aluminum  structures  having 
comparable  sonic  fatigue  resistance  are,  therefore,  of  interest  to 
potential  users  of  the  design  method  in  Section  IV. 5.  In  order  to  provide 
an  estimate  of  the  sonic  fatigue  resistance  of  graphite  relative  to 
aluminum,  comparisons  were  made  between  the  method  in  Section  IV. 5  and 
those  in  References  2  and  5.  One  difficulty  in  making  these  direct 
comparisons  is  that  References  2  and  5  utilize  acoustic  spectrum  levels, 
whereas  Section  IV. 5  here  utilizes  overall  sound  pressure  levels.  In 
order  to  overcome  this,  a  typical  broad  band  sonic  fatigue  design  load 
spectrum  was  used  which  had  an  overall  sound  pressure  level  of  157  dB  and 
a  corresponding  acoustic  spectrum  level  of  132  dB/Hz  in  the  frequency 
range  of  interest.  This  25  dB  difference  between  overall  and  spectrum 
levels  is  compatible  with  the  acoustic  load  in  this  program. 

The  following  example  problem  was  used  for  comparison: 

Required  life  =  1(/  cycles  (using  the  -50%  confidence  level 

from  Reference  2). 


b  =  8 
a/b  =  2 

C  =  .02 

Calculate  required  skin  thicknesses. 
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The  following  results  were  obtained: 

Reference  2  -  Riveted  Aluminum  Skin-Stringer:  t  =  .05-in. 

Reference  5  -  Riveted  Aluminum  Skin-Stringer:  t  -•  .076-in. 

Section  IV. 5  -  Bonded  Graphite  Skin-Stringer:  t  =  .041-in. 

Paragraph  IV. 4. a.  -  Using  the  nonlinear 

response  equation:  t  =  .Q37-in. 

Comparisons  of  sonic  fatigue  resistance  using  References  2  and  5  yield 
similar  results  for  some  configurations  and  very  different  results  for 
others.  Where  differences  ao  occur.  Reference  5  is  the  more  conservative. 

Since  the  density  of  graphite  is  approximately  half  that  of  aluminum,  the 
potential  weight  saving  of  graphite  over  aluminum  for  equivalent  sonic 
fatigue  resistance  is  significantly  large.  Comparing  values  generated  in 
this  program  with  those  in  Reference  2,  the  weight  of  graphite  is  slightly 
less  than  half  that  of  the  equivalent  aluminum  structure. 
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SECTION  V 
CONCLUSIONS 


A  satisfactory  sonic  fatigue  design  method  was  developed  for  curved  and 
flat  graphite-epoxy  skin-stringer  panels.  The  design  method  is  presented 
as  a  self-contained  section  suitable  for  application  to  aircraft 
structural  design.  Design  trade-offs  with  aluminum  structures  indicated 
that  graphite  offers  a  2:1  weight  saving  over  aluminum,  for  comparable 
sonic  fatigue  resistance. 

Analytical  results  indicated  that  ply  stacking  order  may  have  an  effect  on 
sonic  fatigue  life.  However,  the  effects  were  not  quantified  sufficiently 
to  facilitate  the  inclusion  of  this  variable  in  the  design  method.  This 
is  an  area  requiring  further  work. 

The  finite-element  analyses  gave  good  static  strain  distributions  for  the 
test  panel  configurations.  These  computed  strains  displayed  high 
statistical  correlation  with  measured  dynamic  strains.  Element  grid  size 
was  found  to  be  critically  important  at  panel  edges.  Representing  panel 
stiffeners  as  plates  in  three-dimensional  models  resulted  in  more  accurate 
computed  strains  than  when  representing  stiffeners  as  beams  in  two- 
dimensional  models. 

The  finite-element  analyses  did  not  result  in  satisfactory  frequency 
estimates  for  all  the  t.  st  panel  configurations.  On  the  stiffer  panels, 
mode  shapes  were  dominated  by  motion  of  the  substructure,  resulting  in  low 
frequency  estimates.  Test  data  contradicted  this  response  behavior.  The 
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best  analytical  comparisons  with  measured  frequencies  were  obtained  using 
simple  frequency  calculations  assuming  fully-fixed  edges. 

Adhesive  bonding  problems  encountered  during  the  early  stages  of  test 
specimen  fabrication  demonstrated  the  importance  of  assuring  good  bond  and 
laminate  quality  prior  to  fabricating  expensive  sonic  fatigue  test  panels. 
Shaker  testing  of  small  coupons  was  found  to  be  an  excellent  method  of 
evaluating  specimen  quality.  The  modes  of  failure  during  shaker  testing 
showed  good  correlation  with  progressive-wave  tube  failures. 


i 
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The  sonic  fatigue  data  obtained  during  the  progressive-wave  tube  tests 
showed  good  correlation  with  variations  in  panel  configuration  parameters. 
The  data  also  displayed  some  inconsistencies,  characteristic  of  sonic 
fatigue  testing.  Sonic  fatigue  failures  were  generally  observed  to  occur 
over  long  periods  of  time,  often  over  several  million  cycles.  The  first 
signs  of  fatigue  damage  were  fractured  skin  laminate  fibers  in  the 
stiffener-skin  joint  areas.  Isolated  fiber  failures  would  continue  to 
occur  until  the  extent  of  the  skin  damage  resulted  in  separation  from  the 
stiffeners.  This  type  of  slow  progressive  failure  presents  problems  in 
defining  time  to  failure.  The  mode  of  sonic  fatigue  failure  was 
indicative  of  good  quality  structural  panels.  The  fatigue  lives  of  bonded 
and  riveted  joints  were  compared.  Riveted  joints  displayed  slower  rates  of 
progressive  damage. 

Slight  variations  in  ply  orientation  did  not  appear  to  affect  panel 
response  or  fatigue  life.  However,  it  is  expected  that  major  variations 
in  ply  orientation  would  have  significant  effects. 

J  stiffeners  were  observed  to  result  in  significantly  longer  fatigue  lives 
than  did  Z  stiffeners. 

Panel  responses  showed  some  degree  of  nonlinearity,  however,  back-to-back 
strain  measurement  did  not  reveal  any  membrane  strains. 
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Multiple  stepwise  regression  analysis,  relating  rms  strains  to  panel 
configuration  parameters,  provided  the  basis  for  the  recommended  design 
method.  Miles'  equation  del  not  show  good  correlation  with  the  test  data. 
Linear  and  nonlinear  equations  were  developed  to  predict  panel  response, 
tmphasis  was  placed  on  the  importance  of  critically  reviewing  the  test 
data  prior  to  regression  analysis.  The  potential  hazards  of  using 
regression  analysis  were  discussed  in  some  detail.  A  design  nomograph  to 
predict  rms  strains  was  constructed  and  presented  as  part  of  a  sonic  fatigue 
design  method.  A  worked  example  was  also  presented. 
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APPENDIX  A 


FINITE  ELEMENT  MODELS  AND  RESULTS 


Figures  A-l  through  A-32  are  computer  generated  plots  of  the  following 

•  Finite  element  models  used  in  the  static  and  dynamic  analyses 

•  Static  deflections  and  stresses  for  panels  b,  d  and  f 
v  Mode  shapes  for  panels  b  and  f 


A-l 
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Figure  A-6.  Finite  Element  Model 
of  4  x  3  Center  Panel 
(Radius  of  Curvature  =  90  In.) 
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Figure  A-ll.  Static  Deformation  for  Configuration  d 
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Figure  A- 12 .  Z-Displacement.  for  Configuration  d 
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Figure  A-13.  Normal  Stress  (y)  for  Configuration  d 
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Figure  A- 14 .  Normal  Stress  (x)  for  Configuration  d 
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Figure  A-16.  Z-Di spl acement  for  Configuration  f 
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Figure  A-18.  Normal  Stress  (x)  for  Configuration  f 


Figure  A-27.  Third  Harmonic  for  Panel  b  (f 


Figure  A-28.  Fourth  Harmonic  for  Panel  b  (r 


Figure  A-29.  First  (Fundamental)  Harmonic  for  Panel  f  (f 
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APPENDIX  B 


ENGINEERING  DRAWINGS  OF  TEST  STRUCTURES 

Figures  B-l  through  B-3  are  engineering  drawings  for  the  sonic  fatigue 
test  panels  and  the  shaker  test  specimens.  The  test  panel  drawings  show 
strain  gauge  locations  and  their  corresponding  nunhers.  Hercules  3501 
graphite/epoxy  system  was  used  for  the  laminates.  Stiffener-skin  bonds 
utilized  the  3M  AF147  adhesive. 
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Figure  B-l.  Fatigue  Test  Panels  (Sheet  1  of  3)  -  Continued 
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Figure  B-l.  Fatigue  Test  Panels  (Sheet  2  of  3)  -  Continued 
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Fatigue  Test  Panels  (Sheet  3  of  3) 
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Figure  B-l.  Fatigue  Test.  Panels  (Sheet  3  of  3)  -  Concluded 
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APPENDIX  C 

TEST  DATA  USED  IN  THE  DEVELOPMENT  OF  THE  DESIGN  METHOD 

This  section  contains  acoustic  and  strain  spectra  and  associated  overall 
rms  levels  that  are  pertinent  to  the  development  of  the  design  method. 
Figures  C-l  through  C-7  show  microphone  spectra  for  a  sine  sweep  and  for 
each  test  acoustic  load  level.  Figures  C-8  through  C- 101  show  strain 
spectra  and  their  corresponding  overall  rms  strain  levels  for  strain 
gauge  number  10  (located  at  the  center  of,  and  normal  to,  the  longest 
side  of  the  center  bay  —  see  Figure  B-l).  Spectra  are  shown  for  sine 
sweeps  and  broadband  random  acoustic  loading  from  140  dB  up  to  (in  most 
cases)  165  dB  in  5  dB  increments  for  the  following  panel  configurations: 
al,  b2,  d,  f2,  g2,  h,  i,  j,  kl,  1,  n,  o,  p,  q,  r  and  s. 

All  spectra  have  an  effective  filter  bandwidth  of  2-16  Hz  with  an 
aliasing  filter  cutoff  of  960  Hz.  Twenty-second  samples  were  used  for 
all  spectral  analysis. 
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Figure  C-3.  Microphone  Spectrum  -  145  dB  Random 
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Figure  C-4.  Microphone  Spectrum  -  150  dB  Random 
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Figure  C-5.  Microphone  Spectrum  -  155  dB  Random 
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Figure  C-6.  Microphone  Spectrum  -  160  dB  Random 
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Figure  C-8.  Strain  Spectrum  for  Panel  al 
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Figure  C-10.  Strain  Spectrum  for  Panel  al 
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Figure  C- 1 3.  Strain  Spectrum  for  Panel  b2 
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Figure  C-14.  Strain  Spectrum  for  Panel  b2 
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Figure  C  15.  Strain  Spectrum  for  Panel  bZ 
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Figure  C-17.  Strain  Spectrum  for  Panel  b2 
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Figure  C-18.  Strain  Spectrum  for  Panel  d 
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Figure  C-19.  Strain  Spectrum  for  Panel  d 


1000 


400 


n'vtr*. 


PANEL  CONFIGURATION: 
TRANSDUCER: 

OVERALL  R.M.S.  LEVEL: 
INPUT  SPECTRUM: 

INPUT  LEVEL: 


f2 

G10 
22.3  lie 
RANDOM 
140  dB 


Figure  C-25.  Strain  Spectrum  for  Panel  f2 
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Figure  C-29.  Strain  Spectrum  for  Panel  f2 
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Figure  C-30.  Strain  Spectrum  for  Panel  f2 
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Figure  C-3 2.  Strain  Spectrum  for  Panel  g2 
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Figure  C-34.  Strain  Spectrum  for  Panel  g2 
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Figure  C-35.  Strain  Spectrum  for  Par.el  gZ 
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Figure  C-37.  Strain  Spectrum  for  Panel  g2 
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Figure  C-38.  Strain  Spectrum  for  Panel  h 
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Figure  C-42.  Strain  Spectrum  for  Panel  h 
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Figure  C-43.  Strain  Spectrum  for  Panel  h 
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Figure  C-47.  Strain  Spectrum  for  Panel  i 
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Figure  C-49.  Strain  Spectrum  for  Panel  i 


C-50 


1000: 


PANEL  CONFIGURATION: 
TRANSDUCER: 

OVERALL  R.M.S.  LEVEL: 
INPUT  SPECTRUM: 

INPUT  LEVEL: 


i 

G10 

168.2  ue 
RANDOM 
160  dB 


FREQUENCY  (Hz) 


Figure  C-50.  Strain  Spectrum  for  Panel  i 
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Figure  C- 51 .  Strain  Spectrum  for  Panel  i 
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Figure  C-52.  Strain  Spectrum  for  Panel  j 
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Figure  C-54,  Strain  Spectrum  for  Panel  j 
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Figure  C-55.  Strain  Spectrum  for  Panel  j 
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Figure  C-56.  Strain  Spectrum  for  Panel  j 
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Figure  C-57.  Strain  Spectrum  for  Panel  j 
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Figure  C-58.  Strain  Spectrum  for  Panel  j 
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Figure  C-62.  Stiain  Spectrum  for  Panel  kl 
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Figure  C-64.  Strain  Spectrum  for  Panel  kl 
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Figure  C-66.  Strain  Spectrum  for  Panel  1 
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Figure  C-67.  Strain  Spectrum  for  Panel  1 
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Figure  C-68.  Strain  Spectrum  for  Panel  1 
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Figure  C-71.  Strain  Spectrum  for  Panel  1 
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Figure  C-74.  Strain  Spectrum  for  Panel  n 
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Figure  C - 7 6 .  Strain  Spectrum  for  Panel  n 
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Figure  C-83.  Strain  Spectrum  for  Panel  q 
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Figure  C-84.  Strain  Spectrum  for  Panel  q 
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Figure  C-85.  Strain  Spectrum  for  Panel  q 
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Figure  C~86.  Strain  Spectrum  for  Panel  q 
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Figure  C-87.  Strain  Spectrum  for  Panel  q 
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Figure  C-92.  Strain  Spectrum  for  Panel  r 
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Figure  C-97.  Strain  Spectrum  for  Panel  s 
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Figure  C-98.  Strain  Spectrum  for  Panel  s 
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Figure  C-99.  Strain  Spectrum  for  Panel  s 
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The  following  corrections  are  applicable  to  AFWAL-TR-80-3019,  Sonic  Fatigue 
Design  Techniques  for  Advanced  Composite  Aircraft  Structures. 

Delete  pages  177  through  186  (paragraphs  5  and  6)  and  substitute  the 
attached  pages  177  through  188. 
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5.  DESIGN  METHOD 

This  section  contains  a  semiempirical  design  method  for  estimat¬ 
ing  rms  strains  and  natural  frequencies  for  curves  and  flat  cfrp 
stiffened-skin  panels,  subjected  to  random  acoustic  loading.  The  method 
comprises  an  equation  relating  rms  strain  to  panel  configuration 
parameters  and  overall  sound  pressure  levels,  a  procedure  for  estimating 
the  natural  frequencies  of  the  fundamental  fully-fixed  panel  mode,  and  a 
random  fatigue  curve  for  estimating  sonic  fatigue  life. 

The  method  can  be  directly  applied  to  2  and  J  section  stiffeners  and  is 
valid  for  quasi-isotropic  and  most  orthotropic  skin  laminates  typically 
used  for  airplane  skin  panels.  Because  the  response  equation  uses 
overall  sound  pressure  level  as  the  applied  load,  a  technique  for 
utilizing  fan  noise  spectra  from  high  bypass  ratio  jet  engines  is  also 
included.  A  worked  example  is  presented  at, the  end  of  the  section. 

Range  of  Application 

The  use  of  regression  analysis  techniques  in  developing  the  strain 
response  equation  limits  its  application  to  within  the  following  range  of 
panel  configuration  parameters: 

Stiffener  spacing:  4  to  8  inches  (102  to  204  mm) 

Skin  laminate  thickness:  0.Q33  to  0.065  inch  (0.8  to  1.8  mm) 

Strinaer  laminate  thickness:  0.04  to  0.08  inch  (1.0  to  1.0  mm) 

‘  (1.2  times  the  skin  thickness) 

Radius  of  curvature:  flat  down  to  30  inches  (760  mm) 

Aspect  ratio:  1  to  3  (assumed  valid  for  all  values  above  unity) 
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Within  this  parameter  envelope,  the  estimated  strains  showed  an  averaoe 
deviation  of  9  percent  from  measured  values.  A  progressive  deteriorate 
of  this  accuracy  can  be  expected  for  panel  configurations  outside  tne 
given  range.  The  90  percent  confidence  interval  for  estimated  rms 
strains  approximates  a  ±22  percent  variation  in  accuracy. 

CFRP  Skin  Laminates  --  Quasi-isotropic  and  orthotropic  laminates  were 
used  in  the  development  of  the  strain  response  equation  and  the  fatigue 
curve.  Moderate  variations  in  skin  ply  orientation  and  stacking  order  dc 
not  have  a  significant  effect  on  the  accuracy  of  this  design  method. 

The  following  table  gives  the  ply  orientation  and  elastic  modulus  values 
for  the  skin  laminates  used  in  this  programme: 


Table  40.  Skin  Laminate  Plv  Orientations 


ancLila Stic  Modulus  Values 
laminate  density  c  =  0.055  lb/in^ 

(1,522  Kg/m.3) 


PLY  ORIENTATION 

0°  -  Transverse  to  Stringers 
90°  -  Parallel  to  Stringers 

Elastic  Modulus 

(E„) 

Transverse  to 
stringers : 
b  -  direction 

(Ex> 

Parallel  to 
stringers : 
a  -  direction 

106lb/in2 

MN/m2 

1 061 b/i n2 

MN/m2 

(0,  ±4 5)^ 

mm 

51 ,711 

.  3.3 

22,753 

(0,  ±a5,  90),. 

mm 

46,195 

6.7 

46,195 

(02,  ±45}s 

9.7 

66,879 

3.2 

22,063 

(0,  ±45)25 

mm 

51  ,711 

3.3 

22,753 

(0,  90)s 

wa 

on 

m 

64,811 

Stiffener  Design  and  Method  of  Attachment  —  The  response  equation  is 
based  on  data  from  skins  reinforced  with  z  section  stiffeners,  adhesively 
bonded  to  the  skins.  However,  the  equation  can  be  applied  to  J-stiffened 
panels  by  multiplying  calculated  strains  by  0.75.  This  factor  could  also 
be  applied  to  hat  section  stiffeners. 
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RMS  strain  levels  are  not  significantly  affected  by  the  method 
of  skin-stiffener  attachment.  Consequently,  secondary  bonding, 
co-curing,  integral  skin-stiffener  layup  and  rivPted  attachments 
can  be  analyzed.  The  fatigue  curve,  however,  which  is  based  on 
secondarily  bonded  stiffeners,  may  not  be  valid  for  the  other 
methods  of  attachment. 

Stiffener  cross-section  properties  should  provide  for  effective 
panel  edge  restraint.  The  upstanding  webs  of  stringers  should 
be  attached  to  intersecting  frames  in  order  to  provide  continuity 
of  stiffness  and  to  prevent  the  stringers  from  rotating. 


Structural  Damping  --  The  effects  of  damping  have  been  absorbed  into 
the  empirical  factors  in  the  strain  response  equation.  The  equation 
is  valid  for  damping  ratios  in  the  0.017  to  0.03  range.  The  average 
measured  value  was  0.025.  This  range  of  values  may  be  considered 
typical  for  this  class  of  structure.  If  special  damping  factors  need 
to  be  considered,  such  as  the  use  of  damping  treatments,  highly  dam 
resin  systems,  discontinuous  carbon  fibres,  etc.,  then  the  estimate 
rms  strain  should  be  appropriately  factored  down.  In  the  absence  of 
alternative  relevant  data,  multiply  rms  strain  by  the  square  root  of 
the  ratio  (0.025/Actual  Damping  Ratio). 

Nonlinear  Response  Effects  --  Nonlinear  response  effects  have  been 
taken  into  account  to  a  limited  extent,  by  deriving  a  lower  rate  of 
change  in  rms  strains,  with  respect  to  overall  sound  pressure  levels, 
than  is  associated  with  linear  response  behaviour.  These  effects  have 
been  averaged  over  the  data  base  and  do  not  take  into  account  individ¬ 
ual  variations  in  the  degree  of  nonlinear  response  characteristics. 
The  degree  of  nonlinear  response  approximates  to  a  7-dB  increase  in 
sound  pressure  level  resulting  in  a  doubling  of  rms  strain.  Linear 
behaviour  results  in  a  doubling  of  rms  strain  for  a  6-dB  increase  in 
sound  pressure  level. 
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Units  of  Measurement  --  The  equation  for  estimating  rms  strain  uses 
panel  width,  panel  length  and  skin  laminate  thickness  in  non-dimensional 
form.  Consequently,  any  coherent  system  of  units  may  be  used  to  estimate 
rms  strains  for  flat  panels.  The  radius  of  curvature  (R)  is  expressed  in 
inches,  with  an  alternate  expression  provided  for  R  in  millimetres. 

Notation 

The  following  notation  is  used  in  this  section: 

t __r  -  RMS  strain  x  106,  located  at  the  centre  of,  and  normal  to, 

rms 

the  longer  side  of  the  skin  panel 

a  -  length  of  longer  side  of  panel  (frame  or  longeron  spacing) 

b  -  length  (arc  length)  of  shorter  side  of  panel  (stringer 

spacing) 

t  -  skin  laminate  thickness 

R  -  Radius  of  curvature  in  the  b-direction  (inches  or  mm) 

SPL  -  Overall  sound  pressure  level  (dB) 

“  Youngs  modulus  for  the  skin  laminate  material  in  the 

direction  of  the  shorter  side  of  the  panel  (lb/in2  or  N/m2). 
For  laminates  used  in  this  programme,  obtain  values  from 
Table  40. 

p  -  density  of  skin  laminate  material:  (lb/in^)  *  386.4  or 

kg/m^. 

f  -  fundamental  natural  frequency  of  skin  panel  assuming  all 
edges  to  be  fixed.  (Hz) 


180 


V 


velocity  parameter  for  the  skin  laminate  material.  Equals 
1 12 

(E  Ip)  4  200,000  when  E  and  p  are  expressed  in  lb  and  in, 
y  2/2 

and  (E y/p )  4  5,080  when  E^  and  p  are  expressed  in  Kg  and 

m. 


Calculation  of  rms  Strain 

Equation  25  is  used  to  estimate  rms  strain: 


,4/3 


'rms 


=  (|-)  [4  Tanh  (a/b)  -  1] 


10v 


'SPL  -  178 
74 


Tanh 


W) 


(25) 


where  R  is  in  inches. 

When  R  is  in  metres,  the  equation  is  written: 


,4/3 


rms  =  (|)  [4  Tanh  (a/b)  -  1]  10 


/SPL  -  178\ 
~W 


)  Tanh 

10J 


0.43) 


(25a) 


For  large  radii  of  curvature,  (i.e.,  R  greater  than  150  inches  or  greater 
than  4,000  mm)  the  hyperbolic  tangent  of  the  radius  function  is  unity, 
and  the  equation  reduces  to  the  flat  panel  response  equation. 

No  particular  physical  significance  is  attached  to  the  number  (178)  that 
is  subtracted  from  the  sound  pressure  level  in  the  exponent  of  10,  and 
the  equation  is  valid  for  both  positive  and  negative  exponents. 

NOTE  1:  The  equations  were  derived  for  Z  stiffeners.  For  J  stiffeners 

multi  ply  e  by  0.75. 

J  rms 

NOTE  2:  The  equations  are  valid  for  typical  damping  ratios  over  the 

range  0.017  to  0.03.  For  significantly  different  dampinc  ratio 
values  multiply  c  by  the  square  root  of  the  ratio 
(0.025/actual  damping  ratio). 


NOTE  3:  If  a  90S  level  of  confidence  is  required,  then  increase  the 
estimated  rms  strain  by  22%. 
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Calculation  of  Natural  Frequency 


f 


(26) 


where  Vis  defined  under  Notation  and  K  is  obtained  from 

Figure  64  for  given  a/b,  b,  t  and  R  (expressed  in  the  form 
2 

b  /  Rt).  Figure  64  is  expressed  in  both  British  and  S.I.  units. 
Estimation  of  Sonic  Fatigue  Life 

The  estimated  sonic  fatigue  life  is  obtained  by  reading  the  number 

of  cycles  to  failure  (N)  from  Figure  65,  corresponding  to  the  estimated 

rms  strain  (e  ).  The  number  of  cycles  to  failure  is  converted  to  life 
r  ms  ^ 

in  hours  by  the  relationship  (3  g'Q'Qf)»  where  f  is  the  natural  frequency 
calculated  using  Equation  26. 


Worked  Example 

A  six-ply  skin  laminate  having  a  ply  orientation  of  (0,  ±45)s  and 
a  thickness  of  CL033  inch,  has  an  8-inch  stringer  soacing,  with  a  panel 
length  of  12  inches  and  a  radius  of  curvature  of  53  ir.cr.es.  Tr.e  overall 
sound  pressure  level  is  160  dB. 

(i)  For  a  =  12  in,  b  =  8  in,  t  =  0.033  in,  R  =  90  in  and  SPL  =  160, 
equation  25  gives 

4/3  /160-178  \ 

'  Grrax)  r4  Tanh  ¥  -V  1 24 

=  (1,512) (2.62) (0. 178) (0.949) 

=  669  micro-strain 
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Figure  65.  Random  Fatigue  Curve  for  CFRP  Sti f fened-Skin  Panels-RMS  Strain  vs  Cycles  to  Failure 


(ii)  Natural  frequency  is  estimated  from  equation  26: 

f  =  VK  K 
b 

(yp)  1/2 
where  v  =  mmu 

From  Table  40:  E^  =  7.5  x  10^  Ib/ir. 

0.055 
p  =  ISO- 

=  1.148 


From  Figure  64: 


r  b  _  _  8 _  _  n  1  c 

'"or  Rt  "  90  (0.033) 
and  a/b  =  12/8  =  1.5,  then 


K  x  10~6  =  0.6 

k  =  0.6  x  106 

;.  f  =  (1.148)  (0.6  x  106) 
=  355  Hz 


(0.033) 

~~7~ 


( i  i  i )  From  F i gure  65 : 


Number  of  cycles  to  failure  for  crms  =  659  is 
N  =  3.5  x  10^  cycles 

For  a  frequency  of  355  Hz,  the  estimated  sonic  fatigue  life  is 

3.5  x  106 
(3,600) (355 ) 

=  2.74  hours 
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Fan  Noise  from  High  Bypass  Ratio  Jet  Engines 
The  design  method  presented  here  utilizes  the  overall  sound 
pressure  level  as  the  design  load.  While  this  is  adequate  for  the  broa 
band  spectra  that  typify  jet  exhaust  noise,  it  does  not  automatically 
lend  itself  to  designing  for  the  inlet  and  fan  exit  acoustic  spectra  th. 
occur  due  to  fan  noise  on  high  bypass  ratio  engines.  These  spectra  oftt 
have  overall  sound  pressure  levels  that  are  dominated  by  high  acoustic 
spectrum  levels  occurring  at  the  fan  blade  passage  frequency  and  Its 
next  harmonic  (typically  in  the  frequency  range  of  2  to  2.5  KHz  and 
4  to  5  KHz  respectively).  These  frequencies  are  above  the  usual 
frequency  range  of  panel  resonant  response.  This,  including  their 
corresponding  acoustic  spectrum  peaks  in  the  overall  sound  pressure  leve 
to  be  used  for  design  purposes,  may  result  in  overly  conservative  design 
To  deal  with  this  type  of  acoustic  load,  the  applied  overall  sound 
pressure  level  should  not  include  these  two  acoustic  spectrum 
peaks . 

Since  design  acoustic  loads  are  usually  given  from  0  to  10,000  Hz  in  one- 
third  octave  or  one-octave  band  levels,  the  applied  overall  sound 
pressure  level  can  be  obtained  by  summing  these  levels  up  to  10,000  Hz. 
The  blade  passage  tone  and  its  next  harmonic  can  be  eliminated  (if 
predominant)  by  reducing  the  one-third  or  one-octave  band  levels  con¬ 
taining  these  peaks  to  those  levels  contained  in  adjacent  frequency 
bands.  The  one-third  octave  or  one-octave  levels  can  then  be  summed 
in  the  following  way: 


If  the  difference  (in  dB)  between  two 
band  levels  is: 


Then  add  to  the 
larger  level  (dE): 


Example:  If  there  were  four  bands  with  the  following  levels  (d5): 


0 

1 

2 

3 

4 

5 

_ 

6 

7 

8 

9 

10 

>10 

3 

2.5 

2.1 

00 

r—i 

1.5 

CXI 

rH 

1.0 

o 

00 

0.6 

0.5 

0.4 

0 

140,  134, 

\  / 

difference  =  6  dB 
add  140  +  1.0 


=  141dB 


139  and  139, 

\  ./ 


then 


difference  =  0  dB 
then  add  139  +  3 
=  142  dB 


difference  +  1  dB 


add  142  +  2.5 
s  144.5  dB 
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This  144.5  dB  is  the  overall  sound  pressure  level  to  be  used  as 
the  design  load. 

6.  COMPARISONS  WITH  DESIGN  METHODS  FOR  ALUMINUM  PANELS 
CFRP  skin-stringer  structures,  of  the  type  evaluated  in  this 
programme,  are  primarily  in  competition  with  similarly  configured 
aluminum  structures  for  application  on  both  military  and  civil 
airplanes.  Cost/weight  tradeoffs  between  cfrp  and  aluminum 
structures,  having  comparable  sonic  fatigue  resistance,  are  there¬ 
fore  of  interest  to  potential  users  of  the  design  method  in  Section  5. 
A  typical  broadband  sonic  fatigue  design  load  spectrum  was  used 
which  had  an  overall  sound  pressure  level  of  157  dB  and  a  correspond¬ 
ing  acoustic  spectrum  level  of  132  dB/ Hz  in  the  frequency  range  of 
interest.  This  25  dB  difference  between  overall  and  spectrum  levels 
is  compatible  with  the  acoustic  load  in  this  programme. 


The  following  example  problem  was  used  for  comparison: 

7 

Required  life  =  10  cycles  (using  a  50%  confidence  level) 

b  =  8 
a/b  =  2 
C  =  0.02 

Calculate  required  skin  thickness. 

The  following  results  were  obtained: 

Riveted  Aluminiun  Skin-Stringer:  t  =  0.05-in. 

Riveted  Aluminium  Skin-Stringer:  t  *  0.076-in. 

Bonded  cfrp  Skin-Stringer :  t  =  0.037-in. 
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An  8-ply  skin  laminate  having  a  thickness  of  0.044  inch  would  fall 
within  the  90-percent  confidence  interval  for  equation  25. 

Since  the  material  density  of  cfrp  is  approximately  half  that  of 
aluminum,  the  difference  in  weight  between  cfrp  panels  and 
aluminum  panels  having  comparable  sonic  fatigue  resistance, 
approaches  2-1/2  to  1. 
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